Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,340,627 Δ1 bp 9.5% coding (244/1587 nt) ybiT ← ABC‑F family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,340,6270T. 9.5% 257.1 / 21.6 84coding (244/1587 nt)ybiTABC‑F family ATPase
Reads supporting (aligned to +/- strand):  ref base T (35/41);  new base . (3/5);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

GTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGATC  >  NC_002947/3340482‑3340775
                                                                                                                                                  |                                                                                                                                                    
gtgtAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACg                                                                                                                                                  >  2:947664/1‑150 (MQ=255)
gtgtAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACg                                                                                                                                                  <  1:630429/150‑1 (MQ=255)
    aGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGa                                                                                                                                              <  1:283476/150‑1 (MQ=255)
        gTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGatcat                                                                                                                                          <  1:2378913/150‑5 (MQ=255)
        gTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTc                                                                                                                                          <  2:2051856/150‑1 (MQ=255)
           cATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCATATCTGTCCGTTGCCCATGATCACCGTGTCGATCACGGTGAATTCCTc                                                                                                                                       <  1:2226054/150‑1 (MQ=255)
             tttCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGt                                                                                                                                     <  1:412206/150‑1 (MQ=255)
               tCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAg                                                                                                                                   <  2:2151510/150‑1 (MQ=255)
                 ggCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGc                                                                                                                                 >  2:836459/1‑150 (MQ=255)
                  gCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGa                                                                                                                                              <  1:2163255/136‑1 (MQ=255)
                  gCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGa                                                                                                                                              >  2:2163255/1‑136 (MQ=255)
                      aTTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACt                                                                                                                            >  2:2224555/1‑150 (MQ=255)
                       ttCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTg                                                                                                                           >  1:446712/1‑150 (MQ=255)
                       ttCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTg                                                                                                                           >  2:316885/1‑150 (MQ=255)
                       ttCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTCAGGAACTTGCCGGCCGACAGCAGCTGCGGTCACGCTCGGCTTTGACTTTCCACAGCTCGCCGTGGCCCATGCTCACCGTCCCGATCACGGTGAATTCCTCGTAGGCGACCTg                                                                                                                           >  2:84301/1‑150 (MQ=255)
                           gTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTg                                                                                                                           <  2:849166/146‑1 (MQ=255)
                           gTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTg                                                                                                                           >  1:849166/1‑146 (MQ=255)
                              tCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTgg                                                                                                                    <  1:1593830/150‑1 (MQ=255)
                                    tCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCt                                                                                                                                        <  1:1537486/124‑1 (MQ=255)
                                    tCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCt                                                                                                                                        >  2:1537486/1‑124 (MQ=255)
                                          aCGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTc                                                                                                                                       <  2:221185/119‑1 (MQ=255)
                                          aCGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTc                                                                                                                                       >  1:221185/1‑119 (MQ=255)
                                                tggCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCgg                                                                                                  <  2:1873413/148‑1 (MQ=255)
                                                aTGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCTCCCCATGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCgg                                                                                                  <  1:1620153/150‑1 (MQ=255)
                                                aTGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACg                                                                                                                                                  >  2:2336243/1‑102 (MQ=255)
                                                aTGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACg                                                                                                                                                  <  1:2336243/102‑1 (MQ=255)
                                                           ttcGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTc                                                                                       >  1:520738/1‑150 (MQ=255)
                                                              ggTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAg                                                                                    <  2:2340380/150‑1 (MQ=255)
                                                               gTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAg                                                                                                                                   >  1:733886/1‑102 (MQ=255)
                                                               gTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAg                                                                                                                                   <  2:733886/102‑1 (MQ=255)
                                                                     tCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGAcc                                                                             <  2:190843/150‑1 (MQ=255)
                                                                      ccGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCt                                                                            <  1:2282371/150‑1 (MQ=255)
                                                                        ggCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTgg                                                                          >  1:1636088/1‑150 (MQ=255)
                                                                         gCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAATGTGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                         <  1:173583/150‑1 (MQ=255)
                                                                           aGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCt                                                                       >  2:1887465/1‑150 (MQ=255)
                                                                              gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGt                                                                                                                                                >  1:125/1‑74 (MQ=255)
                                                                              gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGt                                                                                                                                                <  2:125/74‑1 (MQ=255)
                                                                              gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCg                                                                    <  1:947664/150‑1 (MQ=255)
                                                                                gTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGa                                                                  <  1:836459/150‑1 (MQ=255)
                                                                                gTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGa                                                                  <  1:1887465/150‑1 (MQ=255)
                                                                                gTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTACTCGGCTCGAGCATGACCTGGCCTGCGGa                                                                  <  1:1647033/150‑1 (MQ=255)
                                                                                 tAGATACGGTCACGTTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCTTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGAc                                                                 <  2:82431/150‑1 (MQ=255)
                                                                                 tAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                         <  1:1167065/142‑1 (MQ=255)
                                                                                 tAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                         >  2:1167065/1‑142 (MQ=255)
                                                                                 tAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGAc                                                                 >  2:2052371/1‑150 (MQ=255)
                                                                                 tAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGAc                                                                 <  1:554510/150‑1 (MQ=255)
                                                                                              gCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTAtt                                                                                             <  1:1413859/109‑1 (MQ=255)
                                                                                              gCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTAtt                                                                                             >  2:1413859/1‑109 (MQ=255)
                                                                                              gCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCAc                                                    >  1:1206638/1‑150 (MQ=255)
                                                                                                tCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCAccg                                                  >  1:1607111/1‑150 (MQ=255)
                                                                                                    cTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGa                                              <  2:1165267/150‑1 (MQ=255)
                                                                                                    cTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGa                                              <  2:1867271/150‑1 (MQ=255)
                                                                                                           tttCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGAtttt                                       >  1:1787311/1‑150 (MQ=255)
                                                                                                           tttCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGAtttt                                       >  1:2321120/1‑150 (MQ=255)
                                                                                                            ttCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTc                                      <  1:383043/150‑1 (MQ=255)
                                                                                                              ccACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCAt                                    <  1:2052371/150‑1 (MQ=255)
                                                                                                                acaGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGa                                  >  2:1809538/1‑150 (MQ=255)
                                                                                                                     tGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGAGGGTATTCGGCTCGAGCATGAGCTCGACTTGGGACTGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTc                             <  1:2224555/150‑1 (MQ=255)
                                                                                                                           tGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGATCTGGCCTGCGGACGGCGCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTg                       <  2:446712/150‑1 (MQ=255)
                                                                                                                               ccATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGc                   >  2:2120/1‑150 (MQ=255)
                                                                                                                                cATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCgg                                                                                                  <  1:533642/70‑1 (MQ=255)
                                                                                                                                cATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCgg                                                                                                  >  2:533642/1‑70 (MQ=255)
                                                                                                                                 aTGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCaa                 <  1:2120/150‑1 (MQ=255)
                                                                                                                                 aTGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCaa                 >  2:1158981/1‑150 (MQ=255)
                                                                                                                                   gATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGGCGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAAcc               >  2:545104/1‑150 (MQ=255)
                                                                                                                                   gATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCACAGCTTGCCCAGGCGGGTCTACGGCTCGAGCATGACCTGGCCTGCGAACGGCTCCAAGTCACCGCCGAGGATTTTCATGAAGGTCGCCTTGCCGCAAcc               >  2:345244/1‑150 (MQ=255)
                                                                                                                                   attCACCGTGATCGACACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAg                                                         <  2:1473352/106‑1 (MQ=255)
                                                                                                                                   attCACCGTGATCGACACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAg                                                         >  1:1473352/3‑108 (MQ=255)
                                                                                                                                    ttCACCGTGATCGACACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGa                                                                                     >  1:738969/2‑79 (MQ=38)
                                                                                                                                    ttCACCGTGATCGACACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGa                                                                                     <  2:738969/78‑1 (MQ=38)
                                                                                                                                    aTCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGAc                                                                 >  2:1252623/1‑99 (MQ=255)
                                                                                                                                    aTCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGAc                                                                 <  1:1252623/99‑1 (MQ=255)
                                                                                                                                    aTCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCg              >  1:259742/1‑150 (MQ=255)
                                                                                                                                     tCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGTTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGGCGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGCCTTGCCGCAACCGt             >  1:2295951/1‑150 (MQ=255)
                                                                                                                                      cACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGtt            <  2:1787311/150‑1 (MQ=255)
                                                                                                                                      cACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGtt            >  1:2250741/1‑150 (MQ=255)
                                                                                                                                        ccGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCACAAATCGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTgg          <  2:676200/150‑1 (MQ=255)
                                                                                                                                         cGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGgc         <  2:711513/150‑1 (MQ=255)
                                                                                                                                           tgtCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGgcgc       <  2:259742/150‑1 (MQ=255)
                                                                                                                                           tgtCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGgcgc       >  1:81254/1‑150 (MQ=255)
                                                                                                                                               cGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGa    <  1:1158981/150‑1 (MQ=255)
                                                                                                                                               cGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGa    <  1:316885/150‑1 (MQ=255)
                                                                                                                                                 aTCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGATc  >  1:1340383/1‑150 (MQ=255)
                                                                                                                                                 aTCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGATc  >  2:114772/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCCGATC  >  NC_002947/3340482‑3340775

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