Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 635,170 C→T 6.0% L492L (CTC→CTT PP_0546 → sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947635,1700CT 6.0% 198.1 / 2.3 66L492L (CTC→CTTPP_0546sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base T (2/2);  total (32/34)
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

TGCGGGTGCTGGAGGAGCGCCAGGTGGTGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCGGGCGCGGCAGGCGCTGCTCGCATTTAACT  >  NC_002947/635028‑635312
                                                                                                                                               |                                                                                                                                              
tgcGGGTGCTGGAGGAGCGCCAGGTGGGGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGCCCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTTAATGGc                                                                                                                                          >  6:205944/1‑149 (MQ=255)
tgcGGGTGCTGGAGGAGCGCCAGGTGGCGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTTAATGGc                                                                                                                                          >  6:201384/1‑149 (MQ=255)
 gcGGGTGCTGGAGGAGCGCCAGGTGGTGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCGTCAATGGcc                                                                                                                                         >  8:93291/1‑149 (MQ=255)
   gggTGCTGGAGGAGCGCCAGGTGGTGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTg                                                                                                                                       <  4:90502/149‑1 (MQ=255)
       gCTGGAGGAGCGCCAGGTGGTGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAgg                                                                                                                                   >  5:17675/1‑149 (MQ=255)
        cTGGAGGAGCGCCAGGTGGTGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGt                                                                                                                                  <  8:145711/149‑1 (MQ=255)
        cTGGAGGAGCGCCAGGTGGTGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGt                                                                                                                                  <  4:10470/149‑1 (MQ=255)
              gagCGCCAGGTGGTGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCgcgc                                                                                                                             <  2:163145/148‑1 (MQ=255)
                    cAGGTGGTGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCcggc                                                                                                                      >  2:301001/1‑149 (MQ=255)
                    cAGGTGGTGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCcggc                                                                                                                      <  2:3604/149‑1 (MQ=255)
                          gtgCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTAcg                                                                                                                >  1:278316/1‑149 (MQ=255)
                            gCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGGTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACAGCCTCAATGGCCTGGAGGTAGCGCTGCCGGCGTTAc                                                                                                                 >  4:219733/1‑146 (MQ=255)
                               gATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGCGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGt                                                                                                                                  <  1:72022/127‑1 (MQ=255)
                               gATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGt                                                                                                                                  >  2:72022/1‑126 (MQ=255)
                               gATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGaa                                                                                                            >  6:62520/1‑148 (MQ=255)
                               gATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGaa                                                                                                            <  5:62520/148‑1 (MQ=255)
                               gATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAAc                                                                                                           >  1:218019/1‑149 (MQ=255)
                               gATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAAc                                                                                                           >  2:264315/1‑149 (MQ=255)
                                    gggtatagCCGCAGTCGGTTGACTTTCTCATCGTCAGCTCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGATTACCTGTATTACCGCCTCAATTGCATGGAGGTCGCGCTGCCGGCGGTACGCGAAc                                                                                                           <  3:223866/138‑1 (MQ=255)
                                    ggggAGATCCGCAGGCGGTGGACGTGAGCATAGTCAGCGCTACGCACCGTGACCTGCTGGAAAGGGTTGAACATGGCAGTTTAAGCTATGACATGTATTACAGCCTCAATGGCCTGGAGTTCGCTCTGCCGGAGTTACGCGaca                                                                                                           >  4:223866/1‑142 (MQ=255)
                                    ggggAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGt                                                                                                      <  4:219416/149‑1 (MQ=255)
                                      ggAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGa                                                                                                    >  6:103512/1‑149 (MQ=255)
                                       gagaGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGa                                                                                                    >  8:26060/1‑148 (MQ=255)
                                           gCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAgg                                                                                               <  2:278316/149‑1 (MQ=255)
                                             cGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGgcg                                                                                             >  5:117852/1‑149 (MQ=255)
                                              gCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGgcg                                                                                             >  2:193718/1‑148 (MQ=255)
                                              gCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGgcg                                                                                             >  7:259150/1‑148 (MQ=255)
                                                    ggTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCcggcg                                                                                                                     <  1:199237/118‑1 (MQ=255)
                                                    ggTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCcggcg                                                                                                                     >  2:199237/1‑118 (MQ=255)
                                                     gTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCt                                                                                     <  6:117852/149‑1 (MQ=255)
                                                     gTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCt                                                                                     <  7:26060/149‑1 (MQ=255)
                                                        gACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGa                                                                                  <  4:167781/149‑1 (MQ=255)
                                                             gcgcATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTc                                                                             >  6:48912/1‑149 (MQ=255)
                                                              cgcATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTc                                                                             <  1:301001/148‑1 (MQ=255)
                                                               gcATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTTAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGTAGTGACAAGGCGCAGTTGCTGGACTTTctg                                                                           <  5:201384/149‑1 (MQ=255)
                                                               gcATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTTAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGTAGTGACAAGGCGCAGTTGCTGGACTTTctg                                                                           <  5:205944/149‑1 (MQ=255)
                                                                 aTCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTctgct                                                                         >  3:132747/1‑149 (MQ=255)
                                                                     tCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTctgct                                                                         <  2:122605/145‑1 (MQ=255)
                                                                     tCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTctgct                                                                         >  1:122605/1‑145 (MQ=255)
                                                                       aGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCa                                                                   <  1:221527/149‑1 (MQ=255)
                                                                        gcgcTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAg                                                                  <  8:218252/149‑1 (MQ=255)
                                                                                        ccTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAAt                                                  <  2:218019/149‑1 (MQ=255)
                                                                                         ctgctgGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATg                                                 <  1:99643/149‑1 (MQ=255)
                                                                                         ctgctgGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATg                                                 >  2:99643/1‑149 (MQ=255)
                                                                                          tgctgGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATg                                                 >  5:327281/1‑148 (MQ=255)
                                                                                            ctgGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGa                                               <  4:174559/148‑1 (MQ=255)
                                                                                             tgGAACGGGTTGAACAGGGCAGTGTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATa                                             <  8:259150/149‑1 (MQ=255)
                                                                                               gAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGa                                           >  6:153568/1‑149 (MQ=255)
                                                                                               gAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGa                                           <  6:17675/149‑1 (MQ=255)
                                                                                               gAACGGGTTGAACAGGGCAGTTTCCGAGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAg                                            <  4:164695/148‑1 (MQ=255)
                                                                                                     gTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCTGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGa                                     <  3:219733/149‑1 (MQ=255)
                                                                                                             gggCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCggg                             >  7:281315/1‑149 (MQ=255)
                                                                                                                    tttCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCGGGCGCgg                        >  6:198504/1‑148 (MQ=255)
                                                                                                                    tttCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCACGCGCAATCGATAGACATCGAGCCCCGGGCGCgg                        >  5:300808/1‑148 (MQ=255)
                                                                                                                     ttCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACttt                                                                              >  3:170125/1‑93 (MQ=255)
                                                                                                                     ttCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACttt                                                                              <  4:170125/93‑1 (MQ=255)
                                                                                                                         gcgAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCGGGCGCGGCAGgcg                  <  1:193718/149‑1 (MQ=255)
                                                                                                                         gcgAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCGGGCGCGGCAGgcg                  <  5:103512/149‑1 (MQ=255)
                                                                                                                                   gTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCGGGCGCGGCAGGCGCTGCTCGCAt        >  3:254406/1‑149 (MQ=255)
                                                                                                                                   gTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCGGGCGCGGCAGGCGCTGCTCGCAt        >  3:104709/1‑149 (MQ=255)
                                                                                                                                    tATTACCGCCTCACTGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCGGGCGCGGCAGGCGCTGCTCGCAt        <  6:40568/148‑1 (MQ=255)
                                                                                                                                    tATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACa                                         <  7:43243/115‑1 (MQ=255)
                                                                                                                                    tATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACa                                         >  8:43243/1‑115 (MQ=255)
                                                                                                                                        aCCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCGGGCGCGGCAGGCGCTGCTCGCATTTaa    <  4:104709/148‑1 (MQ=255)
                                                                                                                                        aCCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCGGGCGCGGCAGGCGCTGCTCGCATTTaa    <  5:48912/148‑1 (MQ=255)
                                                                                                                                         ccGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCGGGCGCGGCAGGCGCTGCTCGCATTTAACt  >  2:133365/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TGCGGGTGCTGGAGGAGCGCCAGGTGGTGCCGATCGGGGGAGAGCCGCAGGCGGTGGACGTGCGCATCGTCAGCGCTACGCACCGTGACCTGCTGGAACGGGTTGAACAGGGCAGTTTCCGCGAGGACCTGTATTACCGCCTCAATGGCCTGGAGGTCGCGCTGCCGGCGGTACGCGAACGCAGTGACAAGGCGCAGTTGCTGGACTTTCTGCTGCGTCAGGAGACGCAGGGGCAATGGATAGACATCGAGCCCCGGGCGCGGCAGGCGCTGCTCGCATTTAACT  >  NC_002947/635028‑635312

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: