Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,829,727 G→T 5.6% K227N (AAG→AAT 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,7270GT 5.6% 303.1 / 5.2 89K227N (AAG→AATsrlRDNA‑binding transcriptional repressor SrlR
Reads supporting (aligned to +/- strand):  ref base G (38/46);  new base T (3/2);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

TCAATGCGGGCGTAACCACCTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAGGAAGCAAGC  >  NC_000913/2829579‑2829875
                                                                                                                                                    |                                                                                                                                                    
tcAATGCGGGCGTAACCACCTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGc                                                                                                                                                     <  1:992297/150‑1 (MQ=255)
        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                                                                                                               <  1:631296/150‑1 (MQ=255)
                                      gTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAg                                                                                                               <  2:114087/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                                                                                                    >  1:825992/1‑145 (MQ=255)
                                                      gCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAaga                                                                                                    <  2:825992/145‑1 (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:886294/1‑150 (MQ=255)
                                                                   cgAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGAt                                                                                  >  2:1521583/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)
                                                                                                                                                    |                                                                                                                                                    
TCAATGCGGGCGTAACCACCTTTAACGAGGTTTATACCGTCAGTAAGGCAATGTGCAATGCCGCGCGCGAAGTGATTTTGATGGCGGACTCATCAAAGTTTGGCCGTAAAAGCCCCAACGTAGTTTGCAGTCTTGAAAGCGTCGATAAGCTGATTACCGACGCAGGTATCGATCCGGCGTTTCGTCAGGCGCTGGAAGAGAAAGGGATCGATGTGATCATAACCGGAGAGAGCAATGAGTGAAGCACTACTGAACGCGGGACGTCAGACGTTAATGCTGGAGTTGCAGGAAGCAAGC  >  NC_000913/2829579‑2829875

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: