Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 971,678 A→T 6.8% intergenic (+33/‑31) 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,6780AT 6.8% 213.8 / 5.1 73intergenic (+33/‑31)yajC/secDprotein translocase subunit YajC/protein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base A (32/36);  new base T (2/3);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

AAGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGAT  >  NC_002947/971535‑971818
                                                                                                                                               |                                                                                                                                            
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:57245/149‑1 (MQ=255)
                              gtgtCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCt                                                                                                           <  1:6693/149‑1 (MQ=255)
                                 tCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGaa                                                                                                        >  2:69286/1‑149 (MQ=255)
                                 tCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGaa                                                                                                        >  5:222128/1‑149 (MQ=255)
                                 tCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGaa                                                                                                        <  2:258294/149‑1 (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)
                                                                                                                                               |                                                                                                                                            
AAGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGAT  >  NC_002947/971535‑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: