Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,463,609 T→G 8.2% F229C (TTC→TGC) ‡ gspK → Type II secretion system protein GspK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,463,6090TG 8.2% 292.6 / 16.7 98F229C (TTC→TGC) ‡gspKType II secretion system protein GspK
Reads supporting (aligned to +/- strand):  ref base T (50/40);  new base G (4/4);  total (54/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

CAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAGCCACTGGTTGCTCAGGTG  >  NC_000913/3463461‑3463757
                                                                                                                                                    |                                                                                                                                                    
cAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTc                                                                                                                                                     >  1:143612/1‑150 (MQ=255)
  gATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCag                                                                                                                                                   <  2:979884/150‑1 (MQ=255)
  gATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCag                                                                                                                                                   <  1:1923425/150‑1 (MQ=255)
    tGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCagag                                                                                                                                                 >  1:998965/1‑150 (MQ=255)
     gCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATGCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGc                                                                                                                                                >  2:1349816/1‑150 (MQ=255)
     gCTTTTTCTGACTGACGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGATAACGATATTCCGCTCTTCAGAGc                                                                                                                                                <  2:1106836/150‑1 (MQ=255)
        ttttCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGtt                                                                                                                                             >  2:984586/1‑150 (MQ=255)
        ttttCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCCGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCATAATATGCTAACCCAGAACGATATTCCGCTCTTCAGACCgat                                                                                                                                             >  2:977466/1‑148 (MQ=255)
         tttCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTGTGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTa                                                                                                                                            <  2:1290265/150‑1 (MQ=255)
         tttCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTa                                                                                                                                            >  2:1417928/1‑150 (MQ=255)
          ttCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTAt                                                                                                                                           <  1:984586/150‑1 (MQ=255)
          ttCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTAt                                                                                                                                           <  1:977466/150‑1 (MQ=255)
          ttCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTGAAGAGCggaat                                                                                                                                           <  1:1430943/150‑5 (MQ=255)
          ttCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTGAAGAGCggaat                                                                                                                                           <  2:75790/150‑5 (MQ=255)
           tCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTAtt                                                                                                                                          <  2:227372/150‑1 (MQ=255)
             tGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTcc                                                                                                                                        <  2:1976267/150‑1 (MQ=255)
             tGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTcc                                                                                                                                        <  1:1654770/150‑1 (MQ=255)
             tGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTGAAGAGCggaatatc                                                                                                                                        <  2:429915/150‑8 (MQ=255)
                  ggCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAAt                                                                                                                                   >  1:1465307/1‑150 (MQ=255)
                  ggCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAAt                                                                                                                                   >  2:93769/1‑150 (MQ=255)
                          aCGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCAc                                                                                                                           <  2:1465307/150‑1 (MQ=255)
                           cGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACg                                                                                                                          >  2:1466671/1‑150 (MQ=255)
                           cGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACg                                                                                                                          >  1:1335596/1‑150 (MQ=255)
                            gTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACgg                                                                                                                         <  1:907584/150‑1 (MQ=255)
                            gTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACgg                                                                                                                         >  2:1960813/1‑150 (MQ=255)
                            gTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTGAAGAGCggaat                                                                                                                                           <  1:1590532/132‑5 (MQ=255)
                            gTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTGAAGAGCggaat                                                                                                                                           >  2:1590532/1‑128 (MQ=255)
                             tCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGa                                                                                                                        >  1:1593168/1‑150 (MQ=255)
                              cAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGAt                                                                                                                       >  1:127890/1‑150 (MQ=255)
                                taTAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATga                                                                                                                     <  1:481931/149‑2 (MQ=255)
                                  tAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTg                                                                                                                   <  2:1255745/150‑1 (MQ=255)
                                    aaaGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGAc                                                                                                                 >  2:1958910/1‑150 (MQ=255)
                                     aaGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACg                                                                                                                <  2:218455/150‑1 (MQ=255)
                                          aTAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCgg                                                                                                           >  2:772393/1‑150 (MQ=255)
                                           tAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCggg                                                                                                          >  2:802474/1‑150 (MQ=255)
                                            aaCCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGt                                                                                                         >  1:593812/1‑150 (MQ=255)
                                            aaCCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGt                                                                                                         >  2:1251580/1‑150 (MQ=255)
                                                gAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTa                                                                                                     <  1:567418/150‑1 (MQ=255)
                                                gAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTa                                                                                                     >  1:1543241/1‑150 (MQ=255)
                                                   aaCATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTg                                                                                                  <  2:301866/150‑1 (MQ=255)
                                                    aCATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGc                                                                                                 <  1:802474/150‑1 (MQ=255)
                                                    aCATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGc                                                                                                 <  2:593812/150‑1 (MQ=255)
                                                     cATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCa                                                                                                <  1:1369439/150‑1 (MQ=255)
                                                     cATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCa                                                                                                >  2:392676/1‑150 (MQ=255)
                                                      aTCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCaa                                                                                               >  1:709967/1‑150 (MQ=255)
                                                           tCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAg                                                                                          <  1:1646510/150‑1 (MQ=255)
                                                           tCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAg                                                                                          <  1:941224/150‑1 (MQ=255)
                                                            cAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAgg                                                                                         <  2:1550155/150‑1 (MQ=255)
                                                                ggCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCggg                                                                                                          <  2:1175942/129‑1 (MQ=255)
                                                                ggCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCggg                                                                                                          >  1:1175942/1‑129 (MQ=255)
                                                                  cTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAgg                                                                                   >  1:511026/1‑150 (MQ=255)
                                                                  cTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAgg                                                                                   <  2:1712278/150‑1 (MQ=255)
                                                                   ttATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAggg                                                                                  <  1:392676/150‑1 (MQ=255)
                                                                       cccGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAgg                                                                              >  1:544529/1‑150 (MQ=255)
                                                                       cccGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATATTAATATGCTAACGGAGAACGATATTCCGATCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAgg                                                                              >  2:363900/1‑150 (MQ=255)
                                                                        ccGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAggc                                                                             >  1:557554/1‑150 (MQ=255)
                                                                          gTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAggctg                                                                           <  2:998965/150‑1 (MQ=255)
                                                                           tATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAggctgg                                                                          >  1:830477/1‑150 (MQ=255)
                                                                                      tctcCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGa                                                               >  1:162751/1‑150 (MQ=255)
                                                                                       ctcCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCGAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGAt                                                              <  1:1349816/150‑1 (MQ=255)
                                                                                       ctcCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGAt                                                              >  2:1887412/1‑150 (MQ=255)
                                                                                        tcCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATg                                                             <  1:1251580/150‑1 (MQ=255)
                                                                                        tcCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATg                                                             >  1:1555486/1‑150 (MQ=255)
                                                                                          ccGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCa                                                           >  1:1733840/1‑150 (MQ=255)
                                                                                               cACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCt                                                      >  1:368057/1‑150 (MQ=255)
                                                                                                aCGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCtt                                                     >  2:1635907/1‑150 (MQ=255)
                                                                                                aCGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCtt                                                     <  2:742971/150‑1 (MQ=255)
                                                                                                    aaCTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCCGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCGGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACt                                                 <  1:1043308/150‑1 (MQ=255)
                                                                                                      cTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTgg                                               >  1:388283/1‑150 (MQ=255)
                                                                                                           cATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCgca                                          >  2:605598/1‑150 (MQ=255)
                                                                                                           cATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCgca                                          <  1:1635907/150‑1 (MQ=255)
                                                                                                            aTCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCgcag                                         >  1:1634468/1‑150 (MQ=255)
                                                                                                              cAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCgcagca                                       >  2:1834181/1‑150 (MQ=255)
                                                                                                               aaTCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGcagcag                                      >  2:1879680/1‑150 (MQ=255)
                                                                                                                 tCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGa                                    >  2:398966/1‑150 (MQ=255)
                                                                                                                    tAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGAttt                                 >  2:532950/1‑150 (MQ=255)
                                                                                                                    tAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGAttt                                 >  1:248341/1‑150 (MQ=255)
                                                                                                                     aaTATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTtc                                <  2:1534532/150‑1 (MQ=255)
                                                                                                                       tatGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTtctc                              <  1:93769/150‑1 (MQ=255)
                                                                                                                         tGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCg                            >  2:1588950/1‑150 (MQ=255)
                                                                                                                          gCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCgg                           <  2:1335596/150‑1 (MQ=255)
                                                                                                                           cTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCgg                           >  2:370859/1‑149 (MQ=255)
                                                                                                                           cTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCgg                           <  1:370859/149‑1 (MQ=255)
                                                                                                                            tAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGgtg                         <  1:1277888/150‑1 (MQ=255)
                                                                                                                             aaCGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGgtgt                        >  1:852643/1‑150 (MQ=255)
                                                                                                                             aaCGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGgtgt                        >  1:1443154/1‑150 (MQ=255)
                                                                                                                             aaCGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGgtgt                        >  2:579577/1‑150 (MQ=255)
                                                                                                                                 gagaACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAg                    >  1:836657/1‑150 (MQ=255)
                                                                                                                                  agaACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAGc                   >  2:1804261/1‑150 (MQ=255)
                                                                                                                                   gaACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAGcc                  <  1:449852/150‑1 (MQ=255)
                                                                                                                                   gaAAGATATTCCGATCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTACGGTGTGAAGcc                  >  2:838280/1‑150 (MQ=255)
                                                                                                                                      cGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAGCCACt               <  1:1887412/150‑1 (MQ=255)
                                                                                                                                        gaataacGCTCTGAAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAGCCACTgg             >  2:207675/7‑150 (MQ=255)
                                                                                                                                         aataacGCTCTGAAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGACTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAGCCACTGGt            >  1:911173/6‑150 (MQ=255)
                                                                                                                                           taccGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATCCTGACGCTAGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCAAAACACCGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAGCCACTGGTTg          >  2:1014256/3‑150 (MQ=255)
                                                                                                                                            tCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAGCCACTGGTTGc         <  1:772393/150‑1 (MQ=255)
                                                                                                                                             acGCTCTGAAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAGCCACTGGTTGCt        >  1:272241/2‑150 (MQ=255)
                                                                                                                                                 tctTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAGCCACTGGTTGCTCAgg    <  2:368057/150‑1 (MQ=255)
                                                                                                                                                   ttCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCTAAAAAGGCCGAGGGAAGGCTGGCTAACCAGGTATGCATTTCTTTACTGGGCTCAGCAGGATTTCTCCGGTGTGAAGCCACTGGTTGCTCAGGTg  <  2:1593168/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGCTCTTCAGAGCGTTATTCCTCAATAACATCACGGATGCTGACGCTCGGGTCTTATTGCAAAAAAGGCCGAGGGAAGGCTGGCTAACCACGGATGCATTTCTTTACTGGGCGCAGCAGGATTTCTCCGGTGTGAAGCCACTGGTTGCTCAGGTG  >  NC_000913/3463461‑3463757

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: