Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 971,679 A→G 6.8% intergenic (+34/‑30) 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,6790AG 6.8% 210.4 / 4.2 73intergenic (+34/‑30)yajC/secDprotein translocase subunit YajC/protein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base A (32/36);  new base G (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.70e-01

AAGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGATCAGCGGTGC  >  NC_002947/971535‑971827
                                                                                                                                                |                                                                                                                                                    
aaGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGTTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGTTTTTATTTTCCAGTCGGGGCGCGCAAcgcg                                                                                                                                                  <  5:408624/149‑1 (MQ=255)
 aGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACgcgc                                                                                                                                                 <  8:78873/149‑1 (MQ=255)
    tttCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGTTGCGCGcccc                                                                                                                                              <  6:402317/149‑1 (MQ=255)
    tttCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGTTGCGCGcccc                                                                                                                                              <  5:216500/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                                                                                                                       >  4:202057/1‑149 (MQ=255)
                           gAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGAt                                                                                                                       <  2:120580/149‑1 (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                                                                                                                 >  5:222128/1‑149 (MQ=255)
                                 tCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGaa                                                                                                                 >  2:69286/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                                                                                           <  5:128507/148‑1 (MQ=255)
                                                        ccAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAatat                                                                                           <  3:67294/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                                        >  2:368922/1‑146 (MQ=255)
                                                                                                             gAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAAc                                        <  1:368922/146‑1 (MQ=255)
                                                                                                              aGTTACCGGTTTTATTTTCCAGGCGTGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACATATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGTTGGTCGGTTTCATTTATTCCGCTCCCAAc                                        <  8:264359/145‑1 (MQ=255)
                                                                                                                ttACCGGTTTTATTTTCCAGTCGGGGCGTGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTAc                                  >  5:327639/1‑149 (MQ=255)
                                                                                                                    cGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCg                               >  6:116678/1‑148 (MQ=255)
                                                                                                                     ggTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGa                             <  7:285321/149‑1 (MQ=255)
                                                                                                                     ggTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGa                             >  2:54616/1‑149 (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                >  4:23179/1‑145 (MQ=255)
                                                                                                                                      ggggCGCGTTGCGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTa                <  3:23179/145‑1 (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               <  4:89850/141‑1 (MQ=255)
                                                                                                                                           gcgcAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTAc               >  3:89850/1‑141 (MQ=255)
                                                                                                                                                aCGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGATCAGCGGTGc  <  8:284291/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AAGGTTTCCGATGACTTCGTGGTGCTGGAAGTGTCCGACACTGTCGAGCTGAAGTTCCAGAAGGGTGCCATTGCTGCGACCCTGCCAAAAGGTACGCTCAAGGCTATCTGAGTTACCGGTTTTATTTTCCAGTCGGGGCGCGCAACGCGCCCCGCGTCTTGAACGGGCGGCGTGATGCTGAACAAATATCCTCTGTGGAAATATGCACTGATCGTGCTGGTACTGGTGGTCGGTTTCATTTATTCCGCTCCCAACCTCTACCCGGATGACCCGGCGGTACAGATCAGCGGTGC  >  NC_002947/971535‑971827

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: