Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,670 C→A 15.7% intergenic (+315/‑40) rnlB → / → yfjP CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,6700CA15.7% 441.3 / 65.1 171intergenic (+315/‑40)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base C (71/73);  new base A (13/14);  total (84/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

Only 100 of 175 total aligned reads displayed.

GTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGG  >  NC_000913/2767571‑2767769
                                                                                                   |                                                                                                   
gTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgc                                                                                                     >  1:2399274/1‑100 (MQ=255)
 tGCTGTGAGTTCAGAAATATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgcg                                                                                                    <  2:1919763/100‑1 (MQ=255)
  gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGGACCTCGCgg                                                                                                   <  2:1449569/100‑1 (MQ=255)
  gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                                   <  1:3197080/100‑1 (MQ=255)
  gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                                   >  1:2340272/1‑100 (MQ=255)
  gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCgg                                                                                                   >  2:3004875/1‑100 (MQ=255)
    tgtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                                 >  1:1864905/1‑100 (MQ=255)
     gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                                <  2:3660142/100‑1 (MQ=255)
     gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTc                                                                                                >  1:1933347/1‑100 (MQ=255)
       gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                              <  1:585297/100‑1 (MQ=255)
       gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                              >  1:3485329/1‑100 (MQ=255)
       gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTg                                                                                              >  2:2719372/1‑100 (MQ=255)
        aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGc                                                                                             <  1:2795194/100‑1 (MQ=255)
        aGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGc                                                                                             <  1:3170946/100‑1 (MQ=255)
         gTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGcc                                                                                            >  2:1537294/1‑100 (MQ=255)
          ttCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                           <  1:1084212/100‑1 (MQ=255)
           tCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGt                                                                                          <  1:655068/100‑1 (MQ=255)
            cAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtt                                                                                         >  2:1069998/1‑100 (MQ=255)
             aGAAAGATTCTGCAATCTATGTTCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCg                                                                                           <  2:262638/97‑1 (MQ=39)
             aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGttt                                                                                        >  2:393398/1‑100 (MQ=255)
             aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttc                                                                                        <  1:2602368/100‑2 (MQ=255)
             aGAAAGATTCTGCAATCTAAGAAAAATAAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGtt                                                                                         <  2:2835546/99‑1 (MQ=37)
              gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                       >  1:2254025/1‑100 (MQ=255)
              gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                       >  2:101941/1‑100 (MQ=255)
              gAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGtttt                                                                                       >  1:1523478/1‑100 (MQ=255)
                aataTTATGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCGCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTGCTGCCggtcatc                                                                                     <  2:606517/97‑7 (MQ=37)
                   aTTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGTACCTCGCGGTTCTGCCGTTTTTGttt                                                                                  <  2:3788203/100‑1 (MQ=255)
                   aTTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGtt                                                                                   >  2:679016/1‑99 (MQ=255)
                   aTTCTGCAATATATTTGCATTAAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGttt                                                                                  <  2:482712/100‑1 (MQ=255)
                    ttCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGtttt                                                                                 <  1:1683253/100‑1 (MQ=255)
                      cTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                               >  1:3330457/1‑100 (MQ=255)
                       tGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                               >  2:4096995/1‑99 (MQ=255)
                       tGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACt                                                                              >  1:3395877/1‑100 (MQ=255)
                          aaTCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACtttt                                                                           >  1:2784020/1‑100 (MQ=255)
                             cTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCag                                                                        <  1:3863247/100‑1 (MQ=255)
                             cTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCag                                                                        <  1:4154649/100‑1 (MQ=255)
                              tATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCaga                                                                       <  2:813492/100‑1 (MQ=255)
                               aTGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCagag                                                                      >  1:2212852/1‑100 (MQ=255)
                                tgtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                     >  1:2511875/1‑100 (MQ=255)
                                 gtgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagg                                                                     <  2:3469072/99‑1 (MQ=255)
                                  tgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGATCCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgaggag                                                                   <  2:1068171/100‑1 (MQ=255)
                                   gCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgt                                                                  >  1:761012/1‑100 (MQ=255)
                                     aataaaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgt                                                                  <  2:3321214/94‑1 (MQ=255)
                                      ttGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggt                                                               <  1:1096006/100‑1 (MQ=255)
                                      ttGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggt                                                               >  1:2346964/1‑100 (MQ=255)
                                      ttGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAAACCGATGAACGGCAGAACCTAGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGTAgtgg                                                                >  2:1648414/1‑99 (MQ=255)
                                       tGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggtg                                                              >  1:1993350/1‑100 (MQ=255)
                                        gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                             >  1:137500/1‑100 (MQ=255)
                                         aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                            <  1:1621287/100‑1 (MQ=255)
                                         aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCACTGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                            >  2:3937180/1‑100 (MQ=255)
                                              aTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGaaaa                                                       <  1:124314/100‑1 (MQ=255)
                                                cTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAt                                                     >  2:879028/1‑100 (MQ=255)
                                                   cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                  <  1:1568982/100‑1 (MQ=255)
                                                   cAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                  <  2:2391813/100‑1 (MQ=255)
                                                     gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGa                                                <  2:3209923/100‑1 (MQ=255)
                                                     gatagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGATGATGAAAAATTTTGa                                                >  2:1248488/1‑100 (MQ=255)
                                                       tagataTCATTCTACTGAACAATCCGAAGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAg                                              <  1:3769518/100‑1 (MQ=255)
                                                         gataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGtt                                            >  1:2003773/1‑100 (MQ=255)
                                                         gataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGtt                                            >  1:464209/1‑100 (MQ=255)
                                                         gataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGTGTGGTGATGAAAACTTTTGAAGtt                                            >  1:1857000/1‑100 (MQ=255)
                                                          atatCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTc                                           <  1:607840/100‑1 (MQ=255)
                                                            atCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTg                                         <  1:2884380/100‑1 (MQ=255)
                                                              cATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCa                                       <  1:3940051/100‑1 (MQ=255)
                                                               aTTCTCCTGAACAATCCGATGCACGGCAGAACCTCGCGGTTCACCCGTCTTTGTTTTACTTCACAAACGAGCGGTGGGGAAGAATTGTGCAgacgtccag                                      >  2:3982804/1‑92 (MQ=255)
                                                                    ccTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACt                                 >  2:3906254/1‑100 (MQ=255)
                                                                     cTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTa                                >  1:842129/1‑100 (MQ=255)
                                                                      tGAACAATCCGCTGAACGACAGAACCTCGCGGTTCTGTCTTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTAc                               >  2:1806610/1‑100 (MQ=255)
                                                                      tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTAc                               <  2:204427/100‑1 (MQ=255)
                                                                       gAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTTAAGTTCTGCAGCCACTACa                              >  2:3542487/1‑100 (MQ=255)
                                                                        aaCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAg                             <  2:1846398/100‑1 (MQ=255)
                                                                        aaCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAg                             <  1:3004888/100‑1 (MQ=255)
                                                                          cAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGaa                           <  1:3378424/100‑1 (MQ=255)
                                                                           aaTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAAc                          <  2:3648992/100‑1 (MQ=255)
                                                                            aTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTTCTTTTCATAGGAGAGGAGGTGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACt                         >  2:22675/1‑100 (MQ=255)
                                                                            aTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGCGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGacc                          >  2:795769/1‑97 (MQ=255)
                                                                            aTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACt                         >  1:2403653/1‑100 (MQ=255)
                                                                            aTCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGAGGTGATGGAAAATTTTGAAGTTCTGCAGCCACTAc                               >  2:3420056/1‑94 (MQ=255)
                                                                              ccGATGATCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTc                        <  2:1358501/99‑1 (MQ=255)
                                                                               cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                      <  1:1481843/100‑1 (MQ=255)
                                                                               cGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTccc                      >  2:1192733/1‑100 (MQ=255)
                                                                                gATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                     >  1:3260340/1‑100 (MQ=255)
                                                                                caaaaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                     >  2:2124002/5‑100 (MQ=38)
                                                                                   aaaCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTc                  >  1:1832824/2‑100 (MQ=255)
                                                                                    aaCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTTCTTTTCAGAGGAGTAGGGATGAAATATTTTGTAGTTCTGTAGCCAATACAGAACTCACttt                   >  2:846518/1‑98 (MQ=255)
                                                                                     aCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCg                <  1:800649/100‑1 (MQ=255)
                                                                                     aCGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACCCCCTTTCCg                >  2:2884540/1‑100 (MQ=255)
                                                                                      cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg               <  2:2534470/100‑1 (MQ=255)
                                                                                      cGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg               <  1:1313738/100‑1 (MQ=255)
                                                                                      cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCgg               >  1:3722619/1‑100 (MQ=255)
                                                                                      cGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGGGGAGATGAAAAATTTTGAAGTTGTGCAGCCACTAAAGATCTCCCTTTc                  >  2:2037531/1‑97 (MQ=37)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt              >  1:783179/1‑100 (MQ=255)
                                                                                       ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt              >  1:1224602/1‑100 (MQ=255)
                                                                                        gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc             >  1:937585/1‑100 (MQ=255)
                                                                                        gCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc             <  1:4118251/100‑1 (MQ=255)
                                                                                        gCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc             >  1:1385545/1‑100 (MQ=255)
                                                                                         cAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCt            >  1:3666029/1‑100 (MQ=255)
                                                                                         cAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCt            >  2:799514/1‑100 (MQ=255)
                                                                                           gAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTAc          >  2:2896573/1‑100 (MQ=255)
                                                                                                tCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCt      >  1:2646998/1‑99 (MQ=255)
                                                                                                tCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTa     <  1:2555867/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGG  >  NC_000913/2767571‑2767769

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: