Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,473,746 T→G 21.5% pseudogene (3230/3497 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,473,7460TG21.5% 106.3 / 3.2 51pseudogene (3230/3497 nt)ydbApseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (17/20);  new base G (5/8);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATACCGACGTGAAACAGCAGTACCTGGAGTTCCGCAGCGAAGGGGCGAAGACTTTCGAACCGAGCGAAGGACTGAAGGTTACGCCATATGCGGGTGTAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAACTGGAAAGAGGATGGCATCAGCGATA  >  NC_000913/1473520‑1473978
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         
cTGGTACCGACGTGAAACAGCAGTACCTGGAGTTCCGCAGCGAAGGGGCGAAGACTTTCGAACCGAGCGAAGGACTGAAGGTTACGCCATATGCGGGTGTAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTg                                                                                                                                                                                                                                <  2:442481/238‑1 (MQ=255)
          cGTGAAACAGCAGTACCTGGAGTTCCGCAGCGAAGGGGCGAAGACTTTCGAACCGAGCGAAGGACTGAAGGTTACGCCATATGCGGGTGTAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGGGGACAGCATCGTCGGGCTGAAACTTGACTACGCAGGGCAAGACGGCCGGagggccagc                                                                                                                                                                                                                     >  2:736724/1‑238 (MQ=255)
          cGTGAAACAGCAGTACCTGGAGTTCCGCAGCGAAGGGGCGAAGACTTTCGAACCGAGCGAAGGACTGAAGGTTACGCCATATGCGGGTGTAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGGAACTGGCCTCACCCGGTACAGACGGCTGGagcgctagc                                                                                                                                                                                                                     >  2:473572/1‑237 (MQ=255)
                                   cGCAGCGAAGGGGCGAAGACTTTCGAACCGAGCGAAGGACTGAAGGTTACGCCATATGCGGGTGTAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGGCTACGCAGGGAAAGACGGCTGGAGGGCTAGCGCTACCCTGGAAGGCGGGCCGaa                                                                                                                                                                                              >  2:206996/1‑236 (MQ=255)
                                   cGCAGCGAAGGGGCGAAGACTTTCGAACCGAGCGAAGGACTGAAGGTTACGCCATATGCGGGTGTAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGGCGGGCTGAAACTGGACTACGCAGGGAAAAGAGGCCGGGGCGCCAGCGCTAGCCCGGAAGGCGGGCCCAAccc                                                                                                                                                                                           >  2:978076/1‑238 (MQ=255)
                                          aaGGGGCGAAGACTTTCCAACCGAGCGAAGGACTGAAGGTTACGCCATATGGGGGTGTAAAACTGCGTCACTCACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTa                                                                                                                                                                                     <  2:1179251/238‑1 (MQ=255)
                                                        ttggAACCGGGGGAAGGACTGAAGGTTACGCCAGTGGCGGGTGTTAAACTGCGTCACACACTTGAAGGTGGCTATCAGGAGCGCAATGGCGGAGACTTGAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCg                                                                                                                                                                      <  1:528563/236‑1 (MQ=255)
                                                             cccGGGCGAAGGACTGAAGGTTACGCCATATGCGGGGTTAAAACCGCGTCAACCACTGGAAGGTGGCTAGCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACgg                                                                                                                                                                 <  1:206996/238‑1 (MQ=255)
                                                                      aGGACTGAAGGTTACGCCATATGCGGGTGTAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTgg                                                                                                                                                        <  2:1095165/239‑1 (MQ=255)
                                                                         aCTGAAGGTTACGCCATATGCGGGTGTAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGgca                                                                                                                                                      <  2:1116418/238‑1 (MQ=255)
                                                                           tGAAGGTTACGCCATATGCGGGTGTAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGgcaggc                                                                                                                                                   >  2:1128103/1‑239 (MQ=255)
                                                                                               ggTGTAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGCACCTGAGCTACGCGACGAGCCAGCGGACGGCAAACCTGGCAGGGGCAGGCAGTCAGCCCTTTaa                                                                                                                               >  2:172163/1‑239 (MQ=255)
                                                                                                   tAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGGAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCGACGCGAAGAGCCAACGTACGGCAAGCCGTGCAGGCGCAGGCAGGCAGCACTTTAACGTc                                                                                                                           >  2:1097622/1‑239 (MQ=255)
                                                                                                                      ggAAGTGGGCTATCAGGAGCCAAATGCCGGAGTCTTTAACCTGAATATGAACGGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGggcg                                                                                                        <  2:586840/239‑1 (MQ=255)
                                                                                                                          ggTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCCGGTAAAGACGGCTGGGGCGCTAAGCGTACGCTGGAAGGCGGGCCCAACCCGAGCTCACCGAAGAGCCACCCGACGGCAAGCCAGGCAGGGGCAGGCGGCCAGCACATTAAAGTCCATggaggtaaaaaagggggggg                                                                                                    >  1:287455/1‑220 (MQ=255)
                                                                                                                           gTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGGACGGCAAGCCTGGCAGGGGCAGGCAGTCAGCCCCTTAAAGTCGGTGACGGGCAGAAGggcggcggc                                                                                                   >  2:1254202/1‑239 (MQ=255)
                                                                                                                              gCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCCGGCCTTTAACGTCGGTGACGGTCCGAAGGGCGGCGGCATc                                                                                                >  2:598289/1‑239 (MQ=255)
                                                                                                                               tgttgaaaagaggcgagCCGAAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGTACAGCATCGTCGGGCCGAAACTGGACTTCGCAGGGAAAGACGGCTGGAGCGCGAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCATTCAGCACTGTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCa                                                                                               <  2:797855/223‑1 (MQ=255)
                                                                                                                               cTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCCTcc                                                                                               >  2:1047462/1‑238 (MQ=255)
                                                                                                                                      gAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGGAAGCATCGTCGGGCTGAAACTGGACTACGCCGGGAAAGAAGGCTGGAGCGCTAGCGCTACGCTGGGAGGCGGGCCGAAACTGAGCTAAGCGGAGAGCCAGCGGTCGGCAAGCCCGGCCGGGCCCGGCGGGCACCACTTTAACGTGGCGGACGGGCCGAggggggggggagccaaaaaaca                                                                                        >  1:1285445/1‑227 (MQ=255)
                                                                                                                                               ggccgaatcTTTAACCTGAATTGGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGCCGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGt                                                                               <  2:1303625/231‑1 (MQ=255)
                                                                                                                                                    aggttttAACCCGAAAGAGAACCGTGGCAGGGAAACGGCGGTGGACAGCCTCGGCGGGCTGAACCTGGACTACGCAGGTAAAGGAGGCTGGAGGGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGc                                                                           <  1:1097622/234‑1 (MQ=255)
                                                                                                                                                     gaCTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCg                                                                          <  2:287455/238‑1 (MQ=255)
                                                                                                                                                     agctttaacgtaaagtgtgaCAGTGGGGGGGAAAGGGGGGGGGAACGCATGGTCGGGCTGAAACTGGAATACGCAGGGAAAGGCGGCTGGAAGGCCAGCGCTACGCGGGAAGGCGGGCCGAACCTGAGCAACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTAAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCg                                                                          <  1:863303/219‑1 (MQ=255)
                                                                                                                                                                 aaTAAGAACAGGGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCg                                                                          <  2:447382/226‑1 (MQ=255)
                                                                                                                                                                           gTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGGAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGGTGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAACGGTCGGCGTGAAGTACAGCGGCAAAGAAAGt                                                   >  2:696030/1‑239 (MQ=255)
                                                                                                                                                                           gTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGGGGCATCAAAAACCTGACAAGGGTCGGCGGGAAGGAAAGCAGCAAAGAAAGt                                                   >  2:435969/1‑239 (MQ=255)
                                                                                                                                                                            tGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAAAAGCCTGACAAACGTCGGGGTGAAGGACCAGAGCAAAGAAAGt                                                   >  2:670789/1‑238 (MQ=255)
                                                                                                                                                                                ggCGAAACCGGGGGGGGCAGCATGCTGGGGCTGAAACTGGACTCCGCAGGGAAAGACGTGTGGAGGCATAGCGCTACGCTGGAAGGCGGGGCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCt                                              <  1:1254202/238‑1 (MQ=255)
                                                                                                                                                                                 gggAAACGGGGGGGGACAGGAGCGTCGGGGGGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGGGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCt                                              <  1:1128103/236‑1 (MQ=255)
                                                                                                                                                                                 gCGAAACGGCGGTGGACAGCAGCGGCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCt                                              <  2:1264936/238‑1 (MQ=255)
                                                                                                                                                                                        gggggTGGACAGGAGTGGCGGGGTGAAATGTGACTACGGAGGGAAAGGCGGCTGGGGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTgg                                      <  1:172163/236‑1 (MQ=255)
                                                                                                                                                                                          gggtggacctcggggggggtgtGAAACGGGACTCGGAAGGAAAAGACGGCTGGACGGCTAGCGCGAGGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGGGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTgg                                      <  1:940123/216‑1 (MQ=255)
                                                                                                                                                                                           ggTGGACAGCATCTTGGGGCTGAAAGTGGACTACGCCGGGAAAGACGGCTGGGGCACTAGGGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGAt                                    <  1:736724/238‑1 (MQ=255)
                                                                                                                                                                                              ggACAGCAGCGTCGGGCTGAAACTGGGCTGCGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCAATTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCg                                 <  2:1225643/238‑1 (MQ=255)
                                                                                                                                                                                              ggACAGCAGCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCg                                 <  2:1138933/238‑1 (MQ=255)
                                                                                                                                                                                                  cgCATCGTTGGGGTGAAACGGGTCTGCGCAGGTAAAGGCGGCTGGAGTGCCAGCCACACGGTGGAAGGCGGGCCGAACCGGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACaa                            <  1:598289/238‑1 (MQ=255)
                                                                                                                                                                                                    cATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACa                             >  2:452482/1‑236 (MQ=255)
                                                                                                                                                                                                       aaagggacagAAACCGGACTACGCAGGTAAAGCCGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAACTGGa                       <  2:1118311/230‑1 (MQ=255)
                                                                                                                                                                                                        ggggggTTGAAACTGGACTACGCAGGGAAAAGCGGCGGGAGCGCGAGCGCTCCGCGGGAAGGCGGGCCGAACCTGAGCGACGCGGAGAGCCAGCGTACGGCAAGCCTGGGAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAACTGGaa                      <  1:696030/236‑1 (MQ=255)
                                                                                                                                                                                                         tcgGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCCCTGAATCTGGATGCGTACAACTGGaaa                     >  2:473072/1‑239 (MQ=255)
                                                                                                                                                                                                                  aaCTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAACTGGAAAGAGGATGGc            <  2:1018242/239‑1 (MQ=255)
                                                                                                                                                                                                                  aaCGGGACCACGGAGGGAAAGGCGGCTGGAGCGCTAGCGGTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCg                                 <  1:719048/218‑1 (MQ=255)
                                                                                                                                                                                                                     tGGACTTCGCAGGGAAAGACGGGTGGGGCGCTAGCGCTACGAGGGAAGGCGGGCCGAACCTGAGGTACGCGAAGAGCCAGGGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGGCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAACTGGAAAGAGGATGGCATc         <  1:939863/239‑1 (MQ=255)
                                                                                                                                                                                                                     tGGACTCCGAAGGTAAGGCGGGCTGGGCGGCTGGCGCTGCGCTGGAAGGCGGGGCGAACCTGAGCCACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGGGCAGTCAGTCAGCACTTTAACGTCGATGGCGGTCAGAAGGGCGGCTGCATCAATAGCCTGACAAGCGGCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAACTGGAAAGAGGATGGCATc         <  1:473572/239‑1 (MQ=255)
                                                                                                                                                                                                                     tGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAACTGGAAAGAGGATGGCATc         <  2:1285445/239‑1 (MQ=255)
                                                                                                                                                                                                                       gACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGGACAGCAGCAACGAAAGTTCGCTGAATCTGGATGCGTACAACTGGAAAGAGGATGGCATCa        >  2:230922/1‑238 (MQ=255)
                                                                                                                                                                                                                            cGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAACTGGa                       >  2:719048/1‑218 (MQ=255)
                                                                                                                                                                                                                            cGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGAGGCGTACAACTGGAAAGAGGATGGCATCAGCGATa  >  2:433615/1‑239 (MQ=255)
                                                                                                                                                                                                                            cGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGGGAAGTACAGCCGCAAGGAAAGTTCCCGGAATCTGGATGCGTACAACGGGAAAAAGGATGGcctgtctcttt  >  1:1118311/1‑229 (MQ=255)
                                                                                                                                                                                                                            cGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAACCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAACTGGAAAGAGGAGGGCATCAGCGATa  >  2:348314/1‑239 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         
CTGATACCGACGTGAAACAGCAGTACCTGGAGTTCCGCAGCGAAGGGGCGAAGACTTTCGAACCGAGCGAAGGACTGAAGGTTACGCCATATGCGGGTGTAAAACTGCGTCACACACTGGAAGGTGGCTATCAGGAGCGCAATGCCGGAGACTTTAACCTGAATATGAACAGTGGCAGCGAAACGGCGGTGGACAGCATCGTCGGGCTGAAACTGGACTACGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGTCAGAAGGGCGGCGGCATCAATAGCCTGACAAGCGTCGGCGTGAAGTACAGCAGCAAAGAAAGTTCGCTGAATCTGGATGCGTACAACTGGAAAGAGGATGGCATCAGCGATA  >  NC_000913/1473520‑1473978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: