Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,814 A→G 11.0% A632A (GCA→GCG)  rsxC → SoxR iron‑sulfur cluster reduction factor component, putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,8140AG11.0% 147.4 / 2.4 74A632A (GCA→GCG) rsxCSoxR iron‑sulfur cluster reduction factor component, putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base G (3/7);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGA  >  NC_000913/1708577‑1709047
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
tCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATTCGAGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCAGCCACAAACGGCGGTCGCAGCGACTACTGCCCGTGCCAAAGCAc                                                                                                                                                                                                                                           >  2:1033369/1‑239 (MQ=17)
       aaaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGc                                                                                                                                                                                                                                    <  1:1323040/238‑1 (MQ=255)
              ccGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCCGTCGCAGCGGCTCTTGCCCTCGCCAAAGCAC‑CCCAGCTGGcaca                                                                                                                                                                                                                             >  2:478294/1‑238 (MQ=14)
               cGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGACCGGGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAACCGCCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAG‑GCCAGCGTGAACAGc                                                                                                                                                                                                                           >  2:701124/1‑239 (MQ=14)
                  tGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGACCAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCaa                                                                                                                                                                                                                         <  1:208952/238‑1 (MQ=255)
                  tGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCaa                                                                                                                                                                                                                         <  1:240610/238‑1 (MQ=255)
                        agcTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCt                                                                                                                                                                                                                   <  1:1078485/238‑1 (MQ=255)
                           tatgggcgctgacaaaaagcgcAAACGAGAACAGCAACAGGCTCATGCCGAGCCAGAACACCAGCTCGATCCGCGCAACGCCGGCGTGGAAGCCGCGAGTGCGCGAGCCAAAGCGAGCAAACGGGAACAG‑CAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATg                                                                                                                                                                                                               <  2:1254558/221‑1 (MQ=14)
                            aTCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGc                                                                                                                                                                                                              <  1:28777/239‑1 (MQ=255)
                             tcgacgtcgccgaaagccgcAAGCGGGAAGAGCAACAGGTAAATAACGAGCCAGAACACCAGGACGACCGGCGAAAAGCGGCCGCCGAAGCCGCTATCGCCCGTGCAAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGc                                                                                                                                                                                                              <  2:256993/221‑1 (MQ=14)
                                 cagagcaaaaGGGGGCAAGCGGAAACAGCAACAGGAAAATGAGGACCCAGAACAACAGGCAGATGCGCGAAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGc                                                                                                                                                                                                         <  2:1202239/232‑1 (MQ=14)
                                   tggcgaaaGCGCTCAAGCGGGAAATGAGATAGGTAAACGCGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCAAAAGCGGAACCAGAAGACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCa                                                                                                                                                                                                       <  2:182660/234‑1 (MQ=14)
                                   gcgcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCCAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCCCCGTTGACGCGGCTATGGACCGTGCGACGGCCC‑GCACACTGGAGCAGCAACAGGCTacggcggcgccc                                                                                                                                                                                                       >  2:1139505/1‑228 (MQ=14)
                                    tgcCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCCCGCAAACGGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTAGTAGCCGATATTCCCCGTGCCAAAGCAC‑ACAAACTGCAGCAGCAACACGCTAATCACGAGccca                                                                                                                                                                                                      >  2:1022766/2‑237 (MQ=11)
                                    cgcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCAC‑GCAAACTGGAACAGCAACAGGCTAATGCCGACCCAg                                                                                                                                                                                                      >  1:32654/1‑239 (MQ=255)
                                    cgcCAAAGCGCGCAAGCTCGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAACCGCCGTCGACGCCGCTATTCCCCGTGCCCAAACGCGCAAGTTGGACCAG‑CAACAGGCTAATGCGGAACCCGAACAACAGGTCGATCCACGACAAGCCGCCGTTGAGGCAGCCATTGCCTGTGCCAAATCACTCCAAGC‑GGACCAGCAACAgcttcctcatgtttcag                                                                                                                                                                                                      >  2:691633/1‑223 (MQ=11)
                                        aaaGCCGGCAAGAAGGAACAGCAACAGGAACATGAGAAAGCAGAACACCAGGCCGAACCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaaca                                                                                                                                                                                                  <  2:763711/239‑1 (MQ=14)
                                              cgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCAC‑ACAAGATGGAACAGCAACCGGCTACTGCCGACGCAAACCAACAga                                                                                                                                                                                             >  2:1370547/1‑237 (MQ=14)
                                                             cAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAACCAGAAGAACAGGTCGATCCGCCCAAAg                                                                                                                                                                              >  1:266038/1‑238 (MQ=255)
                                                             cAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCCCGCAAACTGGACCAG‑CAGCAGGCTAATGCGCAGCCACAAAAACAGGTCGGTCCGCGCAAACCCGCCGTCGAGGCGGCTATTCCCCGTGCCAAGGCAC‑GCAAGCCGAAAAAGCAACGGTCGAATACCGCGCCCGaaaacacggctaccccacccaaac                                                                                                                                                                              >  2:310571/1‑216 (MQ=2)
                                                                  ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTCGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCCC‑GCAAGCTGGACCAGCAACAGGCTAATGCCGCGCCAGAACAACAGGTCGATCCGCGCAAAGccgccg                                                                                                                                                                        >  1:576212/1‑239 (MQ=255)
                                                                       tagCGGGACCAGAACACCAGGACGCTCGGCGGAAAGCCGCCGGCGAAGCAGCTATCGCCCGGGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGaa                                                                                                                                                                   <  2:1380024/237‑1 (MQ=14)
                                                                        tGCGGACCCGGACAAACGGGTGGATCCGCAGAAAGCCGCCAGTGAAGCCGCAAGTGCCCGGGCCAAAGGGCGGAAGCTGGAACAG‑CAACAGGCGAAGGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAg                                                                                                                                                                  <  2:1405592/239‑1 (MQ=14)
                                                                                   aacaacAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGc                                                                                                                                                        >  1:1214964/1‑238 (MQ=255)
                                                                                   aacaacAGGTCGATCCGCGCAAAGACGCCGTCGAAGCCGCTATTGCTCGTGCCAAAGCGCGCAAGCAGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGACGACGCCGCTATTGCCCGGGCCAAAGCGG‑CAAAGAGGGAAGCGCAACAGGCTGATACAGCGCCGCAACAGCAGACCGAGCCGCGGATAGccgccggcacatacacactcgc                                                                                                                                                        >  2:995044/1‑222 (MQ=9)
                                                                                           gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCGCGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGGACGATCCGCGCAAAGCAGCCGTCGCTGCCGCTATTGCCCgcaccca                                                                                                                                               >  1:1006619/1‑233 (MQ=14)
                                                                                             ggaccgaccaaatactgccgccGAAGCCGCTATGGCCCGAGCAAGAGCGGGAAAGCTGGAACAG‑CAACAGGATAATGCGGAACAAGAAGAACAGGCCGATCCGCGCAAAGCCGCCGGCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaaa                                                                                                                                              <  2:210456/222‑1 (MQ=14)
                                                                                             cGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCaaa                                                                                                                                              >  1:845688/1‑238 (MQ=255)
                                                                                                                 ctgAAGCCGCGATTGCGCGTGCCCAAGCGCTCCGGCT‑GTACAGACGACGGGCTATTGCGGAACCAGCAGACGTGGCCGATCCGCGCAAAGCTGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCaa                                                                                                                         <  2:819935/237‑1 (MQ=14)
                                                                                                                   gAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTCGACGCAGCTATTGCCCGCGCTCAAGCCA‑GCAACCAGGACCAGCAAAAGGCTAATGCGGACCCCGAAGACCAGGGCGATCCGCGCAAAGCCGCCGTAGCCGCGAATATTGCCCGAGCCAAAAGACGCAAGCggcgccagcacac                                                                                                                       >  2:823621/1‑228 (MQ=2)
                                                                                                                        cGCTATTGCCCGCGCAAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCt                                                                                                                   <  1:770740/238‑1 (MQ=255)
                                                                                                                                     gCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCaga                                                                                                     >  1:1117526/1‑239 (MQ=255)
                                                                                                                                     gCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGACCGAGCCAAAGCGCGACAGCCGGAACAGCACACGGCCACATCGGAGCCaga                                                                                                     >  2:624478/1‑239 (MQ=14)
                                                                                                                                        aaaGGGAGCAAGCAGGAACAG‑CAACAGGCAAATAAGGAGCCAGAACAAGAGGCGGATCCGCGAAAAGCCGCCGTCGAAGCCGCTATGGCCCGTGCCAAAGCGC‑GCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                                                                                                  <  2:1292618/239‑1 (MQ=11)
                                                                                                                                        aaaGCGGGAAAGCTGAAAAAG‑CAACAGGAAAATACGGACGCAGAACAACAGGTCGATCCGGGAAAAACCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                                                                                                  <  2:1271610/239‑1 (MQ=14)
                                                                                                                                         aaGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                                                                                                  >  2:692520/1‑238 (MQ=255)
                                                                                                                                          aGCGCGCAAGCTGGAACAG‑CAACATGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAAc                                                                                                >  1:304023/1‑239 (MQ=255)
                                                                                                                                              cgcAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAACTGGAACAGCAACAGGCTAATGCCGAGCCCGAACAACAGGTCGATCCGCGCAAAACCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCAGCAAACCGGACCAGCACCCGGCTAACACGGAgaccgcacaaccagt                                                                                            >  2:1230272/1‑225 (MQ=14)
                                                                                                                                              cgcAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCAC‑GCAAGCTGGACCAGCAACCGGCTAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGCGCCAAAGACCGCAAACGGGAACAGCAACCGGCTATTGAGACACCAGACGAACAggc                                                                                            >  2:482862/1‑238 (MQ=14)
                                                                                                                                               gcAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGa                                                                                                                                                                                                           <  2:188937/127‑1 (MQ=38)
                                                                                                                                                cAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                                                                                                                                                                                        <  1:1349681/145‑1 (MQ=255)
                                                                                                                                                cAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                                                                                                                                                                                        >  2:1349681/1‑145 (MQ=255)
                                                                                                                                                              cAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGc                                                                             >  1:1020000/1‑239 (MQ=255)
                                                                                                                                                              cAACAGGCTAATGCCCAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAGCTGGAACAGCAACAGTCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCTCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaaga                                                                                                  <  1:122259/218‑1 (MQ=11)
                                                                                                                                                               aaCAGGCTACTGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCt                                                                            <  1:1057591/239‑1 (MQ=255)
                                                                                                                                                                   ggCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTccc                                                                         >  1:1254558/1‑238 (MQ=255)
                                                                                                                                                                       agtacagagcgaaaacaccAGGGCGACCCGCGAAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCAAAAGGAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGa                                                                    <  2:1228943/221‑1 (MQ=14)
                                                                                                                                                                       actGCGGAAGCAGAGGAACAGGCCGATCCGCGCAAAGCGGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGa                                                                    <  1:1317985/237‑1 (MQ=255)
                                                                                                                                                                        aTGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAAc                                                                                                                                                                                                                              <  1:386348/84‑1 (MQ=38)
                                                                                                                                                                        aTGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAGCAGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATGCTCGCAAAGCTGCCGACGAAGCGGTTCTTGCCAGAGCCAAAGCGCGCAGACGGGCACAGGCACCGGCTATTGAGGAGCCCGAAAAAAAGGTAGAtactgcctcttatccatctct                                                                    >  2:122259/1‑218 (MQ=2)
                                                                                                                                                                         tGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGaa                                                                   >  1:351467/1‑238 (MQ=255)
                                                                                                                                                                         tGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGAAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAg                                                                  <  1:758100/239‑1 (MQ=255)
                                                                                                                                                                              aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCTGCGCAAAGCCGCCGTCGAAGCCGCt                                                             >  1:718555/1‑239 (MQ=255)
                                                                                                                                                                                   gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGc                                                        >  1:932648/1‑239 (MQ=255)
                                                                                                                                                                                          aGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAAAACAGGTCGATCCGCGCAAAGCGGCCGTCGAAGCCGCTATTTCCCGAGCCAAAGAGCGCAAAGCGGAACAGCAACCGGCTAATGCGGAGCCAGAAAAACAGGTAGATCCGCCCAAAGCTCACATCGACGCGGCTCCTGCAagccc                                                  >  2:960792/1‑236 (MQ=14)
                                                                                                                                                                                           ggTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACGGGCTAAAGCCGAGCCAGAAAAACAGGTTGATCCGCGCAAAGCCGCCGTCAACGCGGCCATTGCACGCGcc                                                 >  2:1151738/1‑238 (MQ=14)
                                                                                                                                                                                               gcgccgcaaaagcgcgcctgcGAAGACGGTATTACCCGGCACAAAGCAG‑CCAAGCTGGAACACGAACAGGCTAATGCCGAGCGAGAACAAGAGGTCGACCCCCGCAAAGCCGCCGTCGAAGCCGCTATCACACGAGACAAAGCGCGCAAACGGGAACACCACGCGGCTAATGCGGAGCCAGAAGACGAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGc                                            <  2:1303214/219‑1 (MQ=255)
                                                                                                                                                                                                 cccacgcAAAGCTGCCGTCGACGCCGCTATTGGCCGTGCAAAAGAAC‑GACAGCTGGGACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGACGAACCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCa                                           <  2:1020000/234‑1 (MQ=17)
                                                                                                                                                                                                  ccGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCACAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACg                                         <  1:73520/239‑1 (MQ=255)
                                                                                                                                                                                                        aaaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCAAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGc                                    <  1:893171/238‑1 (MQ=255)
                                                                                                                                                                                                                   cGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCaa                         <  1:1014432/238‑1 (MQ=255)
                                                                                                                                                                                                                           ttatcgtccGCGCCAAAGCGC‑GCAAGCTGGAACATCAACATGCTAAGGCGGAACCAGAACAGCAGGTCGAGACGCGCAAAGCCGCCGTGGAAGCCGCTATTGACCGTGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCTTCGa                                                                    <  2:102157/180‑1 (MQ=255)
                                                                                                                                                                                                                           cTATTGCCCGGGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACg                                         <  2:592291/214‑1 (MQ=255)
                                                                                                                                                                                                                                gCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGt           <  1:1033369/239‑1 (MQ=255)
                                                                                                                                                                                                                                gCCCGAGCCAAAGCGC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGCGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGt           <  1:701124/239‑1 (MQ=255)
                                                                                                                                                                                                                                 cccGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGt           <  1:692520/238‑1 (MQ=255)
                                                                                                                                                                                                                                 cccGTGCCAAAGCAC‑GCAAGCTGGAACACCAACAGGCTAATGCCGAGCCAGAACAACAGGACGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGt           <  1:392149/238‑1 (MQ=255)
                                                                                                                                                                                                                                   cGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACACGCCAAAGCACGCAAGCTGGAACAGCAACAGGCAAAAGCGGTAca        >  2:1404639/1‑238 (MQ=255)
                                                                                                                                                                                                                                     tGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCACGCTGGAACAGCAACAGCCTAATACGGTACCag      >  1:1019851/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGaa                               >  1:592291/1‑214 (MQ=255)
                                                                                                                                                                                                                                     cgCCAAAGCGC‑GCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTAtt                                                          >  1:102157/2‑187 (MQ=255)
                                                                                                                                                                                                                                      gCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCaga     >  1:490427/1‑239 (MQ=255)
                                                                                                                                                                                                                                         aaaGCGC‑GCAAGCTGGAACACCAACAGGCTAATGCGGACCCAGAACAGCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCagaaga  <  1:570341/239‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaGAAC‑GAAAGCTGGAACAGCAGCGGGTTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCAGTATTGCACGTGCCAAAGAGCGCAAGCGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCCGCTATTGCACGCGCAAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCagaaga  <  2:266038/239‑1 (MQ=14)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
TCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAG‑CAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGA  >  NC_000913/1708577‑1709047

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: ‑