Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 283,338 G→A M46I (ATG→ATA yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913283,3380GA97.0% 1486.1 / 21.6 441M46I (ATG→ATAyagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (5/8);  new base A (219/208);  total (225/216)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

Only 100 of 441 total aligned reads displayed.

GGACGATAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCGGCATGCAG  >  NC_000913/283190‑283485
                                                                                                                                                    |                                                                                                                                                   
ggACGATAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAAc                                                                                                                                                   >  8:944398/1‑151 (MQ=255)
ggACGATAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAAc                                                                                                                                                   >  6:511437/1‑151 (MQ=255)
   cGATAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCa                                                                                                                                                <  3:943932/151‑1 (MQ=255)
    gATAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAg                                                                                                                                               <  5:5366/151‑1 (MQ=255)
      tAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGa                                                                                                                                              >  7:124490/1‑150 (MQ=255)
         cGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGc                                                                                                                                          >  7:979804/1‑151 (MQ=255)
         cGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGc                                                                                                                                          >  2:792194/1‑151 (MQ=255)
             gACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCggg                                                                                                                                      <  1:1214642/151‑1 (MQ=255)
             gACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCggg                                                                                                                                      <  2:692092/151‑1 (MQ=255)
               ccATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGAt                                                                                                                                    >  1:30155/1‑151 (MQ=255)
               ccATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGAt                                                                                                                                    >  2:632564/1‑151 (MQ=255)
               ccATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGAt                                                                                                                                    <  2:863258/151‑1 (MQ=255)
                cATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATg                                                                                                                                   >  1:172741/1‑151 (MQ=255)
                      gaAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTgg                                                                                                                             <  6:1256003/151‑1 (MQ=255)
                         aaTTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCCACCTGGTCGGCATAACGCAGAACGCCGGGATGGGCTGGGAc                                                                                                                          >  2:1136433/1‑151 (MQ=255)
                            tttCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGcc                                                                                                                       >  2:1181196/1‑151 (MQ=255)
                             ttCACCCCGCAGGACGACGCGGTTTATGCGGTGAGCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGcc                                                                                                                       <  8:467125/150‑1 (MQ=255)
                              tCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCca                                                                                                                      >  5:473060/1‑150 (MQ=255)
                               cACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCcaac                                                                                                                    >  2:1685663/1‑151 (MQ=255)
                               cACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCcaac                                                                                                                    <  4:1336967/151‑1 (MQ=255)
                                 ccccGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCcaaca                                                                                                                   >  6:321033/1‑150 (MQ=255)
                                  cccGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCcaaca                                                                                                                   >  5:1261832/1‑149 (MQ=255)
                                  cccGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAg                                                                                                                 <  4:1086879/151‑1 (MQ=255)
                                   ccGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGc                                                                                                                <  1:869913/151‑1 (MQ=255)
                                     gCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTc                                                                                                              <  7:268141/151‑1 (MQ=255)
                                     gCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTc                                                                                                              <  1:606172/151‑1 (MQ=255)
                                      cAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTc                                                                                                              >  4:898121/1‑150 (MQ=255)
                                       aGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCAc                                                                                                            >  8:294311/1‑151 (MQ=255)
                                         gacgacGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCg                                                                                                          >  8:625986/1‑151 (MQ=255)
                                          ccgtcgcgTTTTATGCGGTGAGCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCgg                                                                                                         <  6:1702573/147‑1 (MQ=255)
                                          ccgacgCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCgg                                                                                                         <  6:842440/150‑1 (MQ=255)
                                            gacgCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGc                                                                                                        >  6:449483/1‑150 (MQ=255)
                                            gacgCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGc                                                                                                        >  1:1153699/1‑150 (MQ=255)
                                            gacgCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCa                                                                                                       <  6:1106984/151‑1 (MQ=255)
                                                cgTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCaa                                                                                                      >  4:1636215/1‑148 (MQ=255)
                                                   tttATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGt                                                                                                >  1:956032/1‑151 (MQ=255)
                                                    ttATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACACGCTCACCGGCAAAGAGGTg                                                                                               >  8:1382824/1‑151 (MQ=255)
                                                         cGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGAt                                                                                          >  8:1650120/1‑151 (MQ=255)
                                                              aTCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGACCGCCGGGATGGGCTGGGACGCCAACAAGCTCCCCGGCAAAGAGGTGCTGATTATCg                                                                                     >  2:1604759/1‑151 (MQ=255)
                                                              aTCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCg                                                                                     <  3:808668/151‑1 (MQ=255)
                                                                cACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCgg                                                                                    <  7:1282157/150‑1 (MQ=255)
                                                                cACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGc                                                                                   >  7:1118712/1‑151 (MQ=255)
                                                                 aCCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGAGGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCa                                                                                  <  1:371916/151‑1 (MQ=255)
                                                                 aCCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGc                                                                                   >  4:398221/1‑150 (MQ=255)
                                                                  cccACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCAc                                                                                 >  3:1490953/1‑151 (MQ=255)
                                                                   ccACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACt                                                                                >  3:1095864/1‑151 (MQ=255)
                                                                    cACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGCTTATCGGCACTc                                                                               >  4:807475/1‑151 (MQ=255)
                                                                      cgcgGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAg                                                                             <  6:901618/151‑1 (MQ=255)
                                                                       gcggcgGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCGCCGGCAAAGAGGTGCTGATTATCGGCACTCAgg                                                                            >  6:636439/1‑151 (MQ=255)
                                                                       gcggcgGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAgg                                                                            >  4:1393807/1‑151 (MQ=255)
                                                                         ggcggGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGggc                                                                          <  8:371487/151‑1 (MQ=255)
                                                                         ggcggGGCCGCAGGGCGCGCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGggc                                                                          <  3:186841/151‑1 (MQ=255)
                                                                          gcggGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGggcg                                                                         <  2:94052/151‑1 (MQ=255)
                                                                           cggGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGggcgg                                                                        >  7:1055320/1‑151 (MQ=255)
                                                                           cggGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGggcgg                                                                        >  5:1646053/1‑151 (MQ=255)
                                                                            ggggCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGggcggc                                                                       <  1:1053356/151‑1 (MQ=255)
                                                                                  gCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCcgc                                                                 <  4:293747/151‑1 (MQ=255)
                                                                                       gcgcTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCgg                                                            <  8:696519/151‑1 (MQ=255)
                                                                                        ggcTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGga                                                           <  4:1528654/150‑1 (MQ=255)
                                                                                        cgcTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACACGCTCACCGGGAAAGAGGTGCTGATTATCGGCCGTCCGGGCGGCCTCCGCGCCggc                                                           >  8:1180438/1‑150 (MQ=255)
                                                                                        cgcTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGga                                                           <  7:499451/151‑1 (MQ=255)
                                                                                            ctGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAgacg                                                       >  5:1556537/1‑151 (MQ=255)
                                                                                              gCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAgccgg                                                      >  3:417922/1‑150 (MQ=255)
                                                                                              gCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAgacgg                                                      <  1:1408503/150‑1 (MQ=255)
                                                                                              gCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAgacg                                                       <  7:1446025/149‑1 (MQ=255)
                                                                                              gCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAgacg                                                       <  3:337348/149‑1 (MQ=255)
                                                                                              gCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAgac                                                        <  8:1296335/148‑1 (MQ=255)
                                                                                                cGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCATCAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAgacggacg                                                   >  6:1573770/1‑151 (MQ=255)
                                                                                                 gCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGc                                                  >  4:704257/1‑151 (MQ=255)
                                                                                                   tGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGccc                                                >  1:516592/1‑151 (MQ=255)
                                                                                                       cccGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATc                                            <  3:1308533/151‑1 (MQ=255)
                                                                                                       cccGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATc                                            <  7:3987/151‑1 (MQ=255)
                                                                                                             gATGCTGATGGAACCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTg                                      <  2:1418690/151‑1 (MQ=255)
                                                                                                               tGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTggg                                    >  2:1270753/1‑151 (MQ=255)
                                                                                                               tGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTggg                                    >  2:1333939/1‑151 (MQ=255)
                                                                                                                   gATGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTa                                 <  5:279357/150‑1 (MQ=255)
                                                                                                                    aTGGAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTAcc                               >  5:663029/1‑151 (MQ=255)
                                                                                                                      ggAATCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACca                              <  1:764620/150‑1 (MQ=255)
                                                                                                                        aaTCTCCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCCCGCAGGGCGGCAACCGCGCCGGAGACGGACGCCCAATCCCGCTGGGCTACCccac                           >  2:376048/1‑151 (MQ=255)
                                                                                                                           ctcCCAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCgg                        >  7:573720/1‑151 (MQ=255)
                                                                                                                             cccAGCGGCAAGCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGGGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCggg                       <  8:1258075/150‑1 (MQ=255)
                                                                                                                             cccAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACANGCTCACNGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGc                      <  6:1383494/151‑1 (MQ=255)
                                                                                                                              ccAGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGc                      >  3:5120/1‑150 (MQ=255)
                                                                                                                                aGCGGCAACCTGTTCGGCATAAGGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCTGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCAtt                   <  6:1508545/151‑1 (MQ=255)
                                                                                                                                aGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCAtt                   <  3:601921/151‑1 (MQ=255)
                                                                                                                                aGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCAtt                   >  8:1610817/1‑151 (MQ=255)
                                                                                                                                aGCGGCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCAtt                   >  5:871568/1‑151 (MQ=255)
                                                                                                                                    gCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGga               <  7:831581/151‑1 (MQ=255)
                                                                                                                                    gCAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGga               <  2:575693/151‑1 (MQ=255)
                                                                                                                                     cAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTggg                >  5:1436145/1‑149 (MQ=255)
                                                                                                                                     cAACCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGgag              <  8:810695/151‑1 (MQ=255)
                                                                                                                                      aaCCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTggg                >  6:751575/1‑148 (MQ=255)
                                                                                                                                      aaCCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTggg                >  3:848995/1‑148 (MQ=255)
                                                                                                                                      aaCCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTggg                >  4:527355/1‑148 (MQ=255)
                                                                                                                                       aCCTGTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTggg                >  5:11632/1‑147 (MQ=255)
                                                                                                                                           gTTCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCgg         <  1:1571155/150‑1 (MQ=255)
                                                                                                                                             tCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCGGCa       <  5:1092821/150‑1 (MQ=255)
                                                                                                                                             tCGGCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCGGCa       <  7:1628631/150‑1 (MQ=255)
                                                                                                                                                gCATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCGGCatgca   <  7:881756/151‑1 (MQ=255)
                                                                                                                                                 cATAACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCGGCatgc    >  2:508922/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGACGATAACGATGACCATTGAGAAAATTTTCACCCCGCAGGACGACGCGTTTTATGCGGTGATCACCCACGCGGCGGGGCCGCAGGGCGCTCTGCCGCTGACCCCGCAGATGCTGATGGAATCTCCCAGCGGCAACCTGTTCGGCATGACGCAGAACGCCGGGATGGGCTGGGACGCCAACAAGCTCACCGGCAAAGAGGTGCTGATTATCGGCACTCAGGGCGGCATCCGCGCCGGAGACGGACGCCCAATCGCGCTGGGCTACCACACCGGGCATTGGGAGATCGGCATGCAG  >  NC_000913/283190‑283485

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