Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 502,123 A→T 7.8% K138* (AAA→TAA)  PP_0413 → polyamine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947502,1230AT 7.8% 207.9 / 6.3 76K138* (AAA→TAA) PP_0413polyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (33/37);  new base T (3/3);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAGGCCAC  >  NC_002947/501981‑502271
                                                                                                                                               |                                                                                                                                                    
aaGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAagag                                                                                                                                                 <  3:274211/147‑1 (MQ=255)
 aGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGc                                                                                                                                               <  3:242410/149‑1 (MQ=255)
   ccAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAagag                                                                                                                                                 >  3:76306/1‑145 (MQ=255)
   ccAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAagag                                                                                                                                                 <  4:76306/145‑1 (MQ=255)
     aGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTg                                                                                                                                           >  3:140576/1‑149 (MQ=255)
     aGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTg                                                                                                                                           <  1:57159/149‑1 (MQ=255)
       gATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCa                                                                                                                                         <  4:94519/149‑1 (MQ=255)
            ggCAGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGAATGCCCGTGCGCTGCACTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCCTTGCGGGCTg                                                                                                                                    >  2:174900/1‑149 (MQ=255)
             gCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTg                                                                                                                                    <  3:27014/148‑1 (MQ=255)
               cGAGTCGCGCAAGAACCAGCCCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTTTAAGAGGCGGTGCAGGCTGTTg                                                                                                                                 >  4:210001/1‑149 (MQ=255)
                   tCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGa                                                                                                                             <  5:337853/149‑1 (MQ=255)
                       gcAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGt                                                                                                                         >  4:364572/1‑149 (MQ=255)
                        cAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGtt                                                                                                                        >  3:295580/1‑149 (MQ=255)
                            aaCCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTgggg                                                                                                                    <  3:210001/149‑1 (MQ=255)
                                           gACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTg                                                                                                     >  4:225163/1‑149 (MQ=255)
                                              cTTTCCAAACGCTTGAACATGGAGCTGGCTCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTAtt                                                                                                      >  1:221456/1‑146 (MQ=255)
                                              cTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGTCTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTAtt                                                                                                      <  2:221456/145‑1 (MQ=255)
                                              cTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCa                                                                                                  <  4:140576/149‑1 (MQ=255)
                                                 tCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAg                                                                                                 <  4:314030/147‑1 (MQ=255)
                                                    aaaCGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGc                                                                                                                           <  8:213438/118‑1 (MQ=255)
                                                                aTGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGa                                                                                                                                <  2:303983/101‑1 (MQ=255)
                                                                aTGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGa                                                                                                                                >  1:303983/1‑101 (MQ=255)
                                                                aTGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCg                                                                                              <  3:133610/135‑1 (MQ=255)
                                                                aTGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCg                                                                                              >  4:133610/1‑135 (MQ=255)
                                                                aTGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATAc                                                                                >  1:205670/1‑149 (MQ=255)
                                                                  ggAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCaa                                                                                   >  8:221319/1‑144 (MQ=255)
                                                                  ggAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCaa                                                                                   <  7:221319/144‑1 (MQ=255)
                                                                     gCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGtt                                                                           >  1:288580/1‑149 (MQ=255)
                                                                      cTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTc                                                                          <  4:295580/149‑1 (MQ=255)
                                                                       tGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTg                                                                                                                       <  5:338187/103‑1 (MQ=255)
                                                                       tGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTg                                                                                                                       >  6:338187/1‑103 (MQ=255)
                                                                        ggCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGCCGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCgg                                                                        >  5:64885/1‑149 (MQ=255)
                                                                         gCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCg                                                                         <  7:218216/147‑1 (MQ=255)
                                                                             ggTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGAc                                                                    >  5:222172/1‑149 (MQ=255)
                                                                                 aCCGCAGCGTGTTGACCTAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGActtc                                                                 <  2:323404/148‑1 (MQ=255)
                                                                                  ccGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTAtt                                                                                                      <  5:323697/110‑1 (MQ=255)
                                                                                  ccGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTAtt                                                                                                      >  6:323697/1‑110 (MQ=255)
                                                                                  ccGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGActtctt                                                               <  2:205670/149‑1 (MQ=255)
                                                                                   cGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                     <  5:28285/94‑1 (MQ=255)
                                                                                   cGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                     >  6:28285/1‑94 (MQ=255)
                                                                                    gCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACttctttc                                                             <  3:364572/149‑1 (MQ=255)
                                                                                     cAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACttctttct                                                            <  1:192327/149‑1 (MQ=255)
                                                                                     cAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACttctttc                                                             >  1:64168/1‑148 (MQ=255)
                                                                                        cTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATAc                                                                                <  4:190229/126‑1 (MQ=255)
                                                                                        cTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATAc                                                                                >  3:190229/1‑126 (MQ=255)
                                                                                           gttgACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCCGGCGATACGCCGCAATACCAGTTCGGTGCCGTCTTTCTACctcc                                                      >  4:277552/1‑147 (MQ=255)
                                                                                            ttgACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACctgct                                                     >  3:230159/1‑149 (MQ=255)
                                                                                              gACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTgg                                                   >  5:58644/1‑149 (MQ=255)
                                                                                               aCCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGActtct                                                                >  4:34767/1‑135 (MQ=255)
                                                                                               aCCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGActtct                                                                <  3:34767/135‑1 (MQ=255)
                                                                                                    gACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCgg                                             >  6:307563/1‑149 (MQ=255)
                                                                                                     aCTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGt                                            <  4:230159/149‑1 (MQ=255)
                                                                                                                  tGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTAcc                                                         >  4:313994/1‑123 (MQ=255)
                                                                                                                  tGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTAcc                                                         <  3:313994/123‑1 (MQ=255)
                                                                                                                  tGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatc                                    <  7:138305/144‑1 (MQ=255)
                                                                                                                  tGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatc                                    >  8:138305/1‑144 (MQ=255)
                                                                                                                      cTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGAt                            >  2:286158/1‑148 (MQ=255)
                                                                                                                         cAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACctgctg                                                    <  8:1309/121‑1 (MQ=255)
                                                                                                                         cAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACctgctg                                                    >  7:1309/1‑121 (MQ=255)
                                                                                                                         cAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGc                                                  <  5:56632/123‑1 (MQ=255)
                                                                                                                         cAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGc                                                  >  6:56632/1‑123 (MQ=255)
                                                                                                                         cAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGAcc                        >  1:178652/1‑149 (MQ=255)
                                                                                                                            ggCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTg                     <  6:64885/149‑1 (MQ=255)
                                                                                                                              cTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGTCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGt                   <  2:288580/149‑1 (MQ=255)
                                                                                                                              cTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGt                   <  6:222172/149‑1 (MQ=255)
                                                                                                                                 gcaacgcctCTTTATAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGTACcatc                                    <  3:188082/124‑1 (MQ=255)
                                                                                                                                 gcaacgcctCTTTATAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatc                                    >  4:188082/6‑129 (MQ=255)
                                                                                                                                   cctgcctCTTATAAAGAGGCGTTGCAGGCTGTTGAGGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGc               <  5:180747/148‑1 (MQ=255)
                                                                                                                                     cgcctCTTTATAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatcgcat                                <  3:151836/128‑1 (MQ=255)
                                                                                                                                     cgcctCTTTATAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatcgcat                                >  4:151836/2‑129 (MQ=255)
                                                                                                                                        ctctTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGa                             >  7:44947/1‑129 (MQ=255)
                                                                                                                                        ctctTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGa                             <  8:44947/129‑1 (MQ=255)
                                                                                                                                         tctTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGATAGcc          >  3:328070/1‑147 (MQ=255)
                                                                                                                                         tctTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGAAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACGTTGGTGAGCTAGCCaa        <  6:58644/149‑1 (MQ=255)
                                                                                                                                         tctTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTAGGGGCGACCCGGCTTATTGGCACGCGATACGCCGCAATACCAGTTCGGTGAGTTCTTTCTACCTGCTGGACTCGGTGGACCCTCGCATCGATGACCTAGGTGAGCTAGcc          >  2:120109/1‑147 (MQ=255)
                                                                                                                                           ttATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCaa        <  3:225163/147‑1 (MQ=255)
                                                                                                                                            tataAAGAGGCGTTGCAGGCTGTTGACGCGCGTTTGGGCGACCAGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTACTTTCTACCTGATGGCCTCGGTGGCACATCGCATAGATGACATTGGTGAGCTAGCCACgg      >  6:187981/1‑148 (MQ=255)
                                                                                                                                               aaaGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAGGCCAc  >  4:333248/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAGGCCAC  >  NC_002947/501981‑502271

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: