Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 281,965 A→T 7.4% W7R (TGG→AGG)  yagA ← CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,9650AT 7.4% 1300.6 / 4.3 478W7R (TGG→AGG) yagACP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (258/136);  new base T (2/49);  total (260/218)
Fisher's exact test for biased strand distribution p-value = 1.67e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 492 total aligned reads displayed.

GGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTT  >  NC_000913/281816‑282113
                                                                                                                                                     |                                                                                                                                                    
ggcACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                                                                      <  2:7109901/150‑1 (MQ=255)
ggcACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                                                                      <  2:791567/150‑1 (MQ=255)
ggcACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                                                                      <  2:5274761/150‑1 (MQ=255)
 gcACCTTCCTGAGCCCAGCGCTGGAGCCACTTGGAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAg                                                                                                                                                     <  2:5803852/150‑1 (MQ=255)
  cACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAgg                                                                                                                                                    >  2:1604109/1‑150 (MQ=255)
   aCCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAggg                                                                                                                                                   >  1:8145760/1‑150 (MQ=255)
   aCCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAggg                                                                                                                                                   >  1:66827/1‑150 (MQ=255)
    ccTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGc                                                                                                                                                  >  1:2713765/1‑150 (MQ=255)
     cTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCa                                                                                                                                                 >  2:1072994/1‑150 (MQ=255)
      ttCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCAt                                                                                                                                                <  1:1980975/150‑1 (MQ=255)
      ttCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCAt                                                                                                                                                <  2:2978915/150‑1 (MQ=255)
      ttCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCAt                                                                                                                                                <  1:3882739/150‑1 (MQ=255)
         cTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAg                                                                                                                                             >  1:2885755/1‑150 (MQ=255)
          tGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGc                                                                                                                                            <  1:7190868/150‑1 (MQ=255)
                cAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTcc                                                                                                                                      >  2:3408461/1‑150 (MQ=255)
                  gcgcTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCAt                                                                                                                                    >  2:2401456/1‑150 (MQ=255)
                   cgcTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATa                                                                                                                                   >  2:7837732/1‑150 (MQ=255)
                      tGGAGCCAGTTGTAGCCGGTGGCAGGTGAAATGCCGACGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCTGGGCATAAGCGACTCCATAAAc                                                                                                                                <  2:4575137/150‑1 (MQ=255)
                        gAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACgg                                                                                                                              >  1:6134073/1‑150 (MQ=255)
                         aGCCACTTGTTGCGGGTGGCTGGTGAAATGCCTAAGCGACGGCTGAGGTAAGGGATGTTCGCCCCGTCCTGCGAGGCTAACAAAACAAACTCAGTACTTAATGACATGGTATCTCTCGCATCCCTGGGAATAAGCGACTCCATAAACggg                                                                                                                             <  2:4312290/150‑1 (MQ=255)
                           ccACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGtt                                                                                                                           >  1:3303066/1‑150 (MQ=255)
                           ccACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGtt                                                                                                                           >  1:2442147/1‑150 (MQ=255)
                           ccACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGtt                                                                                                                           >  1:1830904/1‑150 (MQ=255)
                            cACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTc                                                                                                                          >  1:3482106/1‑150 (MQ=255)
                             aCTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCt                                                                                                                         >  1:1474904/1‑150 (MQ=255)
                             aCTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCt                                                                                                                         >  2:3345285/1‑150 (MQ=255)
                               ttGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGTCATGGTATCTCTCGCATCCCAGGGCATATGCGACTCCATAAACGGGTTCTTa                                                                                                                       >  2:669811/1‑150 (MQ=255)
                                  tAGCCGTTGCCAGGTTAAATGCCCAACCCACGCCAGAGGTAACGGATGTTCGCCCCGACCTCCTAGGCCAAAAAAACAAACTCACAACTTAATCACATGGTATCTATCGAAACCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGc                                                                                                                    <  2:1806741/150‑1 (MQ=255)
                                     ccGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCtt                                                                                                                 >  1:1149776/1‑150 (MQ=255)
                                     ccGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCtt                                                                                                                 >  1:316911/1‑150 (MQ=255)
                                              ggTGAAATGCCGAAGCGACGGCAGAGGGACCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTCTGCCTTAGTTGTAAg                                                                                                        >  1:161877/1‑150 (MQ=255)
                                               tttaaaTGCCGACGCGAGGGCAGTGGTAAGGTTTGTTCGCCCCGTCCTGCGTGGCTAAAAAAAAAAACTCAGTACTTAATGACTTGTTATCTCCCGCATCCCCGGGCAAAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAgt                                                                                                       <  1:442880/147‑1 (MQ=255)
                                               gTGATATGCCGAAGCGACGGCTGTGGGAACGGTTGTTCGCCCGGTCCTGCGAGGCTAAAAAAACAAACTCAGTACGTAATGACTTGGTATCTCTCGCATCCCAGGGAATACGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAgt                                                                                                       <  2:3209235/150‑1 (MQ=255)
                                               gTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAgt                                                                                                       >  1:1543034/1‑150 (MQ=255)
                                               gTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAgt                                                                                                       >  2:4995124/1‑150 (MQ=255)
                                                tGAAATGCCGAAGCGACGGAAGAGGCAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTAAGTACTTAATCACATGGTATATCTCGCATCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAgtg                                                                                                      <  2:5085715/150‑1 (MQ=255)
                                                  aaaTGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTc                                                                                                    >  1:6199551/1‑150 (MQ=255)
                                                      gCCGACGCGAGGGCTGGGGGACCGGTGGTCCGCCCCGTCCTGCGAGCATAAAAAAACAAACTCAGTACGTAATGAAATGGAATCCCGCGCACCCCTGGGAATACGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTAcc                                                                                                <  1:3642343/150‑1 (MQ=255)
                                                       ccGAAGCGACGGCAGCGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACTTGGTATCTCTCGCAGCCCTGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCa                                                                                               <  2:1476422/150‑1 (MQ=255)
                                                       ccGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGGTCTTATGCCTTAGTTGTAAGTGTCTACCa                                                                                               >  2:2208621/1‑150 (MQ=255)
                                                        cGACGCGCCGGCAGTGGGAACGGAGGTGCCCCCCGCCCTGCGCGGCGAACAAAACACACTCCGTACTTAATGACTGGTTATCTCTCGCATCCCGGGGCATACGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCAt                                                                                              <  2:2823169/150‑1 (MQ=255)
                                                         gAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATg                                                                                             >  2:1018450/1‑150 (MQ=255)
                                                              gACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCTAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTcccc                                                                                        <  2:4690746/150‑1 (MQ=255)
                                                              gACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGAATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTcccc                                                                                        <  2:7270692/150‑1 (MQ=255)
                                                                     gtggtaaCGGATGTTCGCCCCGTCCTGCGTGGCTAAAAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAg                                                                                 <  2:6343874/145‑1 (MQ=255)
                                                                     gagGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACACAGTACGTAATGACATGGTATCTCTCGCTTCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAg                                                                                 <  2:6867715/150‑1 (MQ=255)
                                                                       gggAAGGGATTGTCGCCCCGTCCTGCGTGGCTAAAAAAACAAACTCAGTACGTAATGAATTGTTATCTCTCGCATCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAgtg                                                                               <  1:7122842/150‑1 (MQ=255)
                                                                            ggtaTGTTCGCCCCGTGCTGGTTGGGTAACAAAACAAACACTGTACTTAATTACTTGGTATGTGTCGAATCCCTGGTCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCAc                                                                          <  1:4067180/147‑1 (MQ=255)
                                                                                tGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCCTAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATg                                                                      >  2:1977706/1‑150 (MQ=255)
                                                                                       cccGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCgg                                                               >  2:2597157/1‑150 (MQ=255)
                                                                                         cGTCCTGCTAGGCTAACAAAACAAACTCAGTCCGTAATGACATGTTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGAc                                                             <  2:4440431/150‑1 (MQ=255)
                                                                                             cTGCGAGGCTAAAAAAAAAAACTCAGTACGTAATGACATGTTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTa                                                         <  1:190284/150‑1 (MQ=255)
                                                                                             cTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACTTGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTa                                                         <  1:5804076/150‑1 (MQ=255)
                                                                                                cGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACAC                                                      >  1:792952‑M1/1‑150 (MQ=255)
                                                                                                   ggCGAAAAAAACAAACCAATAACTTAATGACATGGTATCTCGCGCATCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCC                                                   <  1:166789‑M1/150‑1 (MQ=255)
                                                                                                    gCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCA                                                  >  2:5407285‑M1/1‑150 (MQ=255)
                                                                                                    gCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCA                                                  >  1:5194719‑M1/1‑150 (MQ=255)
                                                                                                    gCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCA                                                  >  2:3717265‑M1/1‑150 (MQ=255)
                                                                                                     cGAACAAACCAAACTCTGTACGTAATGACATGGTATCTCTCGCAGCCCCGGGCATAAGCGACTCCATACACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAA                                                 <  2:4007712‑M1/150‑1 (MQ=255)
                                                                                                      gAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTGAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAA                                                >  1:2549242‑M1/1‑150 (MQ=255)
                                                                                                        acaaaacaaaCTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGG                                              >  2:3118919‑M1/1‑150 (MQ=255)
                                                                                                        acaaaacaaaCTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGG                                              >  2:3786662‑M1/1‑150 (MQ=255)
                                                                                                          aaaacaaaCTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGG                                            >  1:4814169‑M1/1‑150 (MQ=255)
                                                                                                            aacaaacTCAGTACGTAATGACATGGTATCTCTCGCATCACAGGGCATACGCGACTCCATAAACGGGTTCTTATGCCTGAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCACAGGGgag                                          >  2:7653698/1‑150 (MQ=255)
                                                                                                             acaaacTCAGTACGTAATGACATGGTATCTCTCGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGgaga                                         <  2:6057091/150‑1 (MQ=255)
                                                                                                              caaacTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGagag                                        >  2:5705531/1‑150 (MQ=255)
                                                                                                                   tCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGAc                                   >  2:4122737/1‑150 (MQ=255)
                                                                                                                      gTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGc                                >  2:1250188/1‑150 (MQ=255)
                                                                                                                        aCGTAATGACATGGTATCTCTCGCAGCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCAc                              <  1:7540796/150‑1 (MQ=255)
                                                                                                                           tattGACAGGGTAGCTCCCGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGa                           <  1:2999489/147‑1 (MQ=255)
                                                                                                                           tAATGACTTGGTACCTCGCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGa                           <  1:838422/150‑1 (MQ=255)
                                                                                                                           tAATGACATGGTAACTCGCGCATCCCAGGGCAAAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGa                           <  1:5415788/150‑1 (MQ=255)
                                                                                                                            aaTGACTTGGTATCCCCCGCACCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAg                          <  1:6405964/150‑1 (MQ=255)
                                                                                                                             aTGACTTGGTATCTCTCGCATCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGc                         <  2:4874299/150‑1 (MQ=255)
                                                                                                                              tGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGcc                        >  1:5209167/1‑150 (MQ=255)
                                                                                                                                aCATGGTATCTCTCGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCCGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCAt                      <  1:6863592/150‑1 (MQ=255)
                                                                                                                                 cATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATc                     >  1:6001923/1‑150 (MQ=255)
                                                                                                                                     gTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCtttt                 >  1:7824055/1‑150 (MQ=255)
                                                                                                                                     gTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCtttt                 >  2:277085/1‑150 (MQ=255)
                                                                                                                                     gTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCtttt                 >  2:2052661/1‑150 (MQ=255)
                                                                                                                                        tctctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTccc              >  2:5129851/1‑150 (MQ=255)
                                                                                                                                        tctcgcGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTccc              <  2:2516510/150‑1 (MQ=255)
                                                                                                                                          tctcGCATCCCCGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTccccc            <  2:81674/150‑1 (MQ=255)
                                                                                                                                           ctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCt           >  1:5711145/1‑150 (MQ=255)
                                                                                                                                           ctcGCAGCCCTGGGAATAAGCGACTCCATACACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCt           <  2:4236393/150‑1 (MQ=255)
                                                                                                                                            tcGCATCCCCGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTc          <  2:6152983/150‑1 (MQ=255)
                                                                                                                                            tcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTc          >  2:8027131/1‑150 (MQ=255)
                                                                                                                                            tcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTc          >  2:3525507/1‑150 (MQ=255)
                                                                                                                                            tcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTc          >  2:7551490/1‑150 (MQ=255)
                                                                                                                                              gCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGc        <  2:5593861‑M2/150‑2 (MQ=255)
                                                                                                                                               cATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcc       >  1:5519028‑M2/1‑148 (MQ=255)
                                                                                                                                                aTCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGccc      >  2:650080‑M2/1‑147 (MQ=255)
                                                                                                                                                 tCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccc     <  1:5000757‑M2/150‑5 (MQ=255)
                                                                                                                                                  ccccGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccct    >  2:6403645‑M2/1‑145 (MQ=255)
                                                                                                                                                  cccAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccct    <  2:1924176‑M2/150‑6 (MQ=255)
                                                                                                                                                   ccAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctt   <  1:1280238‑M2/150‑7 (MQ=255)
                                                                                                                                                   ccAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGccccgt   <  2:5904495‑M2/150‑7 (MQ=255)
                                                                                                                                                    cAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGccccttt  >  1:4639150‑M2/1‑143 (MQ=255)
                                                                                                                                                    cAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGccccttt  <  2:6223657‑M2/150‑8 (MQ=255)
                                                                                                                                                    cAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGccccttt  <  1:1053445‑M2/150‑8 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTT  >  NC_000913/281816‑282113

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