Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,671,938 T→C 5.1% intergenic (+78/‑22) ydgU → / → ydgD stationary phase‑induced protein/putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,671,9380TC 5.1% 250.5 / 2.9 80intergenic (+78/‑22)ydgU/ydgDstationary phase‑induced protein/putative peptidase
Reads supporting (aligned to +/- strand):  ref base T (37/38);  new base C (2/2);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGATGAGGTCAGCACCCTGTTTTTTGGTCATGA  >  NC_000913/1671794‑1672084
                                                                                                                                                |                                                                                                                                                  
cTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAAt                                                                                                                                                 >  1:209231/1‑148 (MQ=255)
cTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAAt                                                                                                                                                 >  1:76904/1‑148 (MQ=255)
cTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAAt                                                                                                                                                 >  7:210656/1‑148 (MQ=255)
    gAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTTTTTATGAGGGTTTGCTTTTAATAATCATAATTAcc                                                                                                                                             <  6:74936/148‑1 (MQ=255)
    gAGTTCATTCTGATCATCCTTATATTATGCGCACTGATTACAGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGAGTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTAcc                                                                                                                                             >  6:354266/1‑148 (MQ=255)
    gAGTTCATTCGGATCATCCTTATTTTATGCGCACTGAATCCCCCCCGTTTTTATCTTTCCAGATTGAAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGGTTGCTTTTAATAATCATAATTAcc                                                                                                                                             <  2:122446/148‑1 (MQ=255)
     aGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGGCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCACGTCTATAGTATTTATGAGGGTGTGCTTTTAATAATCATAATTACcca                                                                                                                                           <  5:214286/149‑1 (MQ=255)
         cATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACccacca                                                                                                                                        <  2:209231/148‑1 (MQ=255)
           ttCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCagag                                                                                                                                     <  8:201279/149‑1 (MQ=255)
             cTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAgt                                                                                                                                    <  4:64025/148‑1 (MQ=255)
               gATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAgtgt                                                                                                                                  >  1:307213/1‑148 (MQ=255)
                atcatcCTTATTTTATGCGCACTGATTACCGCCCGTTATTATCTTTCCTGATTGTAGTTATCAGATTTTACTCACACTTTCATCCCGTCCCGACTAAAGTATATATTAGGGTTTGCTATTAATAATCATAATTACCCACCAGAGtctg                                                                                                                                 >  4:204398/1‑148 (MQ=255)
                   atgcTTTTTTTATTCTCACTGATTAGCGCACGTTTTTATGTTTCGTTATTGTAGTTATCTGTTTTTTCTCCCACTTTCATCCCGTCCCGTCTATAGTTTTTATGAGGGTTTGCTTTTAATTTTGATAATTACCCACCAGAGTGTGatat                                                                                                                             <  3:204398/146‑1 (MQ=255)
                      cTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGc                                                                                                                           <  2:17268/148‑1 (MQ=255)
                      cTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCg                                                                                                                          >  8:141551/1‑149 (MQ=255)
                         aTTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGt                                                                                                                         <  6:3414/147‑1 (MQ=255)
                          ttttATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTAc                                                                                                                       <  3:267675/148‑1 (MQ=255)
                          ttatATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTAc                                                                                                                       >  7:189847/4‑148 (MQ=255)
                            ttATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAAc                                                                                                                    >  7:321090/1‑149 (MQ=255)
                                 cgcACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTg                                                                                                               <  2:76904/149‑1 (MQ=255)
                                  gcACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGc                                                                                                              >  8:303335/1‑149 (MQ=255)
                                    aCTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATaa                                                                                                                                                  <  8:235983/111‑1 (MQ=255)
                                    aCTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATaa                                                                                                                                                  >  7:235983/1‑111 (MQ=255)
                                         ttACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAgtg                                                                                                       <  8:321090/149‑1 (MQ=255)
                                          tACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAgtgt                                                                                                      >  6:296563/1‑149 (MQ=255)
                                            ccGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGtt                                                                                                     <  7:303335/148‑1 (MQ=255)
                                             cGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACACACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTgg                                                                                                   >  3:202135/1‑149 (MQ=255)
                                               cccGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCGCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCAt                                                                                                                                                    <  4:84293/98‑1 (MQ=255)
                                               cccGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTTTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAgtgt                                                                                                                                  >  7:97707/1‑116 (MQ=255)
                                               cccGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTTTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAgtgt                                                                                                                                  <  8:97707/116‑1 (MQ=255)
                                               cccGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCAt                                                                                                                                                    >  3:84293/1‑98 (MQ=255)
                                                ccGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTg                                                                                                >  2:167005/1‑149 (MQ=255)
                                                   tttttATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCCCTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCACAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCaa                                                                                             >  3:17396/1‑149 (MQ=255)
                                                   tttttATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCaa                                                                                             >  1:349951/1‑149 (MQ=255)
                                                   tttttATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCACAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCaa                                                                                             >  3:12203/1‑149 (MQ=255)
                                                           tttCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTg                                                                                                               >  8:308829/1‑123 (MQ=255)
                                                           tttCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTg                                                                                                               <  7:308829/123‑1 (MQ=255)
                                                             tCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTAccc                                                                                                                                            <  8:103854/92‑1 (MQ=255)
                                                             tCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGttt                                                                                      <  8:189847/146‑1 (MQ=255)
                                                             tCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCAAAATTAccc                                                                                                                                            >  7:103854/1‑92 (MQ=255)
                                                               cTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGt                                                                                 >  7:329379/1‑149 (MQ=255)
                                                                  aTTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCt                                                                               <  2:349951/148‑1 (MQ=255)
                                                                  aTTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCACAATTACCCACCAGAGTGTGATATGCGTACAACGATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTc                                                                                <  4:12203/147‑1 (MQ=255)
                                                                  aTTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCACAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTc                                                                                <  4:17396/147‑1 (MQ=255)
                                                                         ttATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTtgt                                                                       <  7:141551/149‑1 (MQ=255)
                                                                         ttATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTGGGGTGCAATTAGTTTGACGTCTGCTTTtgt                                                                       <  5:354266/149‑1 (MQ=255)
                                                                            tCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGtt                                                                    >  1:159733/1‑149 (MQ=255)
                                                                               gATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTg                                                                  >  2:232491/1‑148 (MQ=255)
                                                                                aTTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGttt                                                                   <  2:307213/146‑1 (MQ=255)
                                                                                 ttttACTCCCACTTTCATCCCGTCACGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCa                                                                <  4:61754/148‑1 (MQ=255)
                                                                                       tCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGc                                                                                               >  1:303259/1‑111 (MQ=255)
                                                                                       tCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGc                                                                                               <  2:303259/111‑1 (MQ=255)
                                                                                       tCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATaaa                                                          >  3:21266/1‑148 (MQ=255)
                                                                                        cccACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAAc                                                         >  8:64851/1‑148 (MQ=255)
                                                                                         ccACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCt                                                                            <  8:370614/128‑1 (MQ=255)
                                                                                         ccACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCt                                                                            >  7:370614/1‑128 (MQ=255)
                                                                                         ccACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAAc                                                         >  8:397905/1‑147 (MQ=255)
                                                                                               tCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGtt                                                 >  1:232154/1‑149 (MQ=255)
                                                                                                 aTCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGc                                               <  1:232491/149‑1 (MQ=255)
                                                                                                          cgtcTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCgg                                       >  5:14274/1‑148 (MQ=255)
                                                                                                           gtcTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAAt                                                                                            >  8:328159/1‑94 (MQ=255)
                                                                                                            tcTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCaa                                    <  7:64851/149‑1 (MQ=255)
                                                                                                                taGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACACACCAGAGTGTGATATGCGTACAACCATTTCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAAc                                  <  8:210656/147‑1 (MQ=255)
                                                                                                                     tttATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCtttt                                                                         >  3:320996/1‑103 (MQ=255)
                                                                                                                     tttATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCtttt                                                                         <  4:320996/103‑1 (MQ=255)
                                                                                                                     tttATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGATGAg                            <  5:296563/148‑1 (MQ=255)
                                                                                                                       tATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGATGAGGt                          <  4:21266/148‑1 (MQ=255)
                                                                                                                           aGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGt                                                                                                 >  5:61253/1‑73 (MQ=255)
                                                                                                                           aGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGt                                                                                                 <  6:61253/73‑1 (MQ=255)
                                                                                                                            gggTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTgg                                                                                                   >  5:292614/1‑70 (MQ=255)
                                                                                                                            gggTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTgg                                                                                                   <  6:292614/70‑1 (MQ=255)
                                                                                                                            gggTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTg                                                                  >  7:249085/1‑103 (MQ=255)
                                                                                                                            gggTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTg                                                                  <  8:249085/103‑1 (MQ=255)
                                                                                                                              gTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGATGAGGTTAGCAccc                  <  6:14274/149‑1 (MQ=255)
                                                                                                                                  gCTTTTAATAATCATAATTACCCACCAGAGTGTGGTATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGATGAGGTCAGCACCCTGtt              >  3:325788/1‑149 (MQ=255)
                                                                                                                                  gCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGATGAGGTc                         <  3:293110/138‑1 (MQ=255)
                                                                                                                                  gCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTGTGCAGATAAGCCAGACGTTGCCAGATCGGGAAACGATGTGGTc                         >  4:293110/1‑138 (MQ=255)
                                                                                                                                    ttttAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGATGAGGTCAGCACCCTGtttt            >  7:10391/1‑149 (MQ=255)
                                                                                                                                             tCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGATGAGGTCAGCACCCTGTTTTTTGGTCat    >  8:67043/1‑148 (MQ=255)
                                                                                                                                               aTAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGa                                >  1:272368/1‑118 (MQ=255)
                                                                                                                                               aTAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGa                                <  2:272368/118‑1 (MQ=255)
                                                                                                                                               aTAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGATGAGGTCAGCACCCTGTTTTTTGGTCatga  <  2:159733/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CTTTGAGTTCATTCTGATCATCCTTATTTTATGCGCACTGATTACCGCCCGTTTTTATCTTTCCTGATTGTAGTTATCTGATTTTACTCCCACTTTCATCCCGTCCCGTCTATAGTATTTATGAGGGTTTGCTTTTAATAATCATAATTACCCACCAGAGTGTGATATGCGTACAACCATTGCTGTAGTGTTGGGTGCAATTAGTTTGACGTCTGCTTTTGTGTTTGCAGATAAACCAGACGTTGCCAGATCGGCAAACGATGAGGTCAGCACCCTGTTTTTTGGTCATGA  >  NC_000913/1671794‑1672084

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