Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,827,967 C→G 5.5% intergenic (+169/+30) mtgA → / ← rpoH monofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,827,9670CG 5.5% 263.3 / 2.1 81intergenic (+169/+30)mtgA/rpoHmonofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32
Reads supporting (aligned to +/- strand):  ref base C (37/39);  new base G (2/3);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTCCGC  >  NC_002947/5827825‑5828113
                                                                                                                                                |                                                                                                                                                  
aaCCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGc                                                                                                                                              <  3:164574/149‑1 (MQ=255)
  cccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCCGCCTACAGCCAACCGAATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGATCACCCCGTAGGGCCAGCGGCt                                                                                                                                             >  1:91373/1‑149 (MQ=39)
  cccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCtt                                                                                                                                            <  5:32146/149‑1 (MQ=255)
   ccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCttt                                                                                                                                           >  6:151315/1‑149 (MQ=255)
    cgcgCCTCCGGGCAACCGCATCCCCGGCCTCCTCGCGCGCAACCCCG‑CGCCTACAGGCAACCGCAGCCCCTGT‑AGGAGCGGGCTTACCCGCGCAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCGGGCCCTATGGGGTGAGCGGCtttt                                                                                                                                          <  7:300283/149‑1 (MQ=255)
    cGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCtttt                                                                                                                                          <  6:81450/149‑1 (MQ=255)
      cTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCt                                                                                                                                        <  1:215177/149‑1 (MQ=255)
         cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCttt                                                                                                                                      <  6:186284/148‑1 (MQ=255)
         cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCtt                                                                                                                                       <  6:92198/147‑1 (MQ=255)
         cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                                     <  5:237017/149‑1 (MQ=255)
         cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                                     <  7:301417/149‑1 (MQ=255)
             gggCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTATTTACAAc                                                                                                                                 >  2:381962/1‑149 (MQ=255)
               gCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCACCCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACATccc                                                                                                                               >  7:337583/1‑149 (MQ=255)
                cAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCa                                                                                                                              <  2:63261/149‑1 (MQ=255)
                  aCCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGa                                                                                                                            <  5:372594/149‑1 (MQ=255)
                   ccGCATCCCCGGCCACCTCGCGGGCAAACCCG‑CGCCTACAGGCAAACCCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCACGGGCGTTTCTTTACAACCCAgcc                                                                                                                           >  6:230382/1‑147 (MQ=38)
                            cGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAgg                                                                                                                  >  4:185325/1‑149 (MQ=255)
                            cGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAgg                                                                                                                  <  6:40607/149‑1 (MQ=255)
                                ctcttcGCGGGTAAGCCCG‑CTCCTACAGGGGATGCGGTTGCCTGT‑AGGCGCGGGTTTGCCCGCGAGGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGc                                                                                                              >  8:193363/1‑149 (MQ=17)
                                  cctcGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGTGAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACACCCCAGAGTCAATCAGGCAGc                                                                                                              >  1:74922/1‑148 (MQ=255)
                                      gcgGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATc                                                                                                          <  1:24238/147‑1 (MQ=255)
                                          gCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGCCCAGCGGCTGTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGAGc                                                                                                     >  8:236148/1‑148 (MQ=255)
                                            aaaCCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCtt                                                                                                  <  4:1193/149‑1 (MQ=255)
                                            aaaCCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCtt                                                                                                  >  3:73541/1‑149 (MQ=255)
                                               cccG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCGTAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGAc                                                                                               <  2:91373/149‑1 (MQ=255)
                                                ccG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCt                                                                                                                                        <  1:395722/107‑1 (MQ=255)
                                                ccG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCt                                                                                                                                        >  2:395722/1‑107 (MQ=255)
                                                 cg‑cgCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCACTGATCGCTGCCtgattga                                                                                             <  7:357005/149‑7 (MQ=39)
                                                  g‑cgcCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGc                                                                                                                                                 <  3:319499/96‑1 (MQ=255)
                                                  g‑cgcCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGc                                                                                                                                                 >  4:319499/1‑96 (MQ=255)
                                                  g‑cgcCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCtt                                                                                            <  1:336454/149‑1 (MQ=255)
                                                    cgcCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTg                                                                                           <  6:339280/149‑1 (MQ=255)
                                                    cgcCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTg                                                                                           <  5:53604/149‑1 (MQ=255)
                                                        tACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCa                                                                                       >  4:35868/1‑149 (MQ=255)
                                                          cAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATc                                                                                     <  3:14376/149‑1 (MQ=255)
                                                           aGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCg                                                                                    >  7:83144/1‑149 (MQ=255)
                                                           aGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCg                                                                                    >  3:268795/1‑149 (MQ=255)
                                                           aGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCg                                                                                    >  2:384034/1‑149 (MQ=255)
                                                             gCAACCGCATCCCCTGT‑AGGAGCGGCCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATAGCa                                                                                  >  1:163846/1‑149 (MQ=255)
                                                              cAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAt                                                                                 <  2:74922/149‑1 (MQ=255)
                                                                  cGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttctt                                                                             <  8:337583/149‑1 (MQ=255)
                                                                     aTCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTtctc                                                                          <  5:151315/149‑1 (MQ=255)
                                                                        cccTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCaa                                                                       >  1:285586/1‑149 (MQ=255)
                                                                         ccTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACaa                                                                                                                                  <  3:238820/89‑1 (MQ=255)
                                                                         ccTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACaa                                                                                                                                  >  4:238820/1‑89 (MQ=255)
                                                                          cTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACt                                                                     >  3:328174/1‑149 (MQ=255)
                                                                             t‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                                  <  4:328174/149‑1 (MQ=255)
                                                                                ggAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGa                                                                <  5:287042/149‑1 (MQ=255)
                                                                                           aCCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCa                                                       <  1:381962/147‑1 (MQ=255)
                                                                                            cccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGc                                                     <  1:157541/148‑1 (MQ=255)
                                                                                             ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGcc                                                    >  7:133920/1‑148 (MQ=255)
                                                                                               gcgAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGAt                                                 <  8:19454/149‑1 (MQ=255)
                                                                                                cgAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATa                                                <  1:384034/149‑1 (MQ=255)
                                                                                                  aaGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACg                                              <  8:83144/149‑1 (MQ=255)
                                                                                                       gCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCtt                                                                                                                                            >  1:89818/1‑50 (MQ=255)
                                                                                                       gCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCtt                                                                                                                                            <  2:89818/50‑1 (MQ=255)
                                                                                                          ggCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                                      <  4:268795/149‑1 (MQ=255)
                                                                                                           gCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATAAGCTGTACTTg                                     >  3:104053/1‑149 (MQ=255)
                                                                                                                gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGc                                >  2:72448/1‑149 (MQ=255)
                                                                                                                      aagaaaaGCCGCTGGCCCTATGGGGTGAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTgg                                                                   <  5:69389/105‑1 (MQ=37)
                                                                                                                      aagaaaaGCCGCTGGCCCTATGGGGTGAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTgg                                                                   >  6:69389/4‑108 (MQ=37)
                                                                                                                       aaaaaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCa                         >  5:156387/1‑149 (MQ=255)
                                                                                                                        aaaaaGCCGCTCACCCCATAGGGCCAGCTGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCAt                        <  8:226526/149‑1 (MQ=255)
                                                                                                                         aaaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATg                       >  2:99232/1‑149 (MQ=255)
                                                                                                                         aaaaGCCGCTCACCACATAGGGCGAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATg                       <  4:104053/149‑1 (MQ=255)
                                                                                                                          aaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATg                       >  7:290419/1‑148 (MQ=255)
                                                                                                                           aaGCCGCTGGCCCTATGGGGTGACCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCAGTACTTGTCTGCAAGTTCATg                       >  1:182925/1‑147 (MQ=255)
                                                                                                                             gCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAAt                                                              >  7:123789/1‑106 (MQ=255)
                                                                                                                             gCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAAt                                                              <  8:123789/106‑1 (MQ=255)
                                                                                                                               cGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGt                 <  2:49017/149‑1 (MQ=255)
                                                                                                                                gCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCATCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTg                <  7:193363/149‑1 (MQ=255)
                                                                                                                                    aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCt            >  2:269581/1‑149 (MQ=255)
                                                                                                                                      cccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtc          >  4:324429/1‑149 (MQ=255)
                                                                                                                                      cccATAGGGCCAGCGACTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCTGTTCTTGCCGCGTGGACTtc          >  8:331588/1‑149 (MQ=255)
                                                                                                                                       ccATAGGGCCAGCCGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCCATCAGAGCCTTGACCTTGTTCATCGCATACTTCTACAACTGGCGAATACGCTCAGCCGATACGCTGTACTAGTCTGCCAGTACAAGCAGCGTCCCCTtc          >  3:199798/1‑148 (MQ=255)
                                                                                                                                        cATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtctc        >  5:278370/1‑149 (MQ=255)
                                                                                                                                         aTAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTctct       >  3:17311/1‑149 (MQ=255)
                                                                                                                                           aGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTc     >  8:183534/1‑149 (MQ=255)
                                                                                                                                            gggCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTcc    >  5:201188/1‑149 (MQ=255)
                                                                                                                                            gggCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTcc    <  8:133920/149‑1 (MQ=255)
                                                                                                                                            gggCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTcc    >  7:215985/1‑149 (MQ=255)
                                                                                                                                                cAGCGGCTTTTCTTTACACCCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTCCgc  >  6:357935/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCG‑CGCCTACAGGCAACCGCATCCCCTGT‑AGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTCCGC  >  NC_002947/5827825‑5828113

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑