Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,221,130 G→A Q40* (CAA→TAA)  ygjH ← putative tRNA‑binding protein YgjH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,221,1300GA95.7% 192.2 / ‑2.3 69Q40* (CAA→TAA) ygjHputative tRNA‑binding protein YgjH
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (21/45);  total (22/47)
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.27e-01

GGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTC  >  NC_000913/3220883‑3221372
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
ggATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTAccc                                                                                                                                                                                                                                                 <  1:50672/251‑1 (MQ=255)
 gATGGCGATACGACGTCGCTATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTAcccc                                                                                                                                                                                                                                                <  1:330893/251‑1 (MQ=255)
 gATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTAcccc                                                                                                                                                                                                                                                <  1:569593/251‑1 (MQ=255)
       gATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCATGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCAACAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATc                                                                                                                                                                                                                                           >  1:413131/1‑250 (MQ=255)
        aTACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCa                                                                                                                                                                                                                                          >  1:483922/1‑250 (MQ=255)
          acgacgTCTCTATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAAc                                                                                                                                                                                                                                        <  1:362326/250‑1 (MQ=255)
                  gCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGAt                                                                                                                                                                                                                               <  1:106041/251‑1 (MQ=255)
                  gCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGAGTTAACAGCACCCTTCCAATGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGAt                                                                                                                                                                                                                               >  1:39828/1‑251 (MQ=255)
                    cATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGt                                                                                                                                                                                                                             >  1:39058/1‑251 (MQ=255)
                           tCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACCGCtt                                                                                                                                                                                                                      >  2:203100/1‑251 (MQ=255)
                              ttttATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCtt                                                                                                                                                                                                                      <  1:486618/248‑1 (MQ=255)
                              ttttATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGt                                                                                                                                                                                                                    <  1:425014/250‑1 (MQ=255)
                              ttttATACGACGCGTACGCCCGCAGGCATCAGCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGTTTTCCGCACACAACAGCATGCATGCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTTCAGCGTTTTTTACCCCACAGCAACCTGTACGATGTACAGCtt                                                                                                                                                                                                                      <  1:543533/248‑1 (MQ=255)
                                    acgacgCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGtt                                                                                                                                                                                                             <  2:194847/251‑1 (MQ=255)
                                       acgCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTTTCGGCGTTTTc                                                                                                                                                                                                          <  2:457528/251‑1 (MQ=255)
                                         gcgTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCAt                                                                                                                                                                                                        >  1:275202/1‑251 (MQ=255)
                                          cgTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATg                                                                                                                                                                                                       <  1:402393/251‑1 (MQ=255)
                                           gTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGGTTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATgg                                                                                                                                                                                                      >  2:215555/1‑251 (MQ=255)
                                             aCGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGc                                                                                                                                                                                                     <  1:532408/250‑1 (MQ=255)
                                                cccGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGttt                                                                                                                                                                                                 <  2:275202/251‑1 (MQ=255)
                                                 ccGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTc                                                                                                                                                                                                <  2:455922/251‑1 (MQ=255)
                                                   gCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACAAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCAc                                                                                                                                                                                              >  2:617775/1‑251 (MQ=255)
                                                                      acGTTAACAGAACGCTTTCACTGCCGTAATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCACCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGAcg                                                                                                                                                                           >  2:197534/2‑251 (MQ=255)
                                                                        gTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGAcgcg                                                                                                                                                                         >  1:373802/1‑251 (MQ=255)
                                                                             cagcccgttttCACTGCCGTTTTCGGCTTCCGCACACAACCGCCTGCCTTCCGACGTTTCACCCCGCATTTTCGCTTTTTGCAGATTGCATACGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTCCGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCATCCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCAttt                                                                                                                                                                    >  2:213822/9‑251 (MQ=255)
                                                                                      ttCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACg                                                                                                                                                             <  2:177338/249‑1 (MQ=255)
                                                                                            gCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCaaaa                                                                                                                                                      <  1:193959/250‑1 (MQ=255)
                                                                                             ccGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATc                                                                                                                                                    <  2:450432/251‑1 (MQ=255)
                                                                                              cGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCCCCGTTTTCCCCATCAGCTCTTCTTCGCGGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAAGCTTTCCGACGCGCATGTCCAGACGGGCAAAATCa                                                                                                                                                   <  2:264573/251‑1 (MQ=255)
                                                                                                 cATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGc                                                                                                                                                 <  2:316391/250‑1 (MQ=255)
                                                                                                      gTTTCCGCCCACAACAGCATGCATTCCGACGTTTCACCACGCATTATCGCTTTTTGCAGATTGCAGAAGACCACCACCGTTTTCCCCATCAGCGCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTGTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGGTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTaa                                                                                                                                             <  2:42421/249‑1 (MQ=255)
                                                                                                        ttCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCGTTATGCAGATTTCATAAGACCACCCCCGTTTTCCCCATCAGCTCTTCCTCGCTGTAGTACGGCACCAGGCTGGTCACTGGTTGCAGCGTTTTTTACCCCACATCAACCTGGACGATGTACAGCTTGTCGTCGTTTTCATGGCGGTTCACGTCCACAAACTTTCCGACGCGCATTGCCAGACGTGCAAAATAAGCTTAAGcc                                                                                                                                          <  2:254487/250‑1 (MQ=255)
                                                                                                            gCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGt                                                                                                                                     <  2:138183/251‑1 (MQ=255)
                                                                                                            gCACACAACAGCATGCATTCCGACGTGTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCGTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCATCGTAAGCCACGGt                                                                                                                                     <  2:398982/251‑1 (MQ=255)
                                                                                                                     aGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTGCGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGGTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACTGTTTCCATTGc                                                                                                                            <  2:373802/251‑1 (MQ=255)
                                                                                                                     aGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGc                                                                                                                            >  2:165220/1‑251 (MQ=255)
                                                                                                                      gcatgcatTCCGACGTTTCACCACGCATTTTCGCTTTGTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCt                                                                                                                           <  1:46980/251‑1 (MQ=255)
                                                                                                                        atgcatTCCGACTTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTGCTCTGTAGTCCGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAACCAGCGTAAGCCTCGGTTTCCATTGCTtc                                                                                                                         <  2:383466/251‑1 (MQ=255)
                                                                                                                            atTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGTTTTCCATTGCTgcctcc                                                                                                                     <  1:501987/251‑1 (MQ=255)
                                                                                                                               ccGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCtt                                                                                                                  <  1:42832/251‑1 (MQ=255)
                                                                                                                                  aCGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGCCCACCTCCCTTTGCCCCACCAGCTCTTCTTCGCTGTAGTACGGCACCAGTCTGGTCTCGGTTTGCAGCTTTTTTTACCCCACATCAACCTGTACGATGTGCAGCTGGTCGGCGTTTTCATGTCGTTTCGCTTCCACAATCGTTCCGACGCGCATTTCGAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAgg                                                                                                               <  1:37318/251‑1 (MQ=255)
                                                                                                                                  aCGTCTCACCCCGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCAGCAGCTCTTCTTCGCTGTAGTACGGCCCCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAgg                                                                                                               <  1:160979/251‑1 (MQ=255)
                                                                                                                                       tCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTaa                                                                                                           >  1:463242/1‑250 (MQ=255)
                                                                                                                                        caccacGCATTGTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTTTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCTCAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTaaa                                                                                                          <  2:53092/250‑1 (MQ=255)
                                                                                                                                        caccacGCATGTTCGCTTTTTGCAGATTGCATAAGACCCCCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTTGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTCTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCATTTCCACAATCTTTCCGTCGCGCATTTCCAGACGTGCAAAATCAGAGTAAGCCACGTTTTCCATTGCTTCCTCCCTTAggcgaaaa                                                                                                         <  2:605655/251‑7 (MQ=255)
                                                                                                                                         cccacGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAAt                                                                                                        <  2:87424/250‑1 (MQ=255)
                                                                                                                                                    tCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAAc                                                                                              <  1:20700/250‑1 (MQ=255)
                                                                                                                                                           ttGCAGGTTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCa                                                                                       <  1:215555/250‑1 (MQ=255)
                                                                                                                                                            tGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCaa                                                                                      >  1:485443/1‑250 (MQ=255)
                                                                                                                                                            tGCAGATCGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCGCGGTTTGCAGCGTTTTTTACCCCACAGCAACCTGTACGATGTGCAGCTTGTCGGCGTTTTCATGGCGTTGCACTGCCACAATCTTTCCGACGCGCATTGCCAGACGTGCAAAATCAGCGTAAGTCACGGGGTCCAGTGCTTCCTCCCGTATGGTAGAAGTTGTCTAAACTAGAGCaa                                                                                      <  2:413131/250‑1 (MQ=255)
                                                                                                                                                              cAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATc                                                                                                                                                                                    >  1:323680/1‑154 (MQ=255)
                                                                                                                                                              cAGATTGCATAAGACCACCACCGTTGTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTGTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATc                                                                                                                                                                                    <  2:323680/154‑1 (MQ=255)
                                                                                                                                                               aGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGACACGGTTTGCAGCGTTTTTTACCCCACCTCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGt                                                                                  >  1:276812/1‑251 (MQ=255)
                                                                                                                                                                      aTAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTtctc                                                                           <  1:112767/251‑1 (MQ=255)
                                                                                                                                                                        aaGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCaa                                                                         <  1:345923/251‑1 (MQ=255)
                                                                                                                                                                            ccaccaaggTTTTACCCTTCAGCTCTTCTTCGCTGTAGTACGGCACCCGGCTGGTCCCGGTTTGCAGCGTTTTTTACCCCACATCACCCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGCCGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCGTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCt                                                                     >  2:88121/9‑251 (MQ=255)
                                                                                                                                                                                          cccATCAGCGCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATg                                                       <  2:483922/251‑1 (MQ=255)
                                                                                                                                                                                            cATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAACATgg                                                      >  2:344519/1‑250 (MQ=255)
                                                                                                                                                                                               ccgCTCTTCTTCGCTTTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGgaga                                                   <  2:411100/248‑1 (MQ=255)
                                                                                                                                                                                                aGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCATCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTAAATCCTGTAGGCTCAAAATGGAGAAt                                                 >  1:28291/1‑251 (MQ=255)
                                                                                                                                                                                                aGCTCTGCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTTCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAAt                                                 <  1:570362/251‑1 (MQ=255)
                                                                                                                                                                                                 gCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGTCGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATg                                                >  2:569182/1‑251 (MQ=255)
                                                                                                                                                                                                    cttcttCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGc                                               >  2:437837/1‑249 (MQ=255)
                                                                                                                                                                                                          cGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTg                                        <  1:273504/250‑1 (MQ=255)
                                                                                                                                                                                                                                ggTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTCCAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGCCGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGATAAAAATGCAGAATGCAGGCGTGATCACATTCGTGAGCCGCTgtgt                 >  2:6266/1‑251 (MQ=255)
                                                                                                                                                                                                                                ggTCACGGTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCtgtg                  <  1:423848/250‑1 (MQ=255)
                                                                                                                                                                                                                                 gTCACGGTTTGCAGCGTTTTTTACCCCTCATCAACCTGTACGATGCACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCTCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGGGAGCCGCGGTGtt                <  2:463242/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gTTTGCAGCGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGtt          <  1:344519/251‑1 (MQ=255)
                                                                                                                                                                                                                                               cGTTTTTTACCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTc  <  1:41126/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
GGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACACAACAGCATGCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTC  >  NC_000913/3220883‑3221372

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑