Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,173,233 A→G 6.1% intergenic (‑152/+38) PP_5410 ← / ← glmU DeoR family transcriptional regulator/bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,173,2330AG 6.1% 193.7 / 2.5 65intergenic (‑152/+38)PP_5410/glmUDeoR family transcriptional regulator/bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base G (2/2);  total (33/32)
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.79e-01

AAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGC  >  NC_002947/6173090‑6173381
                                                                                                                                               |                                                                                                                                                     
aaGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACAt                                                                                                                                                   >  6:121258/1‑147 (MQ=255)
 aGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAg                                                                                                                                                >  8:112764/1‑149 (MQ=255)
        attcgaAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCtt                                                                                                                                         >  5:133201/1‑149 (MQ=255)
         ttcgaAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCttt                                                                                                                                        <  1:292296/149‑1 (MQ=255)
           cgaAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTcac                                                                                                                                                     <  8:75677/134‑1 (MQ=255)
           cgaAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTcac                                                                                                                                                     >  7:75677/1‑134 (MQ=255)
            gaAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCgg                                                                                                                                            >  4:16672/1‑142 (MQ=255)
            gaAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCgg                                                                                                                                            <  3:16672/142‑1 (MQ=255)
            gaAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGa                                                                                                                                     >  2:106438/1‑149 (MQ=255)
             aaaGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGAt                                                                                                                                    <  1:106438/149‑1 (MQ=255)
             aaaGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGAt                                                                                                                                    <  6:29599/149‑1 (MQ=255)
              aaGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATg                                                                                                                                   <  7:112764/149‑1 (MQ=255)
              aaGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATg                                                                                                                                   >  3:264694/1‑149 (MQ=255)
              aaGTTAATGGTTGAAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGAt                                                                                                                                    <  2:251503/148‑1 (MQ=255)
                gTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATgg                                                                                                                                  >  4:276755/1‑148 (MQ=255)
                    aTGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAgg                                                                                                                             <  6:133201/149‑1 (MQ=255)
                    aTGGTTGCAATAAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGATTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAg                                                                                                                              <  2:207316/148‑1 (MQ=255)
                     tGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAggg                                                                                                                            <  4:264694/149‑1 (MQ=255)
                            aaTCAAAACCAAAGACCCGCGATTCCGAAAGCTCTCGCTGACCAGTGGATTGATTGCTACTGCTTTGAGGCGCTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTAATCACATAAGTCGGCTTTCGATGGAAAGGGCTCTGaa                                                                                                                     <  4:98959/149‑1 (MQ=255)
                             aTCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATGTGAATAAGTCGGCtt                                                                                                                                         <  5:300456/128‑1 (MQ=255)
                             aTCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATGTGAATAAGTCGGCtt                                                                                                                                         >  6:300456/1‑128 (MQ=255)
                                        aGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGctct                                                                                                         >  4:1167/1‑149 (MQ=255)
                                         gACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCtctt                                                                                                        <  3:1167/149‑1 (MQ=255)
                                          aCGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCtcttc                                                                                                       >  6:323175/1‑149 (MQ=255)
                                              gcgATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGa                                                                                                   <  5:121258/149‑1 (MQ=255)
                                                gATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGAtc                                                                                                 >  8:304628/1‑149 (MQ=255)
                                                     gTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTtctc                                                                                            >  2:252298/1‑149 (MQ=255)
                                                     gTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTtctc                                                                                            >  1:270153/1‑149 (MQ=255)
                                                              tCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGgcgc                                                                                     <  7:304628/147‑1 (MQ=255)
                                                                    aCCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGcc                                                                             <  7:92847/149‑1 (MQ=255)
                                                                      cTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTg                                                                                                    >  1:50071/1‑124 (MQ=255)
                                                                      cTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTg                                                                                                    <  2:50071/124‑1 (MQ=255)
                                                                         tgGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCCCTTCCAGCCCGaaa                                                                        <  8:165236/149‑1 (MQ=255)
                                                                           gATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTc                                                                                  <  6:221888/137‑1 (MQ=255)
                                                                           gATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTc                                                                                  >  5:221888/1‑137 (MQ=255)
                                                                           gATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATg                                                                      >  5:69457/1‑149 (MQ=255)
                                                                            aTTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGt                                                                     >  8:52769/1‑149 (MQ=255)
                                                                                          cTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCa                                                                                >  1:132423/1‑124 (MQ=255)
                                                                                          cTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCa                                                                                <  2:132423/124‑1 (MQ=255)
                                                                                          cTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGgcgcgc                                                        <  3:276755/148‑1 (MQ=255)
                                                                                          cTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCa                                                       >  3:275285/1‑149 (MQ=255)
                                                                                               aGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGc                                                  >  7:279519/1‑149 (MQ=255)
                                                                                               aGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGAGGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGc                                                  >  7:283698/1‑149 (MQ=255)
                                                                                                   ttttttattGGCGCGGCTGCGCCAATAAAAAAGCCGACTTATTCACATAAGTAGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTa                                                >  2:85309/6‑147 (MQ=255)
                                                                                                   ggTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACa                                              <  7:52769/149‑1 (MQ=255)
                                                                                                     tCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGc                                            >  4:229864/1‑149 (MQ=255)
                                                                                                     tCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGc                                            <  6:69457/149‑1 (MQ=255)
                                                                                                       cTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGc                                            >  5:265861/1‑147 (MQ=255)
                                                                                                       cTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCa                                          <  1:252298/149‑1 (MQ=255)
                                                                                                        tataGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAg                                         >  6:190597/1‑149 (MQ=255)
                                                                                                          taGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGt                                       <  5:190597/149‑1 (MQ=255)
                                                                                                            ggCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCAGCATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTAc                                               >  1:215214/1‑140 (MQ=255)
                                                                                                            ggCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTAc                                               <  2:215214/139‑1 (MQ=255)
                                                                                                             gcgcAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCAc                                    >  2:103283/1‑149 (MQ=255)
                                                                                                             gcgcAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCAc                                    >  7:241265/1‑149 (MQ=255)
                                                                                                               gcAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTc                                                                                  >  7:282863/1‑101 (MQ=255)
                                                                                                               gcAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTc                                                                                  <  8:282863/101‑1 (MQ=255)
                                                                                                                    cgcgCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTc                             <  6:265861/149‑1 (MQ=255)
                                                                                                                     gcgcCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGATGGCGCGCACGGGCTACAGCCAGGTCAACCGCTTc                             >  7:240766/1‑148 (MQ=255)
                                                                                                                      cgcCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTa                           <  8:241265/149‑1 (MQ=255)
                                                                                                                                aaaGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGccc                                                                            <  1:277707/90‑1 (MQ=255)
                                                                                                                                aaaGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGccc                                                                            >  2:277707/1‑90 (MQ=255)
                                                                                                                                    ccGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATg             <  2:270153/149‑1 (MQ=255)
                                                                                                                                    ccGACTTATGTGAATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTggg                  >  3:188278/1‑144 (MQ=255)
                                                                                                                                    ccGACTTATGTGAATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTggg                  <  4:188278/144‑1 (MQ=255)
                                                                                                                                      gACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGt           >  3:262360/1‑149 (MQ=255)
                                                                                                                                      gACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGCCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGt           <  1:103283/149‑1 (MQ=255)
                                                                                                                                               acaTAAGTCGGATTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATATTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGAGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGCGTGATGGTCGAGCCGGc  >  7:84080/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
AAGTTTCGATTCGAAAGTTAATGGTTGCAATCAAAACCTAAGACGCGCGATTCGTCAAGCTGTCGCTGACCTGTGGATTGTTTGCTTCTGCTTTGAGGCGGTCCTATAGGCGCAGCCGCGCCAATAAAAAAGCCGACTTATTCACATAAGTCGGCTTTCGATGGAAAGGGCTGTGAATAACTCAGCTCTTCTTGATCTTCTCCGGGCGCTTCCAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGC  >  NC_002947/6173090‑6173381

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: