Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,700,758 T→C 8.2% intergenic (‑71/+79) kdpD ← / ← kdpC sensor protein/potassium‑transporting ATPase subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,700,7580TC 8.2% 175.4 / 5.0 62intergenic (‑71/+79)kdpD/kdpCsensor protein/potassium‑transporting ATPase subunit C
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base C (3/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

CACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGAT  >  NC_002947/4700616‑4700905
                                                                                                                                              |                                                                                                                                                   
cacCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCcgc                                                                                                                                               <  4:302712/149‑1 (MQ=255)
    ttCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                            >  7:152313/1‑148 (MQ=255)
       aGCCGGCCCCGGCCGGTTCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATAc                                                                                                                                         <  7:337968/148‑1 (MQ=255)
         ccggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGt                                                                                                                                       >  7:269729/1‑148 (MQ=255)
         ccggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGt                                                                                                                                       >  4:339160/1‑148 (MQ=255)
          cggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCa                                                                                                                                     <  8:152313/149‑1 (MQ=255)
           ggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAg                                                                                                                                    <  5:204276/149‑1 (MQ=255)
           ggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAg                                                                                                                                    <  5:209275/149‑1 (MQ=255)
               ccggccggGGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCCGGCGGAGTCACTCATGTGTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCCGGGAAATCCCCCATAACGTGTCGGAGCCGCTTAACCCGCGATTACGTCAGTgg                                                                                                                                 >  4:281388/1‑148 (MQ=255)
                cggccggTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGca                                                                                                                               >  4:93985/1‑149 (MQ=255)
                               aGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTcg                                                                                                                >  7:18555/1‑149 (MQ=255)
                                  ttGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCgg                                                                                                             <  6:17151/149‑1 (MQ=255)
                                  ttGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCgg                                                                                                             <  1:125622/149‑1 (MQ=255)
                                  ttGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCgg                                                                                                             <  1:209051/149‑1 (MQ=25)
                                    ggCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGt                                                                                                           >  7:201695/1‑149 (MQ=255)
                                     gCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGCGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGt                                                                                                           <  3:20359/148‑1 (MQ=255)
                                      ccAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTg                                                                                                          >  3:285947/1‑148 (MQ=255)
                                       cAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGaa                                                                                                        >  4:267954/1‑149 (MQ=255)
                                        aaCAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAAc                                                                                                       >  1:366860/1‑149 (MQ=255)
                                         aCAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAAcc                                                                                                      <  8:18555/149‑1 (MQ=255)
                                             cgcATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCgg                                                                                                             >  6:29472/1‑138 (MQ=21)
                                             cgcATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCgg                                                                                                             <  5:29472/138‑1 (MQ=21)
                                                aTCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTcc                                                                                                >  3:257006/1‑148 (MQ=255)
                                                    gcgcGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAg                                                                                           <  7:354793/149‑1 (MQ=255)
                                                           cggAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCt                                                                                    >  6:36742/1‑149 (MQ=255)
                                                           cggAGTCACTCATGGTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCt                                                                                    <  7:129375/149‑1 (MQ=255)
                                                               gTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCa                                                                                <  8:269729/149‑1 (MQ=255)
                                                                  aCTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAgcg                                                                             >  5:203465/1‑149 (MQ=255)
                                                                    tCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTgg                                                                                                                                 <  5:5903/95‑1 (MQ=255)
                                                                    tCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTgg                                                                                                                                 >  6:5903/1‑95 (MQ=255)
                                                                         tttCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGAt                                                                      <  6:58223/149‑1 (MQ=255)
                                                                            cATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTAc                                                                   >  8:279219/1‑149 (MQ=255)
                                                                             aTCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGt                                                                                                                                       >  5:290502/1‑80 (MQ=255)
                                                                             aTCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGt                                                                                                                                       <  6:290502/80‑1 (MQ=255)
                                                                             aTCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTAc                                                                   <  6:338040/148‑1 (MQ=255)
                                                                                cTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGa                                                               >  2:279213/1‑149 (MQ=255)
                                                                                 tttACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGa                                                               >  5:69442/1‑148 (MQ=255)
                                                                                     cTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATg                                                          >  8:192141/1‑149 (MQ=255)
                                                                                        tGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGc                                                        <  6:203465/148‑1 (MQ=255)
                                                                                          tCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAg                                                     >  4:352598/1‑149 (MQ=255)
                                                                                            cAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAAt                                                           <  3:310556/141‑1 (MQ=255)
                                                                                            cAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAAt                                                           >  4:310556/1‑141 (MQ=255)
                                                                                             aaCAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCtt                                                  >  3:84012/1‑149 (MQ=255)
                                                                                                     gccgcTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAgg                                          >  7:230940/1‑149 (MQ=255)
                                                                                                     gccgcTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAgg                                          >  8:83998/1‑149 (MQ=255)
                                                                                                           tCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGt                                    >  6:141284/1‑149 (MQ=255)
                                                                                                                gAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAAc                               <  4:84012/149‑1 (MQ=255)
                                                                                                                    gCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCa                           <  8:230940/149‑1 (MQ=255)
                                                                                                                    gCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCa                           <  7:279219/149‑1 (MQ=255)
                                                                                                                    gCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCa                           <  3:93985/149‑1 (MQ=255)
                                                                                                                      ccccTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAg                          >  6:76270/1‑148 (MQ=255)
                                                                                                                       cccTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCa                        >  6:251669/1‑149 (MQ=255)
                                                                                                                           tAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTcc                                                                                                <  4:256707/73‑1 (MQ=255)
                                                                                                                           tAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTcc                                                                                                >  3:256707/1‑73 (MQ=255)
                                                                                                                            aaCTTGTAGGTGGGGTTTTACCCGCGAATACGTCAGTGGCACCACCCCCGCATTCGCGGGTGAACCCGGTCCTACAGGCTTCCTGCCAGCTCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCTCGtt                    <  8:196495/148‑1 (MQ=38)
                                                                                                                             aCGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTg                   <  8:201695/148‑1 (MQ=255)
                                                                                                                             aCGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGa                  <  1:15943/149‑1 (MQ=255)
                                                                                                                             aCGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGa                  <  3:267954/149‑1 (MQ=255)
                                                                                                                                    ggAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcgg           >  6:224697/1‑149 (MQ=255)
                                                                                                                                     gAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTAGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcggc          >  6:60126/1‑149 (MQ=255)
                                                                                                                                          ggTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCGGGTAAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGcc     >  5:208515/1‑149 (MQ=21)
                                                                                                                                           gTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCGGGTAAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCg    >  2:366383/1‑149 (MQ=21)
                                                                                                                                            tttACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGa   >  2:60970/1‑149 (MQ=255)
                                                                                                                                              tACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGAt  <  7:192141/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGAT  >  NC_002947/4700616‑4700905

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: