Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 971,677 C→T 7.1% intergenic (+32/‑32) yajC → / → secD protein translocase subunit YajC/protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947971,6770CT 7.1% 198.2 / 4.4 71intergenic (+32/‑32)yajC/secDprotein translocase subunit YajC/protein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base T (2/3);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

TCGTCAAGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGAT  >  NC_002947/971530‑971818
                                                                                                                                                   |                                                                                                                                             
tCGTCAAGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGgcgcgc                                                                                                                                               <  3:345115/148‑1 (MQ=255)
     aaGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGTTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGTTTTTATTTTCCAGTCGGGGCGCGCAAcgcg                                                                                                                                         <  5:408624/149‑1 (MQ=255)
      aGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACgcgc                                                                                                                                        <  8:78873/149‑1 (MQ=255)
         tttCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGTTGCGCGcccc                                                                                                                                     <  5:216500/149‑1 (MQ=255)
         tttCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGTTGCGCGcccc                                                                                                                                     <  6:402317/149‑1 (MQ=255)
         tttCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGcccc                                                                                                                                     >  1:258294/1‑149 (MQ=255)
             cGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGt                                                                                                                                 >  8:285321/1‑149 (MQ=255)
              gATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTc                                                                                                                                >  5:235363/1‑149 (MQ=255)
              gATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGAACGCTCAAGGCTATCTGAGTTACCGGTTTAATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTc                                                                                                                                <  1:350291/149‑1 (MQ=255)
               aTGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCTAGGCTATCTGAGTTACCGGTTTTATTTTCCTGTCGGGGCGCGCAACGCGCCCCGCGTCt                                                                                                                               <  1:225400/149‑1 (MQ=255)
                tGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCt                                                                                                                               <  8:19975/148‑1 (MQ=255)
                               ggAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGa                                                                                                               >  2:59661/1‑149 (MQ=255)
                                gAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGAt                                                                                                              <  2:120580/149‑1 (MQ=255)
                                gAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGAt                                                                                                              >  4:202057/1‑149 (MQ=255)
                                gAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGAt                                                                                                              >  3:230637/1‑149 (MQ=255)
                                gAAGGGTCCGACACTGTAGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCAAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGAt                                                                                                              <  4:233115/149‑1 (MQ=255)
                                 aaGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATg                                                                                                             <  1:59661/149‑1 (MQ=255)
                                   gtgtCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCt                                                                                                           <  1:6693/149‑1 (MQ=255)
                                   gtgtCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCt                                                                                                           <  1:57245/149‑1 (MQ=255)
                                      tCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGaa                                                                                                        <  2:258294/149‑1 (MQ=255)
                                      tCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGaa                                                                                                        >  2:69286/1‑149 (MQ=255)
                                      tCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGaa                                                                                                        >  5:222128/1‑149 (MQ=255)
                                        cGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAAc                                                                                                       >  8:279056/1‑148 (MQ=255)
                                             cTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAatat                                                                                                 <  1:69286/149‑1 (MQ=255)
                                              tGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATc                                                                                                >  1:186830/1‑149 (MQ=255)
                                                  gAGCTTAAGTGCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAACGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCCCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCctct                                                                                            <  4:230637/149‑1 (MQ=255)
                                                   aGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCtg                                                                                           >  1:44340/1‑149 (MQ=255)
                                                     cTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCtgtg                                                                                         >  7:221140/1‑149 (MQ=255)
                                                      tGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCtgtg                                                                                         >  1:134208/1‑148 (MQ=255)
                                                        aaGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCtgtg                                                                                         <  7:260523/146‑1 (MQ=255)
                                                        aaGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCGGTGGaa                                                                                      >  8:355340/1‑149 (MQ=255)
                                                         aGTTCCATATGGGTGCCATTGCTGCTACCCTGCAAAAAGGTAGGCTCAAGGCTATCTCAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGaa                                                                                      <  6:26997/148‑1 (MQ=255)
                                                          gTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGaaa                                                                                     >  4:67294/1‑148 (MQ=255)
                                                             ccAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAatat                                                                                  <  3:67294/148‑1 (MQ=255)
                                                             ccAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAatat                                                                                  <  5:128507/148‑1 (MQ=255)
                                                                        ccATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTg                                                                      <  7:279056/149‑1 (MQ=255)
                                                                           ttGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTg                                                                   <  2:134208/149‑1 (MQ=255)
                                                                              ctgcGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTa                                                                >  6:194153/1‑149 (MQ=255)
                                                                                     ccTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCATGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACtggtgg                                                         >  7:264359/1‑149 (MQ=255)
                                                                                         ccAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCg                                                      <  5:106931/148‑1 (MQ=255)
                                                                                         caaaaaGGTACGCCCACGCCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCg                                                      <  8:139295/146‑1 (MQ=255)
                                                                                          aaaacGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCTGGCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGt                                                    <  2:186830/148‑1 (MQ=255)
                                                                                            aaaGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGtt                                                   >  2:291716/1‑148 (MQ=255)
                                                                                              aGGTACGCTCAAGGCTATGTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACtggtg                                                          <  7:201184/139‑1 (MQ=255)
                                                                                              aGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACtggtg                                                          >  8:201184/1‑139 (MQ=255)
                                                                                              aGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCa                                                <  4:67257/149‑1 (MQ=255)
                                                                                               ggTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCAt                                               <  8:221140/149‑1 (MQ=255)
                                                                                                   cgcgCATGGGTATGTGAGTTACCGGTTTTTTTTACCAGTCGGGGCGCGAAACGAGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATAAGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTAt                                           <  6:235363/149‑1 (MQ=255)
                                                                                                   cGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCATCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTa                                            <  5:177946/148‑1 (MQ=255)
                                                                                                        aaGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTcc                                        <  5:194153/147‑1 (MQ=255)
                                                                                                            cTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTcc                                   >  5:58067/1‑148 (MQ=255)
                                                                                                               tCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAAc                               >  7:284291/1‑149 (MQ=255)
                                                                                                                  gAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAAc                               <  1:368922/146‑1 (MQ=255)
                                                                                                                  gAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAAc                               >  2:368922/1‑146 (MQ=255)
                                                                                                                   aGTTACCGGTTTTATTTTCCAGGCGTGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACATATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGTTGGTCGGTTTCATTTATTCCGCTCCCAAc                               <  8:264359/145‑1 (MQ=255)
                                                                                                                     ttACCGGTTTTATTTTCCAGTCGGGGCGTGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTAc                         >  5:327639/1‑149 (MQ=255)
                                                                                                                         cGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCg                      >  6:116678/1‑148 (MQ=255)
                                                                                                                          ggTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGa                    >  2:54616/1‑149 (MQ=255)
                                                                                                                          ggTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGa                    <  7:285321/149‑1 (MQ=255)
                                                                                                                            ttttATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCGTTTATTCCGCTCCCAACCTCTACCCGGAt                   >  1:38721/1‑148 (MQ=255)
                                                                                                                            ttttATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATg                  >  6:129639/1‑149 (MQ=255)
                                                                                                                            ttttATTTTCCAGTCGGGGAGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGGCGGTTTCATTTATTCCGCTCCCAACCTCTACACGGAt                   >  5:193071/1‑148 (MQ=255)
                                                                                                                                aTTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGAccc              >  1:16672/1‑149 (MQ=255)
                                                                                                                                      cAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCcggcgg         <  6:222128/148‑1 (MQ=255)
                                                                                                                                      cAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGAACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCcggcgg         <  2:44340/148‑1 (MQ=255)
                                                                                                                                       aGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTa       <  2:104867/149‑1 (MQ=255)
                                                                                                                                           ggggCGCGTTGCGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTa       <  3:23179/145‑1 (MQ=255)
                                                                                                                                           ggggCGCGTTGCGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTa       >  4:23179/1‑145 (MQ=255)
                                                                                                                                           ggggCGCGTTGCGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTa       >  4:236495/1‑145 (MQ=255)
                                                                                                                                            gggCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTa       >  2:46212/1‑144 (MQ=255)
                                                                                                                                            gggCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGAt  <  3:202057/149‑1 (MQ=255)
                                                                                                                                            gggCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACACGGATGACCCGGCGGTACAGa   >  4:360683/1‑148 (MQ=255)
                                                                                                                                                gcgcAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTAc      >  3:89850/1‑141 (MQ=255)
                                                                                                                                                gcgcAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTAc      <  4:89850/141‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TCGTCAAGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGAT  >  NC_002947/971530‑971818

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