Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,095,945 C→A 5.9% intergenic (‑54/+35) hisQ ← / ← argT histidine/lysine/arginine/ornithine ABC transporter permease/lysine/arginine/ornithine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,095,9450CA 5.9% 220.9 / 3.6 68intergenic (‑54/+35)hisQ/argThistidine/lysine/arginine/ornithine ABC transporter permease/lysine/arginine/ornithine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (27/37);  new base A (2/2);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GGCCGAGCACGATGGCCAGGGCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAATGGCGTCGATCGCCGCGTTGA  >  NC_002947/5095798‑5096077
                                                                                                                                                   |                                                                                                                                    
ggCCGAGCACGATGGCCAGGGCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTctgct                                                                                                                                     >  4:24309/1‑149 (MQ=255)
  ccGAGCACGATGGCCAGGGCCATCGACGACTGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGt                                                                                                                                   <  4:246434/149‑1 (MQ=255)
  ccGAGCACGATGGCCAGGGCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGt                                                                                                                                   <  6:276411/149‑1 (MQ=255)
   cGAGCACGATGGCCAGGGCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGtt                                                                                                                                  <  3:24309/149‑1 (MQ=255)
    gAGCACGATGGCCAGGGCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGttt                                                                                                                                 <  3:190913/149‑1 (MQ=255)
      gCACGATGGCCAGGGCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGc                                                                                                                               <  1:177442/149‑1 (MQ=255)
        aCGATGGCCAGGGCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCAc                                                                                                                             <  3:77901/149‑1 (MQ=255)
         cGATGGCCAGGGCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACAGCAGAGGCTCTGAGGTTTGCAcc                                                                                                                            <  7:317352/149‑1 (MQ=39)
          gATGGCCAGGGCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGc                                                                                                                               <  6:165818/145‑1 (MQ=255)
          gATGGCCAGGGCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGc                                                                                                                               >  5:165818/1‑145 (MQ=255)
                  gggCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAg                                                                                                                   <  6:1768/149‑1 (MQ=255)
                  gggCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAg                                                                                                                   <  2:301019/149‑1 (MQ=255)
                    gCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGt                                                                                                                 <  1:362622/149‑1 (MQ=255)
                       aTCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGttt                                                                                                                                 >  6:186293/1‑130 (MQ=255)
                       aTCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGttt                                                                                                                                 <  5:186293/130‑1 (MQ=255)
                          gacgacAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAAc                                                                                                            <  3:78756/148‑1 (MQ=255)
                              acAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGttt                                                                                                                                 >  5:377812/1‑123 (MQ=255)
                              acAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGttt                                                                                                                                 <  6:377812/123‑1 (MQ=255)
                                  ggCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCt                                                                                                   >  2:316275/1‑149 (MQ=255)
                                                tCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTgg                                                                                     <  4:340090/149‑1 (MQ=255)
                                                  aGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGAc                                                                                   <  3:356006/149‑1 (MQ=255)
                                                     cACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTCCGAGTCTGGACCGt                                                                                <  7:320687/149‑1 (MQ=255)
                                                         ccccGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGAt                                                                            <  4:312431/149‑1 (MQ=255)
                                                            cGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTc                                                                         <  5:146642/149‑1 (MQ=255)
                                                            cGTCGAGGATGACTGCCCCGTAGCCTCTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTCCACACCTCAGAGCCTCTGCTGTTTGCCCCATTGCACAGGGGAAACCCGACGTCTTACGTCGAGTCTGGACCGTAGATGTc                                                                         >  5:102093/1‑149 (MQ=255)
                                                               cGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGa                                                                       >  2:43210/1‑148 (MQ=255)
                                                                gAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGaaa                                                                     >  3:112334/1‑149 (MQ=255)
                                                                  ggATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAAtt                                                                   >  1:98014/1‑149 (MQ=255)
                                                                     tGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGaa                                                                <  1:316275/149‑1 (MQ=255)
                                                                     tGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGaa                                                                <  4:112334/149‑1 (MQ=255)
                                                                        cTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGt                                                              >  6:207549/1‑148 (MQ=255)
                                                                              cGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCGAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCt                                                       >  5:330996/1‑149 (MQ=255)
                                                                                tAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCttt                                                     <  5:207549/149‑1 (MQ=255)
                                                                                     tttCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATg                                                                           <  1:51595/122‑1 (MQ=255)
                                                                                     tttCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATg                                                                           >  2:51595/1‑122 (MQ=255)
                                                                                      ttCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGAtt                                                >  4:122884/1‑148 (MQ=255)
                                                                                         aaCATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGAt                                                                            <  6:104153/117‑1 (MQ=255)
                                                                                          aCATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCttt                                           >  1:26078/1‑149 (MQ=255)
                                                                                               atgaGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACt                                      <  2:345170/149‑1 (MQ=255)
                                                                                               atgaGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACt                                      >  1:345170/1‑149 (MQ=255)
                                                                                                  aGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTg                                    <  3:122884/148‑1 (MQ=255)
                                                                                                    gTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCg                                 >  2:16563/1‑149 (MQ=255)
                                                                                                    gTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCg                                 <  2:98014/149‑1 (MQ=255)
                                                                                                    gTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCGTTGTACTTGCCg                                 >  7:296388/1‑149 (MQ=255)
                                                                                                       cTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTg                              >  2:81276/1‑149 (MQ=255)
                                                                                                       cTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACAGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTg                              >  6:79578/1‑149 (MQ=255)
                                                                                                         tctccTAGAGATTGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTg                              <  3:308794/144‑1 (MQ=255)
                                                                                                             cTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTg                              >  4:308794/1‑143 (MQ=255)
                                                                                                              tAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGaa                                                                >  4:263120/1‑108 (MQ=255)
                                                                                                              tAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGaa                                                                <  3:263120/108‑1 (MQ=255)
                                                                                                              tAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCt                                             >  7:273504/1‑127 (MQ=255)
                                                                                                              tAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCt                                             <  8:273504/127‑1 (MQ=255)
                                                                                                                gagaATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCttt                                           <  8:184375/127‑1 (MQ=255)
                                                                                                                gagaATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCttt                                           >  7:184375/1‑127 (MQ=255)
                                                                                                                   aaTGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAATg                     <  6:184862/146‑1 (MQ=255)
                                                                                                                   aaTGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAATg                     >  5:184862/1‑146 (MQ=255)
                                                                                                                     tGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAATGGCGTc                <  1:43210/149‑1 (MQ=255)
                                                                                                                     tGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGGCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAATGGCGTc                <  3:129715/149‑1 (MQ=255)
                                                                                                                               tGCTAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAATGGCGTCGATCGCcgcg      <  6:330996/149‑1 (MQ=255)
                                                                                                                               tGCAAACAGCAGAGGCTCTGAGGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACGGTAGATGTAGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGAACGAATGTCGTCCATCGCcgcg      >  5:263522/1‑149 (MQ=255)
                                                                                                                                gCAAACAGCAGAGGCTCTGAGGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCttt                                           <  5:216989/111‑1 (MQ=38)
                                                                                                                                gCAAACAGCAGAGGCTCTGAGGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCttt                                           >  6:216989/1‑111 (MQ=38)
                                                                                                                                 cAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTAGGCACGAATGGCGTCGATCGCCGCGtt    <  1:16563/149‑1 (MQ=255)
                                                                                                                                  aaaCCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTTCTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAATGGCGTCGATCGCCGCGtt    >  6:266923/1‑148 (MQ=255)
                                                                                                                                  aaaCCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAATGGCGt                 >  7:219508/1‑135 (MQ=255)
                                                                                                                                  aaaCCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAATGGCGt                 <  8:219508/135‑1 (MQ=255)
                                                                                                                                   aaCCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAATGGCGTCGATCGCCGCGTTGa  >  2:5942/1‑149 (MQ=255)
                                                                                                                                      cTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAAt                      >  4:271715/1‑126 (MQ=255)
                                                                                                                                      cTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAAt                      <  3:271715/126‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                    
GGCCGAGCACGATGGCCAGGGCCATCGACGACAGGGCGAGCTGTAGCGTCAGCCACGCCCCGTCGAGGATGACTGCCCCGTAGCCTTTCAACATGATGAGGTCCTCGACCTAGAGAATGAAAAATGGTGCAAACCTCAGAGCCTCTGCTGTTTGCACCATTGCACAGGTGAAACCCGACGTCTTACTTCGAGTCTGGACCGTAGATGTCGAAATTGAAGTACTTGCTTTCGATTTCTTTGTACTTGCCGTTGGCACGAATGGCGTCGATCGCCGCGTTGA  >  NC_002947/5095798‑5096077

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