Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,829,731 A→C 5.6% I229L (ATT→CTT)  srlR → DNA‑binding transcriptional repressor SrlR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,829,7310AC 5.6% 284.1 / 7.4 90I229L (ATT→CTT) srlRDNA‑binding transcriptional repressor SrlR
Reads supporting (aligned to +/- strand):  ref base A (38/47);  new base C (3/2);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

GGCGTAACCACCTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAGGAAGCAAGCCGT  >  NC_000913/2829587‑2829878
                                                                                                                                                |                                                                                                                                                   
ggCGTAACCACCTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTAcc                                                                                                                                                >  1:1161326/1‑150 (MQ=255)
ggCGTAACCACCTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTAcc                                                                                                                                                <  1:1342900/150‑1 (MQ=255)
 gCGTAACCACCTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGTCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCg                                                                                                                                               >  2:578261/1‑150 (MQ=255)
      accaccTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCa                                                                                                                                          >  1:701390/1‑150 (MQ=255)
      accaccTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCa                                                                                                                                          <  1:54585/150‑1 (MQ=255)
        caccTTTAACGATGTTTATACCTTCATTAAGGAAATGTGCAATGCCGCGCGCTAAGTGATTTTTATGGCGGACTAATCAAAGTTTGTCCTTAAAAGCCCCAACGTAGTTTGCATTCTTGAAAGAGTCTATAAGCTGATTACCGACGAAgg                                                                                                                                        <  2:4570/150‑1 (MQ=255)
        caccTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAgg                                                                                                                                        <  1:978507/150‑1 (MQ=255)
          ccTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGGTAATCAGCTTATCGACGctttca                                                                                                                                      <  2:710325/150‑6 (MQ=255)
           cTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGGTAATCAGCTTATCGACGctttcaa                                                                                                                                     <  1:1512600/150‑7 (MQ=255)
             ttAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCg                                                                                                                                   <  2:1117034/150‑1 (MQ=255)
                aCGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATc                                                                                                                                >  1:1329632/1‑150 (MQ=255)
                 cGAGGTTTATACCGTCAGTAAGTCAATGTGCAATGTCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATcc                                                                                                                               <  1:1042164/150‑1 (MQ=255)
                     gTTTATACCGTCAGTAAGGCAATGTGCCATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAACGCCCCAACGTAGTTTGCAGTCTTGAAAGCGGCGATAAGCTGATTACCGACGCAGGTAGCGATCCGGCg                                                                                                                           >  2:1577834/1‑150 (MQ=255)
                       ttATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGtt                                                                                                                         <  1:1422490/150‑1 (MQ=255)
                        tataCCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTc                                                                                                                       >  1:1381785/1‑150 (MQ=255)
                         ataCCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTc                                                                                                                       >  2:94300/1‑150 (MQ=255)
                          taCCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCg                                                                                                                      >  1:874660/1‑150 (MQ=255)
                            ccGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGTGGACTCATCAAAGTTTGGCCTTAAAAGCCCCATCGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTAGCGACGCAGGTATCGATCCGGCGTTTCGTc                                                                                                                    <  2:505095/150‑1 (MQ=255)
                            ccGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTc                                                                                                                    <  2:792373/150‑1 (MQ=255)
                              gTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAg                                                                                                                  <  2:114087/150‑1 (MQ=255)
                              gTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAg                                                                                                                  <  1:631296/150‑1 (MQ=255)
                               tCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAgg                                                                                                                 <  2:794447/150‑1 (MQ=255)
                               tCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAgg                                                                                                                 <  1:707460/150‑1 (MQ=255)
                                cAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGgc                                                                                                                >  1:1181775/1‑150 (MQ=255)
                                    aaGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTg                                                                                                            >  1:411642/1‑150 (MQ=255)
                                      ggCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGa                                                                                                                                              >  2:556852/1‑114 (MQ=255)
                                      ggCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGa                                                                                                                                              <  1:556852/114‑1 (MQ=255)
                                       gCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGaa                                                                                                         <  1:1580628/150‑1 (MQ=255)
                                       gCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGaa                                                                                                         <  2:40398/150‑1 (MQ=255)
                                        cAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAag                                                                                                        <  2:439313/150‑1 (MQ=255)
                                              gCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAaga                                                                                                       <  2:825992/145‑1 (MQ=255)
                                              gCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAaga                                                                                                       >  1:825992/1‑145 (MQ=255)
                                              gCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCAGTAAAAGCCCCAACGTAGTTTGCAGTGTTGAAAGCGTCGATAAGCTGAGTACCGACGCAGGTATCGCTCCGGAGTTTCGTCAGGCGCTGGAAGAGAAAg                                                                                                  >  2:1211857/1‑150 (MQ=255)
                                                aaTGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAggg                                                                                                <  2:1633308/150‑1 (MQ=255)
                                                   gCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATc                                                                                             <  2:113780/150‑1 (MQ=255)
                                                     cgcgcgcgAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGa                                                                                           <  2:874660/150‑1 (MQ=255)
                                                      gcgcgcgAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGAt                                                                                          <  1:1400010/150‑1 (MQ=255)
                                                      gcgcgcgAAGTGATTTTGATGGCGGACTCATCAAAGTTTGACCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGAt                                                                                          >  1:960258/1‑150 (MQ=255)
                                                       cgcgcgAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGAtg                                                                                         <  1:369955/150‑1 (MQ=255)
                                                        gcgcgAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGAtgt                                                                                        >  1:1318999/1‑150 (MQ=255)
                                                         cgcgAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAATAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGAtgtg                                                                                       <  1:1054695/150‑1 (MQ=255)
                                                           cgAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGAt                                                                                     >  2:1521583/1‑150 (MQ=255)
                                                           cgAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGAt                                                                                     >  2:886294/1‑150 (MQ=255)
                                                            gAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATc                                                                                    >  1:1309071/1‑150 (MQ=255)
                                                             aaGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCa                                                                                   >  2:985207/1‑150 (MQ=255)
                                                             aaGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCa                                                                                   >  1:1589120/1‑150 (MQ=255)
                                                             aaGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCa                                                                                   <  1:657278/150‑1 (MQ=255)
                                                                 gATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAAc                                                                               >  1:1630636/1‑150 (MQ=255)
                                                                 gATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAAc                                                                               >  1:746420/1‑150 (MQ=255)
                                                                     ttGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGga                                                                           >  1:1336965/1‑150 (MQ=255)
                                                                      tGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGGGTTTCGTCAGGCGCTTGAAGAGAAAGGGATCGATGTTATCTTAActtgag                                                                          <  2:701390/150‑6 (MQ=255)
                                                                      tGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGgag                                                                          <  2:1630636/150‑1 (MQ=255)
                                                                      tGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGgag                                                                          >  2:1584437/1‑150 (MQ=255)
                                                                          ggCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGagagag                                                                      <  1:1601487/150‑1 (MQ=255)
                                                                           gCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGc                                                                     <  2:120305/150‑1 (MQ=255)
                                                                               aCTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGGTGGCAGTCTTGAAAGCGTCGATAAGCTGAGTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATg                                                                 <  2:1529482/150‑1 (MQ=255)
                                                                                cTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGa                                                                <  2:1161326/150‑1 (MQ=255)
                                                                                  catcaAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGt                                                              >  2:38268/1‑150 (MQ=255)
                                                                                   atcaAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCaa                                                                   >  1:1013436/1‑144 (MQ=255)
                                                                                   atcaAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCaa                                                                   <  2:1013436/144‑1 (MQ=255)
                                                                                   atcaAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTg                                                             >  2:252688/1‑150 (MQ=255)
                                                                                   atcaAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTg                                                             >  2:629290/1‑150 (MQ=255)
                                                                                     caAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGaa                                                           >  2:1377767/1‑150 (MQ=255)
                                                                                       aaGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGc                                                         <  1:1053862/150‑1 (MQ=255)
                                                                                          tttGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCact                                                      <  1:1584437/150‑1 (MQ=255)
                                                                                              gCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTg                                                  >  2:832776/1‑150 (MQ=255)
                                                                                                cGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGaa                                                >  2:204182/1‑150 (MQ=255)
                                                                                                cGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGaa                                                <  1:395195/150‑1 (MQ=255)
                                                                                                        ccccAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGAc                                        <  1:94300/150‑1 (MQ=255)
                                                                                                                      gCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATg                          <  1:163721/150‑1 (MQ=255)
                                                                                                                        aGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCt                        <  2:1329632/150‑1 (MQ=255)
                                                                                                                         gTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTg                       <  2:960258/150‑1 (MQ=255)
                                                                                                                             tacctgCGTCGGTAATCAGCTTATCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGt                   >  2:725616/6‑150 (MQ=255)
                                                                                                                             tGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGt                   >  2:44976/1‑150 (MQ=255)
                                                                                                                               aaaGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTg                 <  2:1336965/150‑1 (MQ=255)
                                                                                                                               aaaGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTg                 >  1:1540985/1‑150 (MQ=255)
                                                                                                                                aaGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGc                >  2:793855/1‑150 (MQ=255)
                                                                                                                                 tgCGTCGGTAATCAGCTTATCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCa               >  1:348555/2‑150 (MQ=255)
                                                                                                                                 tgCGTCGGTAATCAGCTTATCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCa               >  1:601003/2‑150 (MQ=255)
                                                                                                                                 aGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCa               <  2:1381785/150‑1 (MQ=255)
                                                                                                                                  gCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAg              <  2:1181775/150‑1 (MQ=255)
                                                                                                                                    gTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAGGa            <  1:832776/150‑1 (MQ=255)
                                                                                                                                      cGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAGGaag          >  1:1307605/1‑150 (MQ=255)
                                                                                                                                       gATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAGGaagc         >  2:1588794/1‑150 (MQ=255)
                                                                                                                                          aaGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAGGaagcaag      >  1:830169/1‑150 (MQ=255)
                                                                                                                                           aGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGa                                                            >  1:908693/1‑95 (MQ=255)
                                                                                                                                           aGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGa                                                            <  2:908693/95‑1 (MQ=255)
                                                                                                                                           aGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAGGaagcaagc     <  2:1540985/150‑1 (MQ=255)
                                                                                                                                              tGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAGGAAGCAAGCCGt  <  1:886294/150‑1 (MQ=255)
                                                                                                                                              tGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAGGAAGCAAGCCGt  <  1:1521583/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGCGTAACCACCTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAGGAAGCAAGCCGT  >  NC_000913/2829587‑2829878

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: