Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,239,103 A→C 7.9% intergenic (‑78/‑297) argF ← / → PP_1080 ornithine carbamoyltransferase/molybdopterin‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,239,1030AC 7.9% 176.5 / 5.5 63intergenic (‑78/‑297)argF/PP_1080ornithine carbamoyltransferase/molybdopterin‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base C (2/3);  total (31/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 = 7.61e-01

CGTGGCGGATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTTGCTACAGA  >  NC_002947/1238956‑1239251
                                                                                                                                                     |                                                                                                                                                    
cGTGGCGGATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAg                                                                                                                                                     <  4:347246/149‑1 (MQ=255)
 gTGGCGGATCACCCTGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGc                                                                                                                                                    <  4:12504/149‑1 (MQ=255)
      ggATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CTTACGCCGCag                                                                                                                                               <  4:112005/149‑2 (MQ=255)
        aTCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CTTACGCCGCag                                                                                                                                               <  8:217056/147‑2 (MQ=255)
        aTCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CTTACGCCGCag                                                                                                                                               >  7:217056/1‑146 (MQ=255)
           aCCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGc                                                                                                                                          <  2:16325/149‑1 (MQ=255)
             cccGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTGTGTGACAGCGGGGGTGCCCTGCGG‑‑CGGAAGCCGGCGGGGGCgg                                                                                                                                        >  2:332789/1‑149 (MQ=255)
              ccGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAACGTGCCTAGCGCTCATGAGTGACTACCTTATCTCCAACGGAATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGa                                                                                                                                       <  6:292832/149‑1 (MQ=255)
              ccGAGCAATGCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCGTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGGGACGTGCGG‑‑CGTAAGCCGCACGGGGCGGa                                                                                                                                       <  6:292827/149‑1 (MQ=255)
               cGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGAc                                                                                                                                      >  5:214232/1‑149 (MQ=255)
                gAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGAc                                                                                                                                      >  2:101988/1‑148 (MQ=255)
                     aTTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGaaa                                                                                                                                  >  5:207689/1‑147 (MQ=255)
                     aTTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGaaa                                                                                                                                  >  2:220335/1‑147 (MQ=255)
                            ggtggtGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCTG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAg                                                                                                                         <  1:220335/149‑1 (MQ=255)
                              tggtgAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGaaa                                                                                                                                  <  3:321394/138‑1 (MQ=255)
                              tggtgAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGaaa                                                                                                                                  >  4:321394/1‑138 (MQ=255)
                                gtgAAGTCCAGCAGGGAGAGATAGTGCCTAGCGCTCATGATTGACTACCTTATCTGCAACTGTATGCAGGTCGACCGTATGGTTTTTTTTGAGAACGGGGGTGACCTGCTG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGgcg                                                                                                                      <  7:78824/148‑1 (MQ=255)
                                gtgAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGgcg                                                                                                                      >  2:222027/1‑148 (MQ=255)
                                   aaGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCaa                                                                                                                   <  1:299801/148‑1 (MQ=255)
                                     gTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATAc                                                                                                                <  7:354356/149‑1 (MQ=255)
                                     gTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATAc                                                                                                                >  8:133754/1‑149 (MQ=255)
                                        ccgcagGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACtat                                                                                                             <  6:214232/147‑1 (MQ=255)
                                              ggAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGTAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTa                                                                                                        >  4:257133/1‑148 (MQ=255)
                                                agagaAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATa                                                                                                                                >  6:192158/1‑122 (MQ=255)
                                                agagaAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATa                                                                                                                                <  5:192158/122‑1 (MQ=255)
                                                 gagaAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATg                                                                                                    >  5:45516/1‑149 (MQ=255)
                                                  agaAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGTGGTGACCTGCGG‑‑CGTAAGCAGAAAGGGGCGGACGAAATAGTGAAAGGCGAAATACTATAAGTAAATg                                                                                                    >  5:345076/1‑148 (MQ=255)
                                                     aaGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATTCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCCCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTc                                                                                                  <  1:101988/149‑1 (MQ=255)
                                                       gTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCt                                                                                              >  1:164291/1‑149 (MQ=255)
                                                                 cTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAgg                                                                                    >  6:153502/1‑149 (MQ=255)
                                                                     ttaGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGTTTTTTTTTTACAACGGGGGTTACCTGCGG‑‑CTTAAGCCGCACGGGGCGTACTAAATAGTGAAAGGCGCAATATTATAATTATATGTCGCCTGAAGCCAAgggggg                                                                                <  7:234796/147‑1 (MQ=255)
                                                                      gAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAgggg                                                                                  >  1:411/1‑146 (MQ=255)
                                                                        gTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGaa                                                                              >  8:150955/1‑148 (MQ=255)
                                                                             tACCTTTTATGCAACGGTAGGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGCCATGCGG‑‑CGTAAGCCGCCGGGGGCGGCGGAAATAGGGAAGGGGGCATTACTTTATGTAAATGTCGCTGGAGGCAAGGGGGGGAAAAgggg                                                                        >  2:344650/1‑149 (MQ=255)
                                                                             tACCTTATCTGCAACGGTATGCAGGTCGCCCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCGGAAGCCAAGGGGGGAAAAgcgg                                                                        >  8:149155/1‑149 (MQ=255)
                                                                             tACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAgcgg                                                                        >  3:120875/1‑149 (MQ=255)
                                                                               ccTAATCTGCAACGGTATGCAGGTCGACCGTATGCTTTTTTTTTACAACCGGGGTGACCTGCGG‑‑CGTAAGCAGCACGGGGAGGACGAAATAGGGAAAGGCGCAAAACTATAAGTAAATCTCGCCTGAAGCCAAGGCGGGAAAAgcggcg                                                                      >  6:187399/1‑149 (MQ=255)
                                                                                 ttATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGt                                                                     >  2:63610/1‑148 (MQ=255)
                                                                                 ttATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGAGGCGGACGAAATAGGGAAAg                                                                                                                         >  3:256318/1‑96 (MQ=255)
                                                                                        gttatgTATGCAGGTCGACCGTATGGTTTTTTTTGACATCGGGGGTGACATGCGG‑‑CGTAAGCCGCACGGGGCGTACGAAATAGGGAAAGGCGCAATACTATAAGTTAATGTCTCCTGATGCCAAGGGGGGAAAAGCGGCGTGAttgcaa                                                             <  1:72808/144‑1 (MQ=255)
                                                                                              tATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAACGCGCAATACTATAAGTAAATGTCGACTGAAGCAAAGTGGGCAAAAGCGGCGTGAttgcaa                                                             >  2:72808/1‑143 (MQ=255)
                                                                                                ttcatgTCTACGGTTTGGTTTTTATTGACAACGGGGGTCAACTGCGG‑‑CGTAAGCCGCACTGGGGGGACTAAATAGGGAAAGCCGCAATACTATAAGTAAATGTCGCCTGAATCGAAGGGGGGAAAAGCGGCGTGAttgcatagtgtgtt                                                     <  5:187399/144‑1 (MQ=255)
                                                                                                     gTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgttgcaaa                                                <  1:222027/149‑1 (MQ=255)
                                                                                                     gTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCATGACGCCAAGGGGTGAAAAGCGGCGTGAttgcaaagggtgttgct                                                  >  4:189644/1‑146 (MQ=255)
                                                                                                        gACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgttgcaaagg                                              <  3:257133/148‑1 (MQ=255)
                                                                                                        gACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgttgcaaagg                                              <  7:133754/148‑1 (MQ=255)
                                                                                                         aCCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgttgcaaagggt                                            <  6:207689/149‑1 (MQ=255)
                                                                                                                ggTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTa                                     >  6:13031/1‑149 (MQ=255)
                                                                                                                           aCAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGccc                          <  4:61780/149‑1 (MQ=255)
                                                                                                                           aCAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGccc                          >  5:67876/1‑149 (MQ=255)
                                                                                                                             aaCGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagg                                                          <  4:217647/115‑1 (MQ=255)
                                                                                                                             aaCGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagg                                                          >  3:217647/1‑115 (MQ=255)
                                                                                                                                  gggTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGt                   >  1:37611/1‑149 (MQ=255)
                                                                                                                                     tGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCt                <  1:332789/149‑1 (MQ=255)
                                                                                                                                      gACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCtt               >  5:51151/1‑149 (MQ=255)
                                                                                                                                       ccgtGCGG‑‑CTTACGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTgg                             >  7:279539/4‑134 (MQ=255)
                                                                                                                                       ccgtGCGG‑‑CTTACGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTgg                             <  8:279539/131‑1 (MQ=255)
                                                                                                                                        ccTGCGG‑‑CGTAAGCCGCACTGGGCGGACGAAATAGGTAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCttt              <  6:345076/148‑1 (MQ=255)
                                                                                                                                            cgg‑‑cgTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTAtt           <  6:45516/147‑1 (MQ=255)
                                                                                                                                              g‑‑cgTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTTGCt       <  1:63610/149‑1 (MQ=255)
                                                                                                                                                 cgTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTTGCTa      <  5:153502/149‑1 (MQ=255)
                                                                                                                                                 cgTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTTGCTa      <  1:344650/149‑1 (MQ=255)
                                                                                                                                                 cgTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGGCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTTGCTa      <  7:149155/149‑1 (MQ=255)
                                                                                                                                                     aGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTTGCTACAGa  <  2:164291/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGTGGCGGATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGG‑‑CGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTTGCTACAGA  >  NC_002947/1238956‑1239251

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