Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,844 C→G 10.8% A642A (GCC→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,8440CG10.8% 171.9 / 4.6 63A642A (GCC→GCG) rsxCSoxR iron‑sulfur cluster reduction factor component, putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (32/23);  new base G (3/5);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAG  >  NC_000913/1708559‑1709120
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                      
cgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGa                                                                                                                                                                                                                                                                                      <  2:582325/288‑1 (MQ=255)
     aagcggtgagtGATCGACGAATAACTGACGTTGACGCAGCTATCGCCGGGGCCAACGCG‑GCGAGTCTGGACCAGAGACAGGCTAATGCGTAAGCGGAACACCTGGCCGATCCGCGTAAAGCTGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAGGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAg                                                                                                                                                                                                                                                                                 <  1:349672/278‑1 (MQ=17)
      aaCCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                                                                                                                                                                                                                                                                               <  2:412003/289‑1 (MQ=255)
        ccGTTGCTGATCCACGTAAAACTGACGTTGAAGCAGCTATCGACCGCGCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaaca                                                                                                                                                                                                                                                                             <  2:766247/289‑1 (MQ=255)
             gCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAgg                                                                                                                                                                                                                                                                        >  2:853410/1‑289 (MQ=255)
              aagaacgaggTAAAAGTGCCGTGGATGCGGCTATTGACGAGGCCACAGCGCGCAAG‑CCGGCACAGCCACAGGCTAAGGAGGAACCAGAAGAGAAGGTCGGTCCGCGAAAAGCTGCGGTCGAAGCCGCTAGTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAgg                                                                                                                                                                                                                                                                        <  1:243464/279‑1 (MQ=14)
                  tCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                                                                                                                                                                                                                                                                   <  2:736891/289‑1 (MQ=255)
                   ccACGTATAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGt                                                                                                                                                                                                                                                                       >  2:349672/1‑284 (MQ=255)
                     aCGAAAAACGACCGTCGAAGAAGCGAGCGCCCGCTCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACGGGTCGATCCGCGCAAAACCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCACGCAAGCTGGAACAGCAANAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcg                                                                                                                                                                                                                                                                <  2:70466/289‑1 (MQ=14)
                      cGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGCGCCACAGCCCCACAGCTGGAACAGCAACCGGCTAATGCCGAGCCAGAACAACAGGTCGATCcg                                                                                                                                                                                                                                                                >  2:781633/1‑288 (MQ=255)
                      cGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCAGGAACAGCAACAAGCTAATGCCGCGCCAGAACAACAGGTCGATcgg                                                                                                                                                                                                                                                                >  2:223549/1‑286 (MQ=255)
                                                   ccGCGCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc                                                                                                                                                                                                                                   <  2:480518/288‑1 (MQ=255)
                                                      cgcCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                                                                                                                                                                               >  2:86790/1‑289 (MQ=255)
                                                                         ggAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCGAGAAGAACAGGTCGAGCCGCGCAAAGCCGCCGTCGAACCCGCTCTTGCCCGAGCCAAAGCGC‑GCACACgg                                                                                                                                                                                                            >  2:176637/1‑289 (MQ=255)
                                                                                    gggCAAATGCGGACCCAGCAACACAGGCCGTACCGCGTAAAGCCGCCGTCGAAGCCGCTATCGCCCGCGCCAAAGCGCGCAAGCGGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCaa                                                                                                                                                                                                   <  2:94443/286‑1 (MQ=14)
                                                                                               gAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCgg                                                                                                                                                                                      >  2:583110/1‑289 (MQ=255)
                                                                                                     gAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCCAAGCCCGCAAACCGGACCAGCAACCGGCTAATGCCGAGCCCGAACAACAGGTCGAACCACCCAAAGCCGCAGTTGCACCCCCTATTGCCCGCACCAAAGCCC‑ACAAACGGGCACCGCAACCGGCCGATGCGGAACCaa                                                                                                                                                                                >  1:817989/1‑288 (MQ=14)
                                                                                                                  aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGa                                                                                                                                                                   >  2:782116/1‑289 (MQ=255)
                                                                                                                   tCCGGGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGATAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGAt                                                                                                                                                                  <  2:597124/289‑1 (MQ=255)
                                                                                                                            aactgccgTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCGTGCCAAAGCACGCAAGCCGGAACAGCAACAGGCTAATGCGGACCCAGAAGAACAGGTCAAACCGCCCAAAGCCGCCGTTGACCCGGCTATTGCCCACGCCAAA‑CGCAGCAAACCGGCCACGCAATCGGCTGTTGCCGAGCCCGAAAACCAggccaaaccaaccaa                                                                                                                                                          >  1:277457/5‑275 (MQ=11)
                                                                                                                                 ccgTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCACGCAAAGCTg                                                                                                                                                     >  2:134888/1‑288 (MQ=255)
                                                                                                                                        agcagcCATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGa                                                                                                                                              <  2:315607/288‑1 (MQ=255)
                                                                                                                                              tATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCCGCt                                                                                                                                       >  2:710184/1‑289 (MQ=255)
                                                                                                                                                     aaggcagaaccggccaaaacgacaCAGCAACAGGCGACTGACGGACACGAAGGCCCGGCCGATGCGGGAAAAGCTGCGGTAGAAGCGGCCATTGCCGGGGCAAAAGCACGCAAGCGGGAACAGCAACAGGCAAATGCCGAGGCAGAACACCAGGTCGATCCGCGAAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCa                                                                                                                                 <  1:623084/265‑1 (MQ=255)
                                                                                                                                                      gggCCAAAGCGCGCAAGCGGGAACAGCAACAGGCAAAAGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCAcg                                                                                                                               <  2:618922/287‑1 (MQ=255)
                                                                                                                                                        gCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAACGCGC‑GCAAACGGGCACCGCAACCGGCTCATGCGGACCCAGATGAACCGGTTGACCCGCGCAAAGCTGCCGTCGCAGCCGCTATCGCAccc                                                                                                                              >  2:865824/1‑286 (MQ=14)
                                                                                                                                                             aGCGCGCAAGCTGGAACAGCAACAGGCTAATGAGGAACCAGAACGACAGGTCGATCCGCGCAAAACCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCaaa                                                                                                                        <  2:77742/289‑1 (MQ=255)
                                                                                                                                                                 cgcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCACGCAAACCGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGGGCAAAGCCGCCGTCCACGCCGCTATTGCCCGCGCCCAAGCCC‑GCAAACCAGAACAGCAACAGGCTGATACGGCTCCAGTAGACCAGGCTGATCCGCGCAAAGCCGCCGtttcctctgcccttgcacgctccacagcc                                                                                                                     >  1:169248/1‑260 (MQ=11)
                                                                                                                                                                                cAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                                                                                                                                                                                                                                                                   <  1:580785/132‑1 (MQ=255)
                                                                                                                                                                                cAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                                                                                                                                                                                                                                                                   >  2:580785/1‑132 (MQ=255)
                                                                                                                                                                                cAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCCGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCGGCTCTTGCCCGAGCCAAATCAC‑GCAAACCCGAACAGCAACCGGCTAATGCGGAGCAAGAAGAACAGGTTGATCCGCCCAAAGCGGCGGTCGCTCCGGCCATTTCCCGCCCCAAAGCATGCCAGcccgacccg                                                                                                      >  1:230582/1‑280 (MQ=14)
                                                                                                                                                                                          aTGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCaa                                                                                                                                                                                                   >  1:480501/1‑185 (MQ=255)
                                                                                                                                                                                          aTGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCaa                                                                                                                                                                                                   <  2:480501/185‑1 (MQ=255)
                                                                                                                                                                                           tGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAACCAGAAGAACAGGTTGATCCGCCCAAACCTGCCGTCCAAGCGGCTATTGCCCGCTCCAAAGCACGCAAGCTGGAACCGCAAAAGGCCaaa                                                                                          >  1:437081/1‑288 (MQ=17)
                                                                                                                                                                                           tGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAAt                                                                                          >  2:105294/1‑289 (MQ=255)
                                                                                                                                                                                                agccAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCCGCTATTACACGCGCCAAAGCAC‑GCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGt                                                                                    >  2:170853/3‑289 (MQ=11)
                                                                                                                                                                                                aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTTGAAGCCGCCATTGCCCGGGCCAAAGCCC‑GCAAACTGGAACAGCAACAGGCTAATACCGAACACGAAGAACAAGTCAATACGCGCAAAGCCGCCGTTTCCGCGGCTCTTTCCCGCACCAAAGCCCTCAAGCCGGAACAGCaaacagctaacacggt                                                                                     >  1:258484/1‑275 (MQ=11)
                                                                                                                                                                                                aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGt                                                                                    >  2:347584/1‑289 (MQ=255)
                                                                                                                                                                                                aaCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCCGTCGACGCGGCTATTGCACGCGCCAAAGCACCACAGCAGGAACAGCAAACGGCTAATGC‑gg                                                                                     >  2:46733/1‑288 (MQ=255)
                                                                                                                                                                                                 aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGt                                                                                    <  2:437081/288‑1 (MQ=255)
                                                                                                                                                                                                       aaaaCAGGCCGACCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGAAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCag                                                                               <  2:728515/285‑1 (MQ=255)
                                                                                                                                                                                                          aCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAgaagaa                                                                          >  2:686878/1‑289 (MQ=255)
                                                                                                                                                                                                                ggcggctgcaacggccgtgtgcGAAGCCGCTATTGACCGTGCCAAAACACGCAAGCTGGAACCGCAACGGGCTAAGGCCGCGACAGAACAACAGGTCGAGCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGtt                                                                    <  1:362647/268‑1 (MQ=17)
                                                                                                                                                                                                                   tCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGAGGAGCCCGAAGCACAGGTTGATCCGCGCAAAGCTGCCCTCGACGCTGCTCATGCCCCCGACAAAGCACTCAAGCCGGAACTACCACAAGCTCATGCTGTTACCGTAAGAACACGtgca                                                                  >  1:764702/1‑286 (MQ=17)
                                                                                                                                                                                                                     cggcaactgccgccGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATcc                                                               >  2:362647/9‑289 (MQ=255)
                                                                                                                                                                                                                      gcgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATcc                                                               >  2:306286/1‑288 (MQ=255)
                                                                                                                                                                                                                           aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCaa                                                         <  2:147513/289‑1 (MQ=255)
                                                                                                                                                                                                                                        aGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAAAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGc                                             >  2:832879/1‑288 (MQ=255)
                                                                                                                                                                                                                                          ccGCTAGTGCCCGTGCCAAAGCAGGAAGGCGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAGGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcg                                           <  2:679663/288‑1 (MQ=255)
                                                                                                                                                                                                                                           ggCTGTTGCCCGGGCCAAAGAGCGCAAGCGGGAACCGGAACCGGCTAATGCGGGGCCAGAAGAACAGGATGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAGGAGACGCGGCAAGGGTC‑GCAAGCGGAAAAAGCAAAAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcg                                           <  1:170853/286‑1 (MQ=14)
                                                                                                                                                                                                                                              tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCt                                      >  2:329880/1‑289 (MQ=255)
                                                                                                                                                                                                                                                       tGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcgc                             >  2:491019/1‑289 (MQ=255)
                                                                                                                                                                                                                                                        gCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGCACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCCGTCGAAGCGGCTACTGCACGCGCCCAAGCAAGCAAGCCGGAACAGCAAAAGGCTATTCC‑GGTACCCGCAGAACCGGTTGCTACGCCCAAATCGGCATTTGCCGCGGCCATCGCCcgc                             >  1:649446/1‑288 (MQ=17)
                                                                                                                                                                                                                                                             aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAg                       >  2:397155/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                  gCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaa                  >  2:632463/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                  gCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCACAAGCAAGCAAGCTGGCACAGCAACAGGCTCATGC‑GGTACCAGAAGAAACGGTCGCTCCGCGCAAACCCGCAGTTGCCGCGGCCAATGCCCGGCCTCAGGCCat                  >  1:751719/1‑288 (MQ=17)
                                                                                                                                                                                                                                                                           accAGCAACAGCCACATGCGGAGCAAGAACAACCGGTCGAGCGCCGCAAAGCCGCCGTCGATGCGGCTATTGCCCGAGCCAAAGCGG‑GGAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTGGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgccg          <  1:105294/286‑1 (MQ=18)
                                                                                                                                                                                                                                                                            aCAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCACGCAAAGCCGCCGTGGAAGCCGCTAGTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAAAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgccgcc        <  2:280264/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aCAGGCTAATGCGGAGCAAGAACAACAGGGCGATCCGCGGAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCcagcag  <  2:169248/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aCAGGCCAAGGCCGAGCCAGAACAACAGGTCGATCCGGGCAAGGCCGCCGTCGAAGCCGCAATGGCCCGGGCCAAAGCAC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCcagcag  <  1:397155/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        tAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGc                                                                                                                                                                                                                                            <  2:142433/51‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        tAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaaa                                                                                                                                                                                                                         >  2:770624/1‑70 (MQ=255)
                                                                                                                                                                                                                                                                                        tAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaaa                                                                                                                                                                                                                         <  1:770624/70‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                      
CGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAG‑CTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGC‑GCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGC‑GGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCGCCCAGCAG  >  NC_000913/1708559‑1709120

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