Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,935,077 T→A 6.1% E112D (GAA→GAT PP_3467 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,935,0770TA 6.1% 196.4 / 2.7 66E112D (GAA→GATPP_3467sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (28/34);  new base A (2/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

CGAGGTAATGTTCGCCGCAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCCAGCAGCAGGCGCAGGTCATGGCCCTGTT  >  NC_002947/3934929‑3935225
                                                                                                                                                      |                                                                                                                                                     
cGAGGTAATGTTCGCCGCAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGACGCGCAGCGATTAGCTACAGCAACTGCCCGGGTGCAGATCGCGAAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCt                                                                                                                                                       >  1:158683/1‑149 (MQ=255)
   ggTAATGTTCGCCGCAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTc                                                                                                                                                     >  6:250785/1‑148 (MQ=255)
    gTAATGTTCGCCGCAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGa                                                                                                                                                   <  8:36219/149‑1 (MQ=255)
    gTAATGTTCGCCGCAGATGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGa                                                                                                                                                   >  1:112767/1‑149 (MQ=255)
         gTTCGCCGCAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGc                                                                                                                                               >  1:95168/1‑149 (MQ=255)
          ttCGCCGCAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGc                                                                                                                                               >  8:313582/1‑147 (MQ=255)
          ttCGCCGCAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCa                                                                                                                                             >  4:33380/1‑149 (MQ=255)
          ttCGCCGCAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCa                                                                                                                                             <  7:31624/149‑1 (MQ=255)
             gccgcAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCa                                                                                                                                          >  1:326815/1‑149 (MQ=255)
               cgcAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCATCGAAGGCCAgaagc                                                                                                                                        <  5:77379/149‑5 (MQ=255)
               cgcAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTGCTGGCCATCGAAGGCCAgaagc                                                                                                                                        <  6:138123/149‑5 (MQ=255)
                   gTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTgg                                                                                                                                    >  3:63079/1‑149 (MQ=255)
                    ttGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTGCTGGCCTTCGATGGCCAGCACGTTGGc                                                                                                                                   <  5:310899/149‑1 (MQ=255)
                        tCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTgg                                                                                                                               <  6:339983/149‑1 (MQ=255)
                        tCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTgg                                                                                                                               <  4:168127/149‑1 (MQ=255)
                                ggccggccTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGcc                                                                                                                       >  7:93926/1‑149 (MQ=255)
                                    ggccTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAg                                                                                                                    >  2:112447/1‑148 (MQ=255)
                                     gccTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGa                                                                                                                  <  2:326815/149‑1 (MQ=255)
                                      ccTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGAt                                                                                                                 >  7:343066/1‑149 (MQ=255)
                                        ttCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACg                                                                                                                                        <  6:57911/124‑1 (MQ=255)
                                        ttCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACg                                                                                                                                        >  5:57911/1‑124 (MQ=255)
                                            aCCACGGCCGTGCCGATGGCATTGGTGCCGCGGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGa                                                                                                            >  1:253692/1‑149 (MQ=255)
                                              cACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGc                                                                                                         <  5:250785/149‑1 (MQ=255)
                                                        ccGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGcc                                                                                               >  7:264412/1‑149 (MQ=255)
                                                        ccGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGcc                                                                                               <  5:281475/149‑1 (MQ=255)
                                                         cGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGccc                                                                                              >  2:95105/1‑149 (MQ=255)
                                                             tGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCGTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTgg                                                                                           <  4:63079/149‑1 (MQ=255)
                                                               gCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGgc                                                                                          <  8:208370/148‑1 (MQ=255)
                                                                 aTTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGgcgcg                                                                                       <  8:96453/149‑1 (MQ=255)
                                                                 aTTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGgcgcg                                                                                       >  3:1371/1‑149 (MQ=255)
                                                                 aTTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGgcgcg                                                                                       >  1:244811/1‑149 (MQ=255)
                                                                     gTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGAc                                                                                   <  3:33380/149‑1 (MQ=255)
                                                                         cgcgCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAgg                                                                               >  8:82602/1‑149 (MQ=255)
                                                                           cgcAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAgg                                                                               >  6:4395/1‑147 (MQ=255)
                                                                             gcAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAg                                                                            <  8:93926/149‑1 (MQ=255)
                                                                                 cGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTc                                                                         >  1:26872/1‑147 (MQ=255)
                                                                                  gATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAg                                                                      <  5:4395/149‑1 (MQ=255)
                                                                                       gCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCgg                                                                 <  1:112447/149‑1 (MQ=255)
                                                                                         tCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCgggg                                                               >  8:132170/1‑149 (MQ=255)
                                                                                            aGCAACTGCCCGGGTGCAGATCGCGCAGGCCTACGCTTTTCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGt                                                             >  6:247014/1‑148 (MQ=255)
                                                                                               aaCTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACg                                                         <  7:82602/149‑1 (MQ=255)
                                                                                                    cccGGGTGCAGCTCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACCTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATc                                                    >  6:103037/1‑149 (MQ=255)
                                                                                                     ccGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCg                                                   <  2:26872/149‑1 (MQ=255)
                                                                                                       gggTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGGATCGAt                                                 <  4:218044/149‑1 (MQ=255)
                                                                                                        ggTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATc                                                <  2:95168/149‑1 (MQ=255)
                                                                                                         gTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCgggg                                                               >  3:167620/1‑133 (MQ=255)
                                                                                                         gTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCgggg                                                               <  4:167620/133‑1 (MQ=255)
                                                                                                         gTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATca                                               <  3:192292/149‑1 (MQ=255)
                                                                                                            cAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATcagc                                             <  4:293410/148‑1 (MQ=255)
                                                                                                                  gcgcAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCgg                                       <  3:108088/148‑1 (MQ=255)
                                                                                                                          cTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTc                               >  8:347385/1‑148 (MQ=255)
                                                                                                                             cgcgTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCca                             <  8:264412/147‑1 (MQ=255)
                                                                                                                              gcgTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGCTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCcagc                           <  1:95105/149‑1 (MQ=255)
                                                                                                                                       caCGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCCAGCAGCAGGCGCa                 <  7:132170/149‑1 (MQ=255)
                                                                                                                                        aCGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCCAGCAGCAGGCGCAg                <  2:112767/149‑1 (MQ=255)
                                                                                                                                         cGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCTGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCCAGCAGCAGGCGCAgg               <  4:1371/149‑1 (MQ=255)
                                                                                                                                         cGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCCAGCAGCAGGCGCAgg               <  2:244811/149‑1 (MQ=255)
                                                                                                                                              tgcTGGCCTTCGATGGCCAGCAGGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGtta                                    <  8:103090/121‑2 (MQ=255)
                                                                                                                                              tgcTGGCCATCGAAGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCCAGCAGCAGGCGCAGGTCATg          >  3:253682/3‑149 (MQ=255)
                                                                                                                                               gcTGGCCATCGAAGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCCAGCAGCAGGCGCAGGTCATgg         >  3:205920/2‑149 (MQ=255)
                                                                                                                                                  ggCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTcc                              >  2:244574/1‑125 (MQ=255)
                                                                                                                                                  ggCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTcc                              <  1:244574/125‑1 (MQ=255)
                                                                                                                                                   gCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTgcgc                                 <  7:266336/121‑1 (MQ=255)
                                                                                                                                                   gCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTgcgc                                 >  8:266336/1‑121 (MQ=255)
                                                                                                                                                   gCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCCAGCAGCAGGCGCAGGTCATGGCCCt     <  5:103037/149‑1 (MQ=255)
                                                                                                                                                     cTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCcagcagca                       >  2:217978/1‑129 (MQ=255)
                                                                                                                                                     cTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCCAGCAGCAGGCGCAGGTCATGGCCCTGt   >  7:117023/1‑149 (MQ=255)
                                                                                                                                                      ttCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAg                                                                                >  4:96065/1‑72 (MQ=255)
                                                                                                                                                      ttCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAg                                                                                <  3:96065/72‑1 (MQ=255)
                                                                                                                                                      ttCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCCAGCAGCAGGCGCAGGTCATGGCCCTGtt  >  3:368713/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CGAGGTAATGTTCGCCGCAGTTGATCAGCGTCGGCCGGCCTTCGACCACGGCCGTGCCGATGGCATTGGTGCCGCGCAGCGATTCGCTCCAGCAACTGCCCGGGTGCAGATCGCGCAGGCCTTCGCGTTGCAGCACGTGCTTCTGGCCTTCGATGGCCAGCACGTTGGCCTGGGCATCGCCGAGGATGACGATGCCAGCCTTGCCCTGGCGCGCGACCAGGTAGTCGAGTTCGGGGGTAACGGCATCGATCAGCAGGCGGTTGCGCTCCAGCAGCAGGCGCAGGTCATGGCCCTGTT  >  NC_002947/3934929‑3935225

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: