Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 880,543 A→C intergenic (‑63/‑184) gstB ← / → dacC glutathione S‑transferase/D‑alanyl‑D‑alanine carboxypeptidase; penicillin‑binding protein 6a

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913880,5430AC99.1% 767.9 / ‑2.5 222intergenic (‑63/‑184)gstB/dacCglutathione S‑transferase/D‑alanyl‑D‑alanine carboxypeptidase; penicillin‑binding protein 6a
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (123/96);  total (125/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 222 total aligned reads displayed.

TTAGTTGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTAACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACGGTTTCT  >  NC_000913/880446‑880639
                                                                                                 |                                                                                                
 tAGTTGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACt                                                                                               <  1:2796686/100‑1 (MQ=255)
  aGTTGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACtt                                                                                              >  1:1594596/1‑100 (MQ=255)
   gTTGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACttt                                                                                             <  1:1389038/100‑1 (MQ=255)
   gTTGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACttt                                                                                             >  1:3245439/1‑100 (MQ=255)
   gTTGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACttt                                                                                             >  1:3052530/1‑100 (MQ=255)
    ttGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACtttt                                                                                            >  2:1490524/1‑100 (MQ=255)
     tGAATTATTCCGACCCCACAGCGTAGTCATACTCACCCCTTTTCCGAAAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTa                                                                                           >  1:4977343/1‑100 (MQ=255)
     tGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTa                                                                                           <  1:3948691/100‑1 (MQ=255)
     tGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTa                                                                                           >  2:3213726/1‑100 (MQ=255)
     tGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTa                                                                                           <  2:2275325/100‑1 (MQ=255)
      gAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTaa                                                                                          >  2:2382361/1‑100 (MQ=255)
       aaTTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAg                                                                                         >  2:2112204/1‑100 (MQ=255)
       aaTTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAg                                                                                         <  1:4211453/100‑1 (MQ=255)
       aaTTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAg                                                                                         >  2:3504828/1‑100 (MQ=255)
       aaTTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAg                                                                                         <  1:3212405/100‑1 (MQ=255)
         ttattCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCa                                                                                       <  1:1258292/100‑1 (MQ=255)
          tattCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCaa                                                                                      <  2:174540/100‑1 (MQ=255)
          tattCCGACCCCACAGCGTAATCATACCCACCCCTTTTCCGACAGGCAGCCCCAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCaa                                                                                      <  2:5243426/100‑1 (MQ=255)
           attCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAAc                                                                                     >  2:4093647/1‑100 (MQ=255)
            ttCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACa                                                                                    <  1:890888/100‑1 (MQ=255)
            ttCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACa                                                                                    >  1:1434356/1‑100 (MQ=255)
             tCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAg                                                                                   <  2:4706012/100‑1 (MQ=255)
             tCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAg                                                                                   >  1:2365998/1‑100 (MQ=255)
              ccGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGc                                                                                  >  2:456193/1‑100 (MQ=255)
                gACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTg                                                                                <  2:3931439/100‑1 (MQ=255)
                gACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTg                                                                                >  1:4085482/1‑100 (MQ=255)
                 aCCCCACAGCGTAATCATACTCTCCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTgg                                                                               <  1:3684012/100‑1 (MQ=255)
                 aCCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTgg                                                                               <  2:3379142/100‑1 (MQ=255)
                  ccccACAGCGTAATCATACTCGCCCCTTTTCCGACAGGCAGCCACGCGCGGAACAGGACGCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGc                                                                              <  2:5077810/100‑1 (MQ=255)
                  ccccACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGc                                                                              <  2:1575128/100‑1 (MQ=255)
                  ccccACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGc                                                                              <  2:2929844/100‑1 (MQ=255)
                     cacaGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCaaa                                                                           <  2:517605/100‑1 (MQ=255)
                          cGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAAt                                                                      <  2:2342148/100‑1 (MQ=255)
                          cGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAAt                                                                      <  1:4048440/100‑1 (MQ=255)
                           gTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATg                                                                     <  1:1988421/100‑1 (MQ=255)
                           gTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATg                                                                     <  2:3145796/100‑1 (MQ=255)
                             aaTCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGac                                                                   >  1:1508171/1‑100 (MQ=255)
                             aaTCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGac                                                                   >  1:987001/1‑100 (MQ=255)
                              aTCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGaca                                                                  >  1:811828/1‑100 (MQ=255)
                                cATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACAcc                                                                >  2:2188183/1‑100 (MQ=255)
                                cATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACAcc                                                                >  1:5240316/1‑100 (MQ=255)
                                   aCTCACCCATTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGAAAAAAAAGAAAAAAAAATGACACCAGa                                                             <  2:3433351/100‑1 (MQ=255)
                                    cTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGaa                                                            >  1:3875357/1‑100 (MQ=255)
                                    cTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGaa                                                            >  1:3305061/1‑100 (MQ=255)
                                    cTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGaa                                                            >  2:418062/1‑100 (MQ=255)
                                      cACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAagc                                                          >  1:2442846/1‑100 (MQ=255)
                                       aCCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAagca                                                         >  2:4741792/1‑100 (MQ=255)
                                         cccTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGCAGTAAAGTTCACTTTTAACCAACAGCTGGCAAAAAAATGACACCAGAagcagc                                                       <  1:5178828/100‑1 (MQ=255)
                                         cccTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAagcagc                                                       <  1:928891/100‑1 (MQ=255)
                                           cTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGa                                                     >  1:1546967/1‑100 (MQ=255)
                                               tCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCg                                                 <  1:3505131/100‑1 (MQ=255)
                                               tCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCg                                                 >  1:106819/1‑100 (MQ=255)
                                                    cAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAgg                                            >  2:557439/1‑100 (MQ=255)
                                                      ggCAGCCACAGGCTGAACAGGACTCCGATGGTGACGTAAAGTTCACTTTTACGCGACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGat                                          >  2:421570/1‑100 (MQ=255)
                                                      ggCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGat                                          >  2:2710225/1‑100 (MQ=255)
                                                      ggCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGat                                          >  2:3066841/1‑100 (MQ=255)
                                                      ggCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGat                                          <  2:875900/100‑1 (MQ=255)
                                                      ggCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGat                                          >  1:2959604/1‑100 (MQ=255)
                                                        cAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGatat                                        <  1:392011/100‑1 (MQ=255)
                                                         aGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATAtt                                       >  1:2540756/1‑100 (MQ=255)
                                                         aGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATAtt                                       >  1:367003/1‑100 (MQ=255)
                                                          gCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTg                                      >  1:1310206/1‑100 (MQ=255)
                                                          gCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTg                                      >  2:4423748/1‑100 (MQ=255)
                                                          gCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTg                                      <  1:3157970/100‑1 (MQ=255)
                                                          gCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTg                                      <  1:461247/100‑1 (MQ=255)
                                                           ccACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGc                                     <  1:116131/100‑1 (MQ=255)
                                                             acaGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCAt                                   >  1:1149670/1‑100 (MQ=255)
                                                             acaGACTGAACAGGACTCCGATGGTGAAGTAAAGGTCACTATTAAGCAACAGCTGGCAAACAAATGACACCAGAAACAGCCAAGCGGCAGGATATTGCAt                                   >  2:4645889/1‑100 (MQ=255)
                                                               aGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATaa                                 >  1:4808370/1‑100 (MQ=255)
                                                                ggCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATaaa                                >  1:3094857/1‑100 (MQ=255)
                                                                 gCTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAAc                               <  1:1530047/100‑1 (MQ=255)
                                                                  cTGAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACAAGAAGCACCGCAGCGGCACGATATTGCATAAACt                              >  2:1907299/1‑100 (MQ=255)
                                                                    gaccAGGACTACGATGGTGAAGTAAAGTTCACTTTCAAGCAACAACAAACAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACttt                            <  2:3812001/97‑1 (MQ=255)
                                                                    gAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACttt                            >  2:2122618/1‑100 (MQ=255)
                                                                    gAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACttt                            >  1:4886036/1‑100 (MQ=255)
                                                                    gAACAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACttt                            >  1:2679660/1‑100 (MQ=255)
                                                                     aaCAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTa                           >  2:648517/1‑100 (MQ=255)
                                                                       cAGGACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTaaa                         >  2:1264431/1‑100 (MQ=255)
                                                                         ggACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTaaaaa                       >  2:1968525/1‑100 (MQ=255)
                                                                         ggACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTaaaaa                       >  1:3429325/1‑100 (MQ=255)
                                                                         ggACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTaaaaa                       >  1:1805822/1‑100 (MQ=255)
                                                                         ggACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTaaaaa                       >  1:546416/1‑100 (MQ=255)
                                                                         ggACTCCGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTaaaaa                       >  1:789571/1‑100 (MQ=255)
                                                                              ccGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAACAACTTTa                  <  2:4825902/100‑1 (MQ=255)
                                                                              ccGATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTa                  >  1:1848827/1‑100 (MQ=255)
                                                                                gATGGTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTACAAACTTTAcc                <  2:255915/100‑1 (MQ=255)
                                                                                   ggTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAAc             <  2:4350223/100‑1 (MQ=255)
                                                                                   ggTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACt            >  1:4377552/1‑100 (MQ=255)
                                                                                   ggTGAAGTAAAGGTCACTTTTAAGCAACAGCTGGCAAATAAATGACACCAGAAGCAGCGAAGCGCCAGGATATTGCATAAACTTTAAACACTTTACGAAc             <  2:5357983/100‑1 (MQ=255)
                                                                                    gTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACt            <  1:588338/100‑1 (MQ=255)
                                                                                    gTGAAGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACt            >  2:3017682/1‑100 (MQ=255)
                                                                                      gAAGTAAAGTTCACTTTTAAGCAACAGCTGGAAAAAAAATGACACCCGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAGCTTa          >  2:5304818/1‑100 (MQ=255)
                                                                                       aaGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTAc         <  1:3343123/100‑1 (MQ=255)
                                                                                        aGTAAAGTTCATTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATACACTTTAAAAACTTTACCAACTTACg        <  2:4383400/100‑1 (MQ=255)
                                                                                        aGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACg        <  1:2713800/100‑1 (MQ=255)
                                                                                        aGTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACg        >  2:4965870/1‑100 (MQ=255)
                                                                                         gTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACgg       <  1:4538740/100‑1 (MQ=255)
                                                                                         gTAAAGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTAACAACTTACgg       >  2:2816224/1‑100 (MQ=255)
                                                                                           aaaGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACGGtt     <  2:1277783/100‑1 (MQ=255)
                                                                                            aaGTTCACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACGGttt    >  2:2663013/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TTAGTTGAATTATTCCGACCCCACAGCGTAATCATACTCACCCCTTTTCCGACAGGCAGCCACAGGCTGAACAGGACTCCGATGGTGAAGTAAAGTTAACTTTTAAGCAACAGCTGGCAAAAAAATGACACCAGAAGCAGCGAAGCGGCAGGATATTGCATAAACTTTAAAAACTTTACCAACTTACGGTTTCT  >  NC_000913/880446‑880639

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: