Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,957,531 A→T 8.3% D312V (GAC→GTC)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,957,5310AT 8.3% 169.3 / 5.9 60D312V (GAC→GTC) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (30/25);  new base T (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

CTGTGGCAATTGTTCAAAGGCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGCCATGGCTGG  >  NC_002947/5957383‑5957679
                                                                                                                                                     |                                                                                                                                                    
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)
                                                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                                                                                                 <  4:257265/148‑1 (MQ=255)
                                                      ctgctAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCtgctgc                                                                                                 <  2:150552/148‑1 (MQ=255)
                                                         ctAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTAcc                                                                                                       <  7:103274/139‑1 (MQ=255)
                                                         ctAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTAcc                                                                                                       >  8:103274/1‑139 (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:60697/1‑148 (MQ=255)
                                                              cGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACgcggc                                                                                         >  8:117432/1‑148 (MQ=255)
                                                              cGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACgcggc                                                                                         >  7:309469/1‑148 (MQ=255)
                                                              cGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACgcggc                                                                                         <  7:117432/148‑1 (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                        <  1:251531/136‑1 (MQ=255)
                                                                                                                                           tcgacctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGa                        >  2:251531/1‑136 (MQ=255)
                                                                                                                                            cgacctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGc           >  6:12007/1‑148 (MQ=255)
                                                                                                                                               cctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACctgc                                                                                                    >  5:259302/1‑56 (MQ=255)
                                                                                                                                               cctcgacGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACctgc                                                                                                    <  6:259302/56‑1 (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)
                                                                                                                                                     acgaGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGCCAtggctgg  <  5:264578/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGTGGCAATTGTTCAAAGGCATCGACTCGCCCTACAAATCGCTGCTCAAACTGCTGCTAACCGAAGCCTACGCCAGCGAACACCCTGCGGTGCGCTGCCTGAGCCTGGACTACAAGCAGGCAGTGTTTGCCAACCAGCTCGACCTCGACGAGCTGGACCCCTACGTGATGGTCTACCGGCGCATCGAACGTTACCTGCTGCAACGCGGCGAAACCGCGCGCCTGGAGCTGGTAAGACGCAGCCTGTACCTGAAGGTGAACAAGAAGCTCAGCGACCAGGGCCGCAGCCATGGCTGG  >  NC_002947/5957383‑5957679

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: