Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,288,561 (C)7→6 pseudogene (354/2525 nt) yejO ← pseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,288,5550C.98.7% 306.0 / ‑2.7 76pseudogene (360/2525 nt)yejOpseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (38/37);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288434‑2288688
                                                                                                                         |                                                                                                                                     
aaTACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAAc                                                                                                                     >  2:394596/1‑139 (MQ=255)
 aTACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACg                                                                                                                    >  2:1433254/1‑139 (MQ=255)
          tGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGcccccca                                                                                                                                 >  2:489865/1‑117 (MQ=255)
          tGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGcccccca                                                                                                                                 <  1:489865/118‑2 (MQ=255)
          tGGCCCGCATGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCGAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           <  1:985606/139‑1 (MQ=255)
           ggCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  1:400201/139‑1 (MQ=255)
            gCCCGCGTGGTTGTTTTAATAGCACCGCTTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTcca                                                                                                                      <  2:992948/126‑3 (MQ=255)
            gCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCaa                                                                                                                      >  1:992948/1‑126 (MQ=255)
               cgcgTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcccgcc                                                                                                      >  1:1539046/1‑139 (MQ=255)
                  gTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                                                                   <  2:1539046/139‑1 (MQ=255)
                        gttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGt                                                                                             >  1:997542/1‑139 (MQ=255)
                         ttttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                                                            <  2:997542/139‑1 (MQ=255)
                         ttttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                                                            <  1:1433254/139‑1 (MQ=255)
                                  cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTc                                                                                   >  1:19672/1‑139 (MQ=255)
                                   aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGTGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCg                                                                                  >  1:279292/1‑139 (MQ=255)
                                   aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCg                                                                                  >  1:1254710/1‑139 (MQ=255)
                                    ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              <  2:1390557/126‑1 (MQ=255)
                                    ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              >  1:1390557/1‑126 (MQ=255)
                                      gCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  2:1030117/139‑1 (MQ=255)
                                      gCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:578842/139‑1 (MQ=255)
                                       ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgc                                                                              >  1:1079609/1‑139 (MQ=255)
                                       ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgc                                                                              >  2:500375/1‑139 (MQ=255)
                                           aagtCTGATCTACCCCAAATGCCAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCa                                                                          <  2:1322617/136‑1 (MQ=255)
                                                     tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTa                                                                >  1:1278687/1‑139 (MQ=255)
                                                     tACCGCAAATGCAAATACCCCCTCTTTAACCTACAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTa                                                                >  2:1186335/1‑139 (MQ=255)
                                                      aCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGTGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  2:279292/139‑1 (MQ=255)
                                                      aCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  2:1418323/139‑1 (MQ=255)
                                                      aCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  1:94292/139‑1 (MQ=255)
                                                       ccGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGc                                                              >  1:1549665/1‑139 (MQ=255)
                                                        cGCAAATGCAAATCCCCCCTCTTTAACCTCCGTTAGCCCTCCTTTTTCGATACGGACAAAGAC‑GCCCCCCATTCTTTTTCAACGATTAGGCGCCGCCCTTAAATCTTGTTACATTTCGCTGCGCCATAAACCTGTGGct                                                             >  1:1322617/1‑138 (MQ=255)
                                                          caaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                           >  2:213023/1‑139 (MQ=255)
                                                             atgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCa                                                        >  1:210967/1‑139 (MQ=255)
                                                                  aatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattatt                                                   >  1:956745/1‑139 (MQ=255)
                                                                     ccccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAg                                                <  1:394596/139‑1 (MQ=255)
                                                                       ccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCg                                              <  2:624115/139‑1 (MQ=255)
                                                                         ccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                             <  1:197345/122‑1 (MQ=255)
                                                                         ccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                             >  2:197345/1‑122 (MQ=255)
                                                                         ccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTa                                            >  2:856395/1‑139 (MQ=255)
                                                                         ccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTa                                            >  2:789804/1‑139 (MQ=255)
                                                                         ccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTa                                            >  1:1635672/1‑139 (MQ=255)
                                                                               tAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTgg                                      >  2:355916/1‑139 (MQ=255)
                                                                                   cTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTg                                  <  2:1079609/139‑1 (MQ=255)
                                                                                      cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                                 >  1:833046/1‑73 (MQ=255)
                                                                                      cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                                 <  2:833046/73‑1 (MQ=255)
                                                                                      cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcg                               >  1:895202/1‑139 (MQ=255)
                                                                                       agtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                           >  1:966540/1‑110 (MQ=255)
                                                                                       agtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcgc                              >  1:1493868/1‑139 (MQ=255)
                                                                                       agtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcgc                              >  1:199088/1‑139 (MQ=255)
                                                                                       agtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATATTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                           <  2:966540/110‑1 (MQ=255)
                                                                                        gtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcgcg                             <  1:1416968/139‑1 (MQ=255)
                                                                                         tagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcccgc                                                                                                       >  1:781449/1‑64 (MQ=255)
                                                                                         tagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcccgc                                                                                                       <  2:781449/64‑1 (MQ=255)
                                                                                           gtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCa                          <  1:1201192/139‑1 (MQ=255)
                                                                                           gtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCa                          <  1:690038/139‑1 (MQ=255)
                                                                                           gtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCa                          <  2:19672/139‑1 (MQ=255)
                                                                                            tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATc                                                                                               <  1:1326187/69‑1 (MQ=255)
                                                                                            tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATc                                                                                               >  2:1326187/1‑69 (MQ=255)
                                                                                             cctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGc                                               >  2:1591692/1‑116 (MQ=255)
                                                                                             cctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGc                                               <  1:1591692/116‑1 (MQ=255)
                                                                                              ctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               >  1:1622748/1‑99 (MQ=255)
                                                                                              ctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  2:1622748/99‑1 (MQ=255)
                                                                                              ctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCAccc                       >  1:125249/1‑139 (MQ=255)
                                                                                              ctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCAccc                       >  2:1281755/1‑139 (MQ=255)
                                                                                                cccTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                                                            >  1:386604/1‑68 (MQ=255)
                                                                                                cccTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                                                            <  2:386604/68‑1 (MQ=255)
                                                                                                     ttCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  2:1111969/1‑76 (MQ=255)
                                                                                                     ttCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:1111969/76‑1 (MQ=255)
                                                                                                      tCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAAc               <  1:856395/139‑1 (MQ=255)
                                                                                                            gCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCa                          >  1:1317327/1‑122 (MQ=255)
                                                                                                            gCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCa                          <  2:1317327/122‑1 (MQ=255)
                                                                                                              gACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGAt       <  1:789804/139‑1 (MQ=255)
                                                                                                               aCAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           >  1:121299/1‑38 (MQ=255)
                                                                                                               aCAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           <  2:121299/38‑1 (MQ=255)
                                                                                                               aCAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                           <  1:1535930/86‑1 (MQ=255)
                                                                                                               aCAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                           >  2:1535930/1‑86 (MQ=255)
                                                                                                                   aGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  2:895202/139‑1 (MQ=255)
                                                                                                                   aGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  2:125249/139‑1 (MQ=255)
                                                                                                                   aGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGAGTAGTCTGGTTTGCCCCGCCACCCATTTGAACAATTTGATAACCg  <  2:956745/139‑1 (MQ=255)
                                                                                                                         |                                                                                                                                     
AATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288434‑2288688

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑