Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,409,250 (T)8→9 intergenic (+217/‑261) ydaN → / → dbpA putative Zn(II) transporter/ATP‑dependent RNA helicase, specific for 23S rRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,409,2421.T97.5% 279.1 / 0.0 79intergenic (+209/‑269)ydaN/dbpAputative Zn(II) transporter/ATP‑dependent RNA helicase, specific for 23S rRNA
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base T (34/43);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑TTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACTTATTTATTTGCAAATTCC  >  NC_000913/1409144‑1409336
                                                                                                   |                                                                                               
gTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑t                                                                                               >  2:1792944/1‑100 (MQ=255)
  cAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑ttt                                                                                             >  2:2387994/1‑100 (MQ=255)
  cAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑ttt                                                                                             >  1:578641/1‑100 (MQ=255)
     gCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑tttttt                                                                                          <  2:186791/100‑1 (MQ=255)
      cGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑ttttttt                                                                                         >  1:2762717/1‑100 (MQ=255)
       gttGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑ttttttn                                                                                         <  1:1013959/97‑2 (MQ=255)
       gATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGTATAGCACAGTTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑tttttttt                                                                                        <  2:1648137/100‑1 (MQ=255)
         tGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑tttttttttt                                                                                      <  2:2436484/100‑3 (MQ=255)
         tGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑tttttttttg                                                                                      <  2:182063/100‑3 (MQ=255)
         tGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑tttttttttg                                                                                      <  1:496162/100‑3 (MQ=255)
         tGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑tttttttttg                                                                                      <  1:1727943/100‑3 (MQ=255)
         tGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑tttttttttg                                                                                      <  1:1649085/100‑3 (MQ=255)
         tGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑tttttttttg                                                                                      >  2:1466471/1‑98 (MQ=255)
          gggCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑tttttttttgc                                                                                     >  2:2329196/1‑97 (MQ=255)
           ggCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑tttttttttgcc                                                                                    <  1:329314/100‑5 (MQ=255)
            gCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑tttttttttgcct                                                                                   >  1:657216/1‑95 (MQ=255)
             cGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTg                                                                                   <  2:1637417/100‑1 (MQ=255)
             cGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTg                                                                                   <  2:76764/100‑1 (MQ=255)
                 gcgcTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGttat                                                                               >  2:69470/1‑100 (MQ=255)
                 gcgcTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGttat                                                                               <  2:1085374/100‑1 (MQ=255)
                 gcgcTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGttat                                                                               >  1:2020666/1‑100 (MQ=255)
                  ctcTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGttatt                                                                              <  1:1430533/98‑1 (MQ=255)
                  cgcTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGttat                                                                               >  1:878952/1‑99 (MQ=255)
                    cTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTTCCTGttattta                                                                            >  2:292691/1‑100 (MQ=255)
                    cTCCAGATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGttattta                                                                            <  2:731488/100‑1 (MQ=255)
                     tCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGttatttat                                                                           <  2:1444994/100‑1 (MQ=255)
                      ccATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCTTGTTATTTATc                                                                          >  2:1571263/1‑100 (MQ=255)
                      ccATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATc                                                                          <  2:426944/100‑1 (MQ=255)
                        atatTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAg                                                                        >  1:1623506/1‑100 (MQ=255)
                          atTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCg                                                                      >  1:1268992/1‑100 (MQ=255)
                               cTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTAc                                                                 <  2:1623493/100‑1 (MQ=255)
                               cTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTAc                                                                 <  2:712036/100‑1 (MQ=255)
                                 tACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCTTCTAcc                                                                >  2:167243/1‑99 (MQ=255)
                                 tACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTAccc                                                               >  1:1894354/1‑100 (MQ=255)
                                 tACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTAccc                                                               <  1:545175/100‑1 (MQ=255)
                                 tACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTAccc                                                               >  1:1674381/1‑100 (MQ=255)
                                  aCTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCt                                                              >  1:92143/1‑100 (MQ=255)
                                   cTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCtt                                                             >  1:656839/1‑100 (MQ=255)
                                   cTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑TTTTTTTTGCCTGTTATTTATCAGCGTCTACCCttt                                                            <  2:1880006/100‑1 (MQ=255)
                                   cTTCCTTTTTCGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCtt                                                             <  1:1743350/100‑1 (MQ=255)
                                    ttCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCttt                                                            <  1:292694/100‑1 (MQ=255)
                                    ttCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCttt                                                            <  2:1171242/100‑1 (MQ=255)
                                      ccTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑TTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAag                                                         <  1:282133/100‑1 (MQ=255)
                                       cTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAag                                                         <  1:325795/100‑1 (MQ=255)
                                       cTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAag                                                         <  2:1644628/100‑1 (MQ=255)
                                       cTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAag                                                         <  2:76933/100‑1 (MQ=255)
                                         tttttGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAagag                                                       <  2:1362562/100‑1 (MQ=255)
                                           tttGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTc                                                     <  1:1376636/100‑1 (MQ=255)
                                           tttGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGACTGTTATTTATCAGTGTCTACCCTTTAACAGt                                                      >  2:2342260/1‑99 (MQ=255)
                                             tGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCa                                                   >  2:247228/1‑100 (MQ=255)
                                              gAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCTTATTATTTATTATCTTCTACCCTTTAAGAGTTCa                                                   >  2:1375030/1‑99 (MQ=255)
                                               aaTTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTTTCTTTTTTTTATTAGCGTCTACCCTTTAAGAGTTCAcc                                                 >  2:1832665/1‑100 (MQ=255)
                                                aTTACTGCATAGCTCAATTGATTCGTACGTCGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCAccc                                                <  1:1062062/100‑1 (MQ=255)
                                                aTTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCAccc                                                <  1:1792958/100‑1 (MQ=255)
                                                aTTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCAccc                                                <  2:1290984/100‑1 (MQ=255)
                                                 ttACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCa                                               >  1:1477116/1‑100 (MQ=255)
                                                 ttACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCa                                               <  2:1560750/100‑1 (MQ=255)
                                                  tACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCaa                                              <  2:1479654/100‑1 (MQ=255)
                                                   aCTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAAt                                             >  1:855760/1‑100 (MQ=255)
                                                    cggCATAGCACAAATGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATg                                            >  1:1872908/3‑100 (MQ=255)
                                                    cTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATg                                            >  1:1955990/1‑100 (MQ=255)
                                                      gCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTTTTTGTTATTTATCAGCGTCTCCCCTTTAAGAGTCCACCCAAt                                             >  2:1306719/1‑97 (MQ=255)
                                                      gCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTTCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGAc                                          >  2:587642/1‑100 (MQ=255)
                                                       cATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGAcc                                         <  2:855754/100‑1 (MQ=255)
                                                        aTAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCa                                        <  2:849242/100‑1 (MQ=255)
                                                        aTAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCa                                        <  2:2569877/100‑1 (MQ=255)
                                                           gCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCagag                                     >  2:858281/1‑100 (MQ=255)
                                                           gCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCagag                                     >  1:493766/1‑100 (MQ=255)
                                                            tcttccgatctCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAgg                                    <  1:1882687/90‑1 (MQ=255)
                                                                aTTGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGaaa                                <  1:2329214/100‑1 (MQ=255)
                                                                 ttGATTCGTACGCCGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGaaa                                <  1:2388013/99‑1 (MQ=255)
                                                                  tGATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTTCTTGTTTTTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAata                              >  2:43558/1‑100 (MQ=255)
                                                                   gATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAatat                             >  1:2548488/1‑100 (MQ=255)
                                                                   gATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAAGGAAatat                             >  2:2398281/1‑100 (MQ=255)
                                                                   gATTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAAAGTCCACCCAATGACCAGAGGGAAatat                             >  2:1663731/1‑100 (MQ=255)
                                                                    aTTCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATg                            <  1:1378445/100‑1 (MQ=255)
                                                                      tCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATgac                          <  2:1950761/100‑1 (MQ=255)
                                                                      tCGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATgac                          <  1:1181929/100‑1 (MQ=255)
                                                                       cGTACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTTGCCTGTTATTTATAAGCGTCTACCCTTTAAGAGACCACCCAATGACCAGAGGGAAATATga                           >  2:17330/1‑99 (MQ=255)
                                                                         tACGACGCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATgacgac                       <  2:1894339/100‑1 (MQ=255)
                                                                          acgacgCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGaca                      <  2:1955974/100‑1 (MQ=255)
                                                                           cgacgCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGacac                     >  1:110562/1‑100 (MQ=255)
                                                                              cgCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACtta                  <  2:2422244/100‑1 (MQ=255)
                                                                              cgCCGACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACtta                  <  1:1466482/100‑1 (MQ=255)
                                                                                ccGATTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACttatt                <  1:1832680/100‑1 (MQ=255)
                                                                                  gACTTTGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACttattta              >  1:2760913/1‑100 (MQ=255)
                                                                                   aCTTTGATGAGTCGGCTTTTTTTTTTTTTTTTTTTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACCCttatttat             >  2:2433868/1‑100 (MQ=255)
                                                                                     tttGATGAGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACTTatttattt           <  2:578637/100‑1 (MQ=255)
                                                                                       tgatgaGTCGGCTTTTTTTTTTTTTTTTTTTTATCAGCGTCTACCCTTTAAGATTCCACTCTATGTTCAGAGGTAAATATTACGACACTTatttattt           >  2:1434644/1‑98 (MQ=255)
                                                                                       tgatgaGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACTTATTTATTTGc         <  2:231909/100‑1 (MQ=255)
                                                                                       tgatgaGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACTTATTTATTTGc         <  1:2025852/100‑1 (MQ=255)
                                                                                          tgaGTCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACTTATTTATTTGCaaa      <  2:92141/100‑1 (MQ=255)
                                                                                             gTCGGCTTTTTTTTTTTTTTTTTTTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACTTATTTATTTGCAAAtt    >  2:2385538/1‑99 (MQ=255)
                                                                                              tCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACTTATTTATTTGCAAATTcc  <  1:1200368/100‑1 (MQ=255)
                                                                                              tCGGCTTTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACTTATTTATTTGCAAATTcc  <  1:167246/100‑1 (MQ=255)
                                                                                                   |                                                                                               
GTCAAGCGATGGGCGTTGCGCTCCATATTGTCTTACTTCCTTTTTTGAATTACTGCATAGCACAATTGATTCGTACGACGCCGACTTTGATGAGTCGGC‑‑TTTTTTTTGCCTGTTATTTATCAGCGTCTACCCTTTAAGAGTCCACCCAATGACCAGAGGGAAATATGACGACACTTATTTATTTGCAAATTCC  >  NC_000913/1409144‑1409336

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