Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,537,952 A→T 6.5% intergenic (+14/+19) PP_1348 → / ← pduO MutT/nudix family protein/thiamine‑phosphate pyrophosphorylase/ATP:cob(I)alamin adenosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,537,9520AT 6.5% 189.8 / 2.9 61intergenic (+14/+19)PP_1348/pduOMutT/nudix family protein/thiamine‑phosphate pyrophosphorylase/ATP:cob(I)alamin adenosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base T (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

GCAGCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTT  >  NC_002947/1537806‑1538095
                                                                                                                                                  |                                                                                                                                               
gcAGCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCa                                                                                                                                                 >  1:257606/1‑147 (MQ=255)
 cAGCGGTTGACTGCTGGGGTGAACAAGCCGGTTTACCTGCTGGGGGGGGTCGGGCCGAGTGAGCGCGAGCCGGCCTGGGAAGCCGGGGCACAGCGTGTAGCCGGGATCCGGGCGGTCTGGCCCGAGGTTTGAGTCACGGGGCTGCAAag                                                                                                                                              <  2:207281/149‑1 (MQ=255)
   gCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAagca                                                                                                                                            <  5:217574/149‑1 (MQ=255)
      gTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGcc                                                                                                                                         >  4:140960/1‑149 (MQ=255)
      gTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGcc                                                                                                                                         <  2:51236/149‑1 (MQ=255)
         gACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCTTTGCAGcccc                                                                                                                                       <  8:37583/148‑1 (MQ=255)
            tGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCtt                                                                                                                                    >  7:8096/1‑148 (MQ=255)
            tGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCtt                                                                                                                                    >  7:139406/1‑148 (MQ=255)
              ccGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAagca                                                                                                                                            <  1:257782/138‑1 (MQ=255)
              ccGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAagca                                                                                                                                            >  2:257782/1‑138 (MQ=255)
              ccGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCttt                                                                                                                                   <  5:191393/147‑1 (MQ=255)
              ccGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCtt                                                                                                                                    <  8:139406/146‑1 (MQ=255)
              ccGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTa                                                                                                                                  <  5:62122/148‑1 (MQ=255)
              ccGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTa                                                                                                                                  <  2:257606/148‑1 (MQ=255)
               cGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAg                                                                                                                                 >  4:99354/1‑148 (MQ=255)
                gggTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTa                                                                                                                               >  2:22181/1‑149 (MQ=255)
                  gTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATc                                                                                                                             <  3:140960/149‑1 (MQ=255)
                      aacaaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGc                                                                                                                          <  5:5289/148‑1 (MQ=255)
                      aacaaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGCGAAGCCGGGGCACAGGGTCTAGCCGGGATCAGCGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGcc                                                                                                                         >  5:135338/1‑149 (MQ=255)
                       acaaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCt                                                                                                                        >  4:247929/1‑149 (MQ=255)
                         aaGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGAACAGGGTGTAGCCGGGATCCGGGCGTTCAGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTATAGTATCAGCCtt                                                                                                                       >  3:49387/1‑148 (MQ=255)
                           gCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAg                                                                                                                    >  1:164400/1‑149 (MQ=255)
                           gCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGAGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAg                                                                                                                    >  1:168679/1‑149 (MQ=255)
                                      tgctgGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTa                                                                                                                                  >  5:75763/1‑124 (MQ=255)
                                      tgctgGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTa                                                                                                                                  <  6:75763/124‑1 (MQ=255)
                                        ctgGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGgggcggg                                                                                                       <  2:164400/149‑1 (MQ=255)
                                           ggtgggGTAGGGCCGTGTGAGGGCGTGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGACGGGATCCGGTCGTTCTGGCCCGAGGTTTTAGTCATGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCt                                                                                                     <  6:135338/148‑1 (MQ=255)
                                           ggtgggGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGgggcgggc                                                                                                      <  3:99354/147‑1 (MQ=255)
                                             tgggGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTTTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGcc                                                                                                  <  5:66484/149‑1 (MQ=255)
                                             tgggGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGc                                                                                                   >  3:122406/1‑148 (MQ=255)
                                                      gCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACtt                                                                                         <  8:8096/149‑1 (MQ=255)
                                                       ccGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGCGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTc                                                                                        >  2:219911/1‑149 (MQ=255)
                                                          aGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGc                                                                                     <  4:122406/149‑1 (MQ=255)
                                                              aGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTTGGGGCGGGCTGCCACAACACTTCTGCAAt                                                                                  <  1:22181/148‑1 (MQ=255)
                                                              aGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAAt                                                                                  <  3:247929/148‑1 (MQ=255)
                                                                        ggCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTgccgcc                                                                       <  2:168679/149‑1 (MQ=255)
                                                                           cTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTAt                                                                                                                              <  1:234348/91‑1 (MQ=255)
                                                                           cTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTAt                                                                                                                              >  2:234348/1‑91 (MQ=255)
                                                                           cTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTg                                                                                                               <  5:168164/106‑1 (MQ=255)
                                                                           cTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTg                                                                                                               >  6:168164/1‑106 (MQ=255)
                                                                              ggAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTgccgccgc                                                                     >  4:79855/1‑145 (MQ=255)
                                                                                   ccGGGGCTCAGGGTGTAGCGGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCATAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACATCACTTCTGCAATGCCTTGCCGCCGCCCGATGATc                                                            <  4:49387/149‑1 (MQ=255)
                                                                                   ccGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATc                                                            >  8:125750/1‑149 (MQ=255)
                                                                                          acaGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGtgctg                                                     >  5:220434/1‑149 (MQ=255)
                                                                                           caGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGtgctgc                                                    <  3:79855/149‑1 (MQ=255)
                                                                                                      cGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAagag                                           >  7:153846/1‑147 (MQ=255)
                                                                                                          aTCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTc                                     <  1:219911/149‑1 (MQ=255)
                                                                                                                 cGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGc                              <  7:138201/149‑1 (MQ=255)
                                                                                                                        gCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCACTCAGGTTTGGGGGCGGGCTGCCACAACACTTATGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGa                       >  1:203027/1‑149 (MQ=255)
                                                                                                                          ccGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGata                     >  5:144290/1‑149 (MQ=255)
                                                                                                                           cGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGata                     >  4:200969/1‑148 (MQ=255)
                                                                                                                               gTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGc                <  7:125750/149‑1 (MQ=255)
                                                                                                                                  tGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGt                                                         >  7:20377/1‑105 (MQ=255)
                                                                                                                                  tGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGt                                                         <  8:20377/105‑1 (MQ=255)
                                                                                                                                     taaagggggCTGCTTTGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTATGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGAc           >  6:196354/6‑148 (MQ=255)
                                                                                                                                     gTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGAc           >  6:17460/1‑148 (MQ=255)
                                                                                                                                      aaagggggCTGCTTTGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAagag                                           >  6:121054/5‑115 (MQ=255)
                                                                                                                                      aaagggggCTGCTTTGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAagag                                           <  5:121054/111‑1 (MQ=255)
                                                                                                                                          ggggCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGc     >  6:245994/1‑149 (MQ=255)
                                                                                                                                            ggCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATAc                  >  7:194904/1‑134 (MQ=255)
                                                                                                                                            ggCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATAc                  <  8:194904/134‑1 (MQ=255)
                                                                                                                                             gCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTAGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCtt  >  2:93585/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCAGCGGTTGACTGCCGGGTTCAACAAGCCGGTTTACCTGCTGGGTGGGGTAGGGCCGAGTGAGCGCGAGCAGGCCTGGGAAGCCGGGGCACAGGGTGTAGCCGGGATCCGGGCGTTCTGGCCCGAGGTTTGAGTCAAGGGGCTGCAAAGCAGCCCCCTTTAGTATCAGCCTTCAGGTTTGGGGGCGGGCTGCCACAACACTTCTGCAATGCCTTGCCGCCGCCCGATGATCCGTGCTGCCACAAAGAGCAGGTCGGAAAGCCGGTTGATATACGCCAGACCAACGCCTT  >  NC_002947/1537806‑1538095

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: