Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,714,808 C→G 6.9% I56M (ATC→ATG)  hisE → phosphoribosyl‑ATP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,714,8080CG 6.9% 207.4 / 4.7 69I56M (ATC→ATG) hisEphosphoribosyl‑ATP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base G (2/3);  total (32/37)
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

TGCTTGAAGAACGCAAGCAAGCGGCGCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGAACGCCGCTTCGGGCTTTCCGGGCATGATGAAAA  >  NC_002947/5714672‑5714951
                                                                                                                                        |                                                                                                                                               
tGCTTGAAGAACGCAAGCAAGCGGCGCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTAc                                                                                                                                     <  4:105430/149‑1 (MQ=255)
  cTTGAAGAACGCAAGCAAGCGGCGCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACa                                                                                                                                    <  3:188513/148‑1 (MQ=255)
  cTTGAAGAACGCAAGCAAGCGGCGCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATAGCGGCCAAGGATGCCCAGATCAGCAAGGATTACa                                                                                                                                    <  6:267372/148‑1 (MQ=255)
              aagcaagcGGCGCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACTATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTa                                                                                                                       >  2:293713/1‑149 (MQ=255)
              aagcaagcGGCGCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCt                                                                                                                        >  4:340031/1‑148 (MQ=255)
              aagcaagcGGCGCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGAGTACAGCGATGTCATCTa                                                                                                                       >  4:270231/1‑149 (MQ=255)
               agcaagcGGCGCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGAGTACAGCGATGTCATCTAc                                                                                                                      >  2:133178/1‑149 (MQ=255)
                gcaagcGGCGCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTAc                                                                                                                      >  1:20924/1‑148 (MQ=255)
                  aagcGGCGCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGaa                                                                                                                   <  4:79919/149‑1 (MQ=255)
                      ggcgCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCAATGCGGCCAAGCATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAAccg                                                                                                               >  2:329679/1‑149 (MQ=255)
                         gcCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAAccgccg                                                                                                            >  3:297222/1‑149 (MQ=255)
                           cTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCGCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGAc                                                                                                          <  2:103489/149‑1 (MQ=255)
                             gACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCt                                                                                                        >  1:322563/1‑149 (MQ=255)
                              aCAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCtg                                                                                                       >  7:123574/1‑149 (MQ=255)
                              aCAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCt                                                                                                        <  3:340031/148‑1 (MQ=255)
                                aGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCtgtg                                                                                                     >  5:223378/1‑149 (MQ=255)
                                aGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCtgtg                                                                                                     >  6:84548/1‑149 (MQ=255)
                                  ctctTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATGAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGt                                                                                                   >  3:40922/1‑149 (MQ=255)
                                  ctctTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATGAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGt                                                                                                   >  3:36274/1‑149 (MQ=255)
                                     ttACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGttt                                                                                                 >  2:58981/1‑148 (MQ=255)
                                       aCGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCGAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTc                                                                                                <  3:334853/147‑1 (MQ=255)
                                         gTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAg                                                                                            >  2:115895/1‑149 (MQ=255)
                                          tGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGc                                                                                           <  5:84548/149‑1 (MQ=255)
                                              cAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGtt                                                                                                  >  4:64238/1‑138 (MQ=255)
                                              cAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGtt                                                                                                  <  3:64238/138‑1 (MQ=255)
                                              cAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTgg                                                                                       >  2:137128/1‑149 (MQ=255)
                                              cAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTgg                                                                                       <  1:137128/149‑1 (MQ=255)
                                              cAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTAGGAGAAGAGTCGGTCGAAACGATAATTGCGGCCAAGGATGCACAGATAAGCAAGGATTACAGCGATGTCATCTACGAAACAGCCTACCTGATGTTTCATAGCGTgg                                                                                       >  8:95210/1‑149 (MQ=255)
                                               aGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGt                                                                                      <  1:5833/149‑1 (MQ=255)
                                                   tGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGa                                                                                    >  1:156131/1‑147 (MQ=255)
                                                   tGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGa                                                                                    <  2:156131/147‑1 (MQ=255)
                                                   tGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGATCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGCATGCCCAGATCAGCAAGGATTACAGCGATCTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTCGTGATg                                                                                  >  1:192336/1‑149 (MQ=255)
                                                      aCCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGc                                                                                                                                  >  7:194370/1‑98 (MQ=255)
                                                      aCCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGc                                                                                                                                  <  8:194370/98‑1 (MQ=255)
                                                       ccACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCa                                                                              <  1:115895/149‑1 (MQ=255)
                                                        cacaAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAAccgc                                                                                                              >  6:73315/1‑116 (MQ=255)
                                                        cacaAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAAccgc                                                                                                              <  5:73315/116‑1 (MQ=255)
                                                           aaGGGCCTGAACAAGATCCTTGAAAATCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATGAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAgcgc                                                                          <  4:40922/149‑1 (MQ=255)
                                                           aaGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATGAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAgcgc                                                                          <  4:36274/149‑1 (MQ=255)
                                                                ccTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGTCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTgg                                                                     <  3:241394/149‑1 (MQ=255)
                                                                    aaCAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGc                                                                   <  8:123574/147‑1 (MQ=255)
                                                                       aaGATCATTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTACTCTACGAAACCGCCGACCTGTGGTTTGCTAGCCTTGTGATGCTCAGCGCGCTGGGACAgga                                                              >  8:338089/1‑147 (MQ=255)
                                                                             cTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGc                                                         <  6:223378/148‑1 (MQ=255)
                                                                                 aaaaGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGcc                                                    >  5:214106/1‑149 (MQ=255)
                                                                                  aaaGCTCGTCGAAGAGTCGGTCGAAACGTTGATTGCGGCCAAGGATTCGCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGGCGACCTGTGGTTTCATTGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGcg                                                    <  2:192336/148‑2 (MQ=255)
                                                                                  aaaGCTCGGCGAAGATTCGGTCGAAACGATCATTGCGGCCCAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGcc                                                    <  7:95210/148‑1 (MQ=255)
                                                                                     gCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGc                                                >  5:229270/1‑149 (MQ=255)
                                                                                      cTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCt                                               <  6:214106/149‑1 (MQ=255)
                                                                                       tCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTg                                              >  7:256316/1‑149 (MQ=255)
                                                                                         ggCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGa                                            >  3:165148/1‑149 (MQ=255)
                                                                                            gAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGa                                         <  2:20924/149‑1 (MQ=255)
                                                                                             aaGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGaa                                        <  4:165148/149‑1 (MQ=255)
                                                                                              agagTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAAc                                       <  8:256316/149‑1 (MQ=255)
                                                                                              agagTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAAc                                       <  2:63451/149‑1 (MQ=255)
                                                                                                agTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTg                                     <  2:154006/149‑1 (MQ=255)
                                                                                                 gtcggtcgAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTg                                     >  2:86509/1‑148 (MQ=255)
                                                                                                    ggtcgAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGaa                                  <  4:297222/148‑1 (MQ=255)
                                                                                                             gATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGAACGCCGCTTCg                        <  2:322563/149‑1 (MQ=255)
                                                                                                               tCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGAACGCCGCTTgggg                      <  1:58981/149‑1 (MQ=255)
                                                                                                               tCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGAACGCCGCTTCggg                      <  6:323773/149‑1 (MQ=255)
                                                                                                               tCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGAACGCCGCTTCggg                      <  5:7782/149‑1 (MQ=255)
                                                                                                                  ttGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGAACGCCGCTTCGGGCtt                   >  5:333621/1‑149 (MQ=255)
                                                                                                                     cgtccAAGGATGCCCAGATGAGCAAGGTTTACAGCGTTGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTTGATGAACTGGAACGCCGCTTCGGGCTTTcc                <  3:66340/146‑1 (MQ=255)
                                                                                                                       gCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTgg                                                                     >  7:283266/1‑94 (MQ=255)
                                                                                                                       gCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTgg                                                                     <  8:283266/94‑1 (MQ=255)
                                                                                                                       gCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGAACGCCGCTTCGGGCTTTCCgg              <  1:329679/149‑1 (MQ=255)
                                                                                                                       gCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGAACGCCGCTTCGGGCTTTCCgg              >  4:19479/1‑149 (MQ=255)
                                                                                                                                  ccAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGAACGCCGCTTCGGGCTTTCCGGGCATGATGaa    >  3:248013/1‑148 (MQ=255)
                                                                                                                                   cAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGAACGCCGCTTCGGGCTTTCCGGGCATGATGaaaa  <  3:270231/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                               
TGCTTGAAGAACGCAAGCAAGCGGCGCCTGACAGCTCTTACGTGGCCAGCCTGTACCACAAGGGCCTGAACAAGATCCTTGAAAAGCTCGGCGAAGAGTCGGTCGAAACGATCATTGCGGCCAAGGATGCCCAGATCAGCAAGGATTACAGCGATGTCATCTACGAAACCGCCGACCTGTGGTTTCATAGCCTGGTGATGCTCAGCGCGCTGGGCCAGCATCCGCAAGCCGTGCTGGATGAACTGGAACGCCGCTTCGGGCTTTCCGGGCATGATGAAAA  >  NC_002947/5714672‑5714951

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: ‑