Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,957,528 T→A 7.9% L311H (CTC→CAC) ‡ cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,957,5280TA 7.9% 191.7 / 5.6 63L311H (CTC→CAC) ‡cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base A (2/3);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CTGTGGCAATTGTTCAAAGGCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGCCATGGC  >  NC_002947/5957383‑5957676
                                                                                                                                                 |                                                                                                                                                     
cTGTGGCAATTGTTCAAAGGCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGctcgacctcga                                                                                                                                                   >  7:141385/1‑149 (MQ=255)
   tgGCAATTGTTCAAAGGCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTcgacga                                                                                                                                                <  1:144432/149‑1 (MQ=255)
      cAATTGTTCAAAGGCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCt                                                                                                                                             <  8:141385/149‑1 (MQ=255)
         ttGTTCAAAGGCATCGACTCTCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGa                                                                                                                                          <  8:200449/149‑1 (MQ=255)
         ttGTTCAAAGGCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGTCGAGGTCGAGCTggt                                                                                                                                          <  6:350618/149‑2 (MQ=255)
             tCAAAGGCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGAcccc                                                                                                                                      >  4:35206/1‑149 (MQ=255)
              cAAAGGCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGAc                                                                                                                                         <  4:250032/145‑1 (MQ=255)
              cAAAGGCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGAc                                                                                                                                         >  3:250032/1‑145 (MQ=255)
                  ggCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGt                                                                                                                                 <  1:208239/149‑1 (MQ=255)
                     aTCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGAt                                                                                                                              >  5:92547/1‑149 (MQ=255)
                                cTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGgc                                                                                                                   >  3:257265/1‑149 (MQ=255)
                                 tACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGgcg                                                                                                                  >  1:150552/1‑149 (MQ=255)
                                     aaTCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGAt                                                                                                                              <  5:155823/133‑1 (MQ=255)
                                     aaTCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGAt                                                                                                                              >  6:155823/1‑133 (MQ=255)
                                         gctgctCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATg                                                                                                                             >  6:160140/1‑130 (MQ=255)
                                         gctgctCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATg                                                                                                                             <  5:160140/130‑1 (MQ=255)
                                          ctgctCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACg                                                                                                         <  4:148614/149‑1 (MQ=255)
                                             ctCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGt                                                                                                        <  5:28189/147‑1 (MQ=255)
                                              tCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGTCGAGGTCGAGCTggt                                                                                                                                          <  8:64378/112‑2 (MQ=255)
                                              tCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGTCGAGGTCGAGCTggt                                                                                                                                          >  7:64378/1‑111 (MQ=255)
                                                aaaCTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGctcgacctc                                                                                                                                                      <  2:208619/99‑1 (MQ=255)
                                                aaaCTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGctcgacctc                                                                                                                                                      >  1:208619/1‑99 (MQ=255)
                                                aaaCTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTcgacga                                                                                                                                                >  3:104291/1‑104 (MQ=255)
                                                aaaCTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTcgacga                                                                                                                                                <  4:104291/104‑1 (MQ=255)
                                                aaaCTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTAc                                                                                                     <  3:35206/147‑1 (MQ=255)
                                                  aCTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACctgc                                                                                                 >  6:20120/1‑149 (MQ=255)
                                                      ctgctAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTg                                                                                                                                            >  6:10303/1‑102 (MQ=255)
                                                      ctgctAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTg                                                                                                                                            <  5:10303/102‑1 (MQ=255)
                                                      ctgctAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCtgctgc                                                                                              <  2:150552/148‑1 (MQ=255)
                                                      ctgctAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCtgctgc                                                                                              <  4:257265/148‑1 (MQ=255)
                                                         ctAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTAcc                                                                                                    >  8:103274/1‑139 (MQ=255)
                                                         ctAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTAcc                                                                                                    <  7:103274/139‑1 (MQ=255)
                                                         ctAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAAc                                                                                           <  3:121632/148‑1 (MQ=255)
                                                          tAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAAcgc                                                                                         >  8:87700/1‑149 (MQ=255)
                                                            aCCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTAc                                                                                                                                   >  4:347588/1‑105 (MQ=255)
                                                              cGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACgcggc                                                                                      >  8:117432/1‑148 (MQ=255)
                                                              cGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACgcggc                                                                                      >  7:309469/1‑148 (MQ=255)
                                                              cGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACgcggc                                                                                      <  7:117432/148‑1 (MQ=255)
                                                              cGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACgcggc                                                                                      >  8:60697/1‑148 (MQ=255)
                                                                 aGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGaaa                                                                                  <  7:87700/149‑1 (MQ=255)
                                                                   ccTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAAcc                                                                                >  3:5949/1‑149 (MQ=255)
                                                                    cTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACcg                                                                               >  3:72791/1‑149 (MQ=255)
                                                                       cGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACcgc                                                                              <  6:92547/147‑1 (MQ=255)
                                                                        gCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGTCGAGGTCGAGCTggt                                                                                                                                          >  5:171894/1‑85 (MQ=38)
                                                                        gCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGTCGAGGTCGAGCTggt                                                                                                                                          <  6:171894/86‑2 (MQ=38)
                                                                              gAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTg                                                                       >  2:113686/1‑147 (MQ=255)
                                                                              gAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTg                                                                       >  6:102107/1‑147 (MQ=255)
                                                                                  acCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTg                                                                 >  3:47982/1‑149 (MQ=255)
                                                                                   cccTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTgg                                                                <  7:60697/149‑1 (MQ=255)
                                                                                    ccTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTgg                                                                >  2:293553/1‑148 (MQ=255)
                                                                                       gCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGt                                                               >  4:230915/1‑146 (MQ=255)
                                                                                           tGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACg                                                         >  7:230171/1‑148 (MQ=255)
                                                                                              gcTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGAACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAg                                                      >  1:96856/1‑149 (MQ=255)
                                                                                              gcTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGc                                                     <  5:20120/149‑1 (MQ=255)
                                                                                                                 cAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAAc                                   >  2:271006/1‑148 (MQ=255)
                                                                                                                 cAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAAc                                   <  4:5949/148‑1 (MQ=255)
                                                                                                                   aGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAAc                                   >  7:93748/1‑146 (MQ=255)
                                                                                                                    gcaggcagTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGaa                                          >  2:84156/1‑138 (MQ=255)
                                                                                                                    gcaggcagTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGaa                                          <  1:84156/138‑1 (MQ=255)
                                                                                                                      aggcagTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACaagaag                             <  1:271006/149‑1 (MQ=255)
                                                                                                                      aggcagTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACaagaag                             <  8:309469/149‑1 (MQ=255)
                                                                                                                      aggcagTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGAAGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCTGCGAAACCGCGCGCCTGGAGCTCGTAAGACGCACCCTCTACCTGAAGGTGAACaag                                >  3:183031/1‑146 (MQ=255)
                                                                                                                          agTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCa                         <  4:72791/149‑1 (MQ=255)
                                                                                                                                 gCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGAc                    >  6:111299/1‑147 (MQ=255)
                                                                                                                                  ccAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCa                  <  3:230915/148‑1 (MQ=255)
                                                                                                                                         gCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCa          >  6:264578/1‑149 (MQ=255)
                                                                                                                                          ctcgacctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAg         <  8:230171/149‑1 (MQ=255)
                                                                                                                                           tcgacctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGa                     >  2:251531/1‑136 (MQ=255)
                                                                                                                                           tcgacctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGa                     <  1:251531/136‑1 (MQ=255)
                                                                                                                                            cgacctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGc        >  6:12007/1‑148 (MQ=255)
                                                                                                                                               cctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACctgc                                                                                                 <  6:259302/56‑1 (MQ=255)
                                                                                                                                               cctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACctgc                                                                                                 >  5:259302/1‑56 (MQ=255)
                                                                                                                                               cctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGCCAtg    >  4:124695/1‑149 (MQ=255)
                                                                                                                                               cctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGCCAtg    >  8:7040/1‑149 (MQ=255)
                                                                                                                                                ctcgacGAGCTGGACCCCTACGTGATGTTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCa                         >  2:158917/1‑127 (MQ=255)
                                                                                                                                                ctcgacGAGCTGGACCCCTACGTGATGTTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCa                         <  1:158917/127‑1 (MQ=255)
                                                                                                                                                ctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGcc       >  3:63167/1‑145 (MQ=255)
                                                                                                                                                ctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGcc       <  4:63167/145‑1 (MQ=255)
                                                                                                                                                 tcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGCCAtggc  >  1:269162/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CTGTGGCAATTGTTCAAAGGCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGCCATGGC  >  NC_002947/5957383‑5957676

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: