Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 11,664 A→C D41E (GAT→GAG yaaI ← UPF0412 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091311,6640AC99.0% 356.0 / ‑4.1 97D41E (GAT→GAGyaaIUPF0412 family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (46/50);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-01

CTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGG  >  NC_000913/11516‑11810
                                                                                                                                                    |                                                                                                                                                  
cttTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTc                                                                                                                                                  >  1:857173/1‑151 (MQ=255)
cttTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTc                                                                                                                                                  >  2:131041/1‑151 (MQ=255)
   ttATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGt                                                                                                                                                <  2:71315/150‑1 (MQ=255)
   ttATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGt                                                                                                                                                <  2:253059/150‑1 (MQ=255)
    tATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGtt                                                                                                                                               >  2:798315/1‑150 (MQ=255)
    tATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTg                                                                                                                                              >  1:512990/1‑151 (MQ=255)
        tCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGtt                                                                                                                                           <  1:179768/150‑1 (MQ=255)
         cTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTc                                                                                                                                         >  2:983121/1‑151 (MQ=255)
          tGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCg                                                                                                                                        >  1:532250/1‑151 (MQ=255)
           gAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTc                                                                                                                                         <  1:292989/149‑1 (MQ=255)
           gAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTc                                                                                                                                         <  1:370535/149‑1 (MQ=255)
                 atatTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTAcgcg                                                                                                                                 >  2:135456/1‑151 (MQ=255)
                  tatTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCgg                                                                                                                                >  1:426696/1‑151 (MQ=255)
                   atTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGc                                                                                                                               >  1:573222/1‑151 (MQ=255)
                      cAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCAtt                                                                                                                            >  1:610837/1‑151 (MQ=255)
                            ctGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATg                                                                                                                      >  1:1083037/1‑151 (MQ=255)
                             tGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATgg                                                                                                                     >  1:405303/1‑151 (MQ=255)
                               aCTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGgcg                                                                                                                   >  1:243905/1‑151 (MQ=255)
                                  ggCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCg                                                                                                                <  2:125343/151‑1 (MQ=255)
                                     gcTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCgagg                                                                                                             >  1:775824/1‑151 (MQ=255)
                                        aCGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGt                                                                                                          <  1:683266/151‑1 (MQ=255)
                                        aCGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGt                                                                                                          <  1:421136/151‑1 (MQ=255)
                                         cGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGtt                                                                                                         >  1:386769/1‑151 (MQ=255)
                                          gggCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTGTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGttt                                                                                                        <  1:623589/151‑1 (MQ=255)
                                           ggcggcTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTa                                                                                                       <  2:621809/151‑1 (MQ=255)
                                           ggcggcTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTa                                                                                                       <  2:595705/151‑1 (MQ=255)
                                            gcggcTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTAt                                                                                                      >  1:843505/1‑151 (MQ=255)
                                            gcggcTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTAt                                                                                                      >  1:1234123/1‑151 (MQ=255)
                                             cggcTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATg                                                                                                     <  2:1284596/151‑1 (MQ=255)
                                                 tttGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTc                                                                                                  >  1:351170/1‑150 (MQ=255)
                                                 tttGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCg                                                                                                 >  2:347642/1‑151 (MQ=255)
                                                    gAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGCGTTTATGGTCGttt                                                                                              >  1:1121701/1‑151 (MQ=255)
                                                    gAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGttt                                                                                              <  1:878919/151‑1 (MQ=255)
                                                     aaaTAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTg                                                                                             >  2:1105486/1‑151 (MQ=255)
                                                       aTAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGcc                                                                                           >  2:975751/1‑151 (MQ=255)
                                                         aaaCCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTg                                                                                         <  2:1198530/151‑1 (MQ=255)
                                                         aaaCCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTg                                                                                         <  2:577514/151‑1 (MQ=255)
                                                               aTGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGt                                                                                   >  2:1268611/1‑151 (MQ=255)
                                                                  cACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTgc                                                                                 >  2:1150550/1‑150 (MQ=255)
                                                                   aCCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTgc                                                                                 >  1:117049/1‑149 (MQ=255)
                                                                    ccGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTgcagc                                                                              >  1:134098/1‑151 (MQ=255)
                                                                     cGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTgcagc                                                                              >  2:1248089/1‑150 (MQ=255)
                                                                        ttAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGcagcac                                                                            >  2:235025/1‑149 (MQ=255)
                                                                          aaCTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGcagcacagca                                                                        <  1:731964/151‑1 (MQ=255)
                                                                              gTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATc                                                                      >  1:933375/1‑149 (MQ=255)
                                                                              gTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCa                                                                     <  1:1004230/150‑1 (MQ=255)
                                                                                aaaTCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGc                                                                  >  1:907646/1‑151 (MQ=255)
                                                                                  aTCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTa                                                                <  2:22220/151‑1 (MQ=255)
                                                                                     gCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATc                                                             <  2:1002100/151‑1 (MQ=255)
                                                                                      ccATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATc                                                             >  1:211808/1‑150 (MQ=255)
                                                                                        aTGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGcc                                                          >  1:981630/1‑151 (MQ=255)
                                                                                        aTGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGcc                                                          >  2:1001995/1‑151 (MQ=255)
                                                                                           gTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAgg                                                       <  2:177134/151‑1 (MQ=255)
                                                                                            tCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGc                                                      >  2:952975/1‑151 (MQ=255)
                                                                                             cGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCt                                                     <  2:214443/151‑1 (MQ=255)
                                                                                              ggCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTg                                                    <  2:236329/151‑1 (MQ=255)
                                                                                               gCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTgg                                                   <  1:840063/151‑1 (MQ=255)
                                                                                                 agagagTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCg                                                 <  1:576587/151‑1 (MQ=255)
                                                                                                   agagTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGa                                               <  2:374777/151‑1 (MQ=255)
                                                                                                     agTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGaa                                              >  1:1169776/1‑150 (MQ=255)
                                                                                                         gTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGt                                         >  2:1004526/1‑151 (MQ=255)
                                                                                                           aTGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTaa                                       <  1:1082671/151‑1 (MQ=255)
                                                                                                             gCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCTGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTaaa                                      <  2:626551/150‑1 (MQ=255)
                                                                                                             gCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTaaa                                      >  2:593513/1‑150 (MQ=255)
                                                                                                             gCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTaaa                                      <  1:492208/150‑1 (MQ=255)
                                                                                                                tttCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACg                                  >  2:912599/1‑151 (MQ=255)
                                                                                                                 ttCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACg                                  <  1:1228910/150‑1 (MQ=255)
                                                                                                                 ttCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACg                                  <  2:791394/150‑1 (MQ=255)
                                                                                                                  tCACAATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACgg                                 <  2:430149/150‑1 (MQ=255)
                                                                                                                   cacaATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGa                                <  1:746324/150‑1 (MQ=255)
                                                                                                                    acaATGCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGAtt                              >  1:572372/1‑151 (MQ=255)
                                                                                                                         gCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATa                         <  2:266400/151‑1 (MQ=255)
                                                                                                                         gCGACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATa                         <  1:600499/151‑1 (MQ=255)
                                                                                                                           gACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATa                         >  1:10070/1‑149 (MQ=255)
                                                                                                                           gACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAg                       <  1:132823/151‑1 (MQ=255)
                                                                                                                           gACAAACAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAg                       <  1:389053/151‑1 (MQ=255)
                                                                                                                               aaCAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGAtt                   >  1:1238458/1‑151 (MQ=255)
                                                                                                                                 cAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATtct                 <  1:813143/151‑1 (MQ=255)
                                                                                                                                 cAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATtct                 >  2:733227/1‑151 (MQ=255)
                                                                                                                                 cAGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATtct                 >  1:721661/1‑151 (MQ=255)
                                                                                                                                  aGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATtctc                <  2:476459/151‑1 (MQ=255)
                                                                                                                                  aGGAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATtctc                <  2:1121771/151‑1 (MQ=255)
                                                                                                                                  aGGAAGTTTCAGCGCCAGCTCGGGGGTTGCGGGACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGAtggaagaaggaggata                <  2:322917/151‑16 (MQ=255)
                                                                                                                                   ggAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTctct               >  1:1270255/1‑151 (MQ=255)
                                                                                                                                   ggAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTctct               <  1:14700/151‑1 (MQ=255)
                                                                                                                                    gAAGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCtt              <  2:443532/151‑1 (MQ=255)
                                                                                                                                      aGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAg            <  1:58855/151‑1 (MQ=255)
                                                                                                                                      aGTTTCAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAg            <  1:949004/151‑1 (MQ=255)
                                                                                                                                           cAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTggg        <  1:758833/150‑1 (MQ=255)
                                                                                                                                           cAGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTgg         >  2:1065740/1‑149 (MQ=255)
                                                                                                                                            aGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTggg        >  1:1159288/1‑149 (MQ=255)
                                                                                                                                            aGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGaa      <  2:553130/151‑1 (MQ=255)
                                                                                                                                            aGCGCCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGaa      <  1:951631/151‑1 (MQ=255)
                                                                                                                                              cgcCAGCGCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGaa      <  2:1276362/149‑1 (MQ=255)
                                                                                                                                               gcCAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAagag   <  2:940369/151‑1 (MQ=255)
                                                                                                                                                ccAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAgg  <  2:166112/151‑1 (MQ=255)
                                                                                                                                                 cAGCTCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAgg  >  1:618753/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCATGGTCGGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGG  >  NC_000913/11516‑11810

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