Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 186,789 G→A 6.5% A43V (GCC→GTC)  PP_0163 ← transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947186,7890GA 6.5% 182.1 / 2.4 64A43V (GCC→GTC) PP_0163transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base A (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

ACTGCATCCTCGCCGATCAGCGGTGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCTTGGCCACGCGACGCCTGCCTCAACGAT  >  NC_002947/186642‑186936
                                                                                                                                                   |                                                                                                                                                   
aCTGCATCCACGCCGATCAGCGGTGCAAAGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACgg                                                                                                                                                     >  5:307837/1‑148 (MQ=255)
   gCATCCTCGCCGATCAGCGGTGCAACGTCGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAg                                                                                                                                                  >  6:357580/1‑148 (MQ=255)
     aTCCTCGCCGATCAGCGGTGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTAc                                                                                                                                               <  8:233154/149‑1 (MQ=255)
        cTCGCCGATCAGCGGTGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGGGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAg                                                                                                                                            >  5:292117/1‑149 (MQ=255)
        cTCGCCGATCAGCGGTGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAg                                                                                                                                            >  1:173524/1‑149 (MQ=255)
              gATCAGCGGTGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGGGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGc                                                                                                                                       >  7:156825/1‑148 (MQ=255)
               aTCAGCGGTGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGc                                                                                                                                       >  5:272549/1‑147 (MQ=255)
               aTCAGCGGTGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCaa                                                                                                                                     <  1:276476/149‑1 (MQ=255)
                  aGCGGTGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCg                                                                                                                                  >  4:39123/1‑149 (MQ=255)
                  aGCGGTGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCg                                                                                                                                  <  5:344833/149‑1 (MQ=255)
                  aGCGGTGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGATCTCGAACGGCAGTAACAGGGGGCaacc                                                                                                                                   >  7:337229/1‑146 (MQ=255)
                      gTGCAACGACGGCGCCCTTGTGCATCTCCACTGTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGCCGCCGTACAGGGTGGCGCGCTCTTGTTGGCGCACGGCAGTACCAGGGGGCACGCGACCa                                                                                                                              >  1:274865/1‑149 (MQ=255)
                       tGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCAt                                                                                                                             >  8:149878/1‑149 (MQ=255)
                        gCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGTAGTACCAGGGGGCAAGCGACCATc                                                                                                                            <  6:272549/149‑1 (MQ=255)
                         cAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCg                                                                                                                           <  4:165163/149‑1 (MQ=255)
                         cAACGACGGCGCCCTAGTGCATCTCCACTTTCACCAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGACATGCGACAGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGTCAGTACCAGGGTGCAAGCGACCATc                                                                                                                            <  6:181435/148‑1 (MQ=255)
                          aaCGACGTGGCCCTTGTGCAGCTCCACTTTCAACAGGGATTGCGGTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGTGAGTACTAGGGTGCAAGCGACCATCga                                                                                                                          <  6:307837/149‑1 (MQ=255)
                          aaCGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGGGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCga                                                                                                                          <  6:292117/149‑1 (MQ=255)
                          aaCGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTGTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATc                                                                                                                            >  2:324587/1‑147 (MQ=255)
                          aaCGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATc                                                                                                                            >  3:170883/1‑147 (MQ=255)
                           acgacgGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATc                                                                                                                            >  7:143498/1‑146 (MQ=255)
                                  cgcCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGgcg                                                                                                                  >  2:366434/1‑149 (MQ=255)
                                    cccTTTTGAATCTCACCTTTCAACAGCGATTGCGCATCGAGCTGGCGCAGGGCTTCGCGAACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCt                                                                                                                <  8:156825/149‑1 (MQ=255)
                                     ccTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTg                                                                                                               <  7:88185/149‑1 (MQ=255)
                                      cTTGTGCAGCTCAACTTTCAACAGGGAGTGCGCTTCGCTCTGGCGCAGGGCTGCGCGCACGGGCATTCGACTGACGCCGTACAGGGTGGCGAGATCTGGTTGGCGAACGGCAGTACCAGTGGGCAAGCGACCATCGATAATGGCGCTgc                                                                                                              <  2:274865/149‑1 (MQ=255)
                                             aTCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCg                                                                                                       >  7:286490/1‑149 (MQ=255)
                                             aTCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCg                                                                                                       >  8:257990/1‑149 (MQ=255)
                                                  cACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCATGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGttctt                                                                                                  <  8:163888/149‑1 (MQ=255)
                                                  cACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGttct                                                                                                   <  8:143498/148‑1 (MQ=255)
                                                  cACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCAAGAATGGCGCTGCGCAGGCGttctt                                                                                                  <  7:149878/149‑1 (MQ=255)
                                                   aCTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGttcttc                                                                                                 <  1:260684/149‑1 (MQ=255)
                                                   aCTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGttcttc                                                                                                 <  2:49623/149‑1 (MQ=255)
                                                   aCTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGttcttc                                                                                                 >  1:49623/1‑149 (MQ=255)
                                                     tttCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCACACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCATGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCAACGAGAATGGCGCTGCGCAGGCGTTCTTCg                                                                                                >  7:57756/1‑148 (MQ=255)
                                                     attCAACAGGGATTGCGCATCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTCACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGAACCAGGCGCCAAGCGACCATCGAGAAAGGCGCAGCGCAGGCGatcttcca                                                                                               >  4:355399/2‑147 (MQ=255)
                                                         aaCAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCga                                                                                                                          <  8:96450/118‑1 (MQ=255)
                                                         aaCAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATcac                                                                                           <  5:357580/149‑1 (MQ=255)
                                                            aGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATcacgg                                                                                         >  2:28985/1‑148 (MQ=255)
                                                                aTTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATcacggcacgg                                                                                    <  1:366434/149‑1 (MQ=255)
                                                                 ttGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACggg                                                                                   <  7:257990/149‑1 (MQ=255)
                                                                   gcgcTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGc                                                                                  <  2:173524/148‑1 (MQ=255)
                                                                   gcgcTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGACAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGcc                                                                                 >  6:151316/1‑149 (MQ=255)
                                                                   gcgcTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGACAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGcc                                                                                 >  8:130996/1‑149 (MQ=255)
                                                                         cGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGAt                                                                             <  3:39123/147‑1 (MQ=255)
                                                                         cGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGAt                                                                             <  8:286490/147‑1 (MQ=255)
                                                                                  gctgggCTTCGGGCACGGGGATGGGACTGACGCCGTAAAGGGTGTCGCGCTCTTGTTGGCGAACGTCAGTACCAGGGGGCAAGAGACCAGCGAGAATGTCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGGCTGATGCGGCGGCAAc                                                                  <  8:337229/146‑1 (MQ=255)
                                                                                  gcAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAAc                                                                  >  7:153451/1‑149 (MQ=255)
                                                                                  gcAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAAc                                                                  >  7:291643/1‑149 (MQ=255)
                                                                                               cACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCa                                                     <  4:170883/149‑1 (MQ=255)
                                                                                                  gggCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGca                                                  <  8:153451/149‑1 (MQ=255)
                                                                                                             tGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAAt                                       >  4:4961/1‑149 (MQ=255)
                                                                                                                 gCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGACAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGa                                   <  7:130996/149‑1 (MQ=255)
                                                                                                                 gCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGACAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGa                                   <  5:151316/149‑1 (MQ=255)
                                                                                                                    gTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGAttt                                >  3:100403/1‑149 (MQ=255)
                                                                                                                     tACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGAttt                                >  8:127018/1‑148 (MQ=255)
                                                                                                                        aGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCtt                            <  8:291643/149‑1 (MQ=255)
                                                                                                                          ggTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGCCCAGCACATTGCTCAATTTGATTTTCTTgg                          <  1:28985/149‑1 (MQ=255)
                                                                                                                           gTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCTTGGc                         >  4:186720/1‑149 (MQ=255)
                                                                                                                                gAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCTTGGCCAcgc                    >  6:158979/1‑149 (MQ=255)
                                                                                                                                         ttgGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCTTGGCCACGCGACgcctgc           >  3:311276/1‑149 (MQ=255)
                                                                                                                                         ttgGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCTTGGCCACGCGACgcctgc           >  4:31415/1‑149 (MQ=255)
                                                                                                                                           ggCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCTTGGCCACGCGACgcctgcc          <  3:174466/148‑1 (MQ=255)
                                                                                                                                              gAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCTTGGCCACGCGACGCCTGCCTCaa      >  2:95945/1‑149 (MQ=255)
                                                                                                                                               aaCGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCTTGGCCACGCGACGCCTGCCTCAAc     <  3:186720/149‑1 (MQ=255)
                                                                                                                                               aaCGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCTTGGCCACGCGACGCCTGCCTCAAc     >  7:39108/1‑149 (MQ=255)
                                                                                                                                                  ggCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCTTGGCCACGCGACGCCTGCCTCAACGAt  <  1:95945/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACTGCATCCTCGCCGATCAGCGGTGCAACGACGGCGCCCTTGTGCATCTCCACTTTCAACAGGGATTGCGCTTCGAGCTGGCGCAGGGCTTCGCGCACGGGCATGCGACTGACGCCGTACAGGGTGGCGAGCTCTTGTTGGCGAACGGCAGTACCAGGGGGCAAGCGACCATCGAGAATGGCGCTGCGCAGGCGTTCTTCGATCACGGCACGGGCCTGATGCGGCGGCAACTGCTCGCTGGCCAGCACATTGCTCAATTTGATTTTCTTGGCCACGCGACGCCTGCCTCAACGAT  >  NC_002947/186642‑186936

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: ‑