Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,851,736 T→A 7.7% intergenic (+83/+36) apt → / ← recR adenine phosphoribosyltransferase/recombination protein RecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,851,7360TA 7.7% 190.6 / 7.3 65intergenic (+83/+36)apt/recRadenine phosphoribosyltransferase/recombination protein RecR
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base A (3/2);  total (34/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 = 2.58e-01

GCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAGC  >  NC_002947/4851587‑4851884
                                                                                                                                                     |                                                                                                                                                    
gcgTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACTAAAGCCGCGGATGTTTTTTAGGCGCGGGtt                                                                                                                                                      >  1:938557/1‑150 (MQ=255)
             gCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATgg                                                                                                                                         <  2:589222/150‑1 (MQ=255)
                       ggCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                                                                               >  2:924524/2‑150 (MQ=255)
                            aCCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTAGTGGGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGccc                                                                                                                          <  2:568989/150‑1 (MQ=255)
                               ttCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCaa                                                                                                                       <  2:488196/150‑1 (MQ=255)
                                tCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAACCCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAAt                                                                                                                      <  1:137536/150‑1 (MQ=255)
                                tCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAAt                                                                                                                      <  1:437404/150‑1 (MQ=255)
                                  tGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGa                                                                                                                    <  1:924524/150‑1 (MQ=255)
                                      tGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCa                                                                                                                <  1:442527/150‑1 (MQ=255)
                                          cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTCCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCGGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGc                                                                                                            >  1:581227/1‑150 (MQ=255)
                                                 tCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATgg                                                                                                     >  2:50816/1‑150 (MQ=255)
                                                 tCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATgg                                                                                                     <  1:298550/150‑1 (MQ=255)
                                                   ccTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCCGAGCGAAATGGGc                                                                                                   >  2:129056/1‑150 (MQ=255)
                                                   ccTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGGAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGc                                                                                                   >  1:613655/1‑150 (MQ=255)
                                                    cTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGGATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGcc                                                                                                  >  1:720564/1‑150 (MQ=255)
                                                     tCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGTCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGAAGCTTAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                                                                                 <  1:699317/150‑1 (MQ=39)
                                                        gCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTAAACCcgcgc                                                                                                                                             <  1:56168/103‑2 (MQ=255)
                                                        gCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTAAACCcgcgc                                                                                                                                             >  2:56168/1‑102 (MQ=255)
                                                          gAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGcc                                                                                            >  1:1027285/1‑150 (MQ=255)
                                                             tATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGGAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATGAGAGCGAAATGGGCCGGCGCCCTg                                                                                         >  1:193868/1‑150 (MQ=255)
                                                              aTTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCaga                                                                                                              >  1:294629/1‑128 (MQ=255)
                                                              aTTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCaga                                                                                                              <  2:294629/128‑1 (MQ=255)
                                                                   gTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGa                                                                                                                    >  2:856854/1‑117 (MQ=255)
                                                                   gTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGa                                                                                                                    <  1:856854/117‑1 (MQ=255)
                                                                    tGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGgcc                                                                                                                                      >  2:411957/1‑98 (MQ=255)
                                                                    tGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGgcc                                                                                                                                      <  1:411957/98‑1 (MQ=255)
                                                                    tGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                                                                               <  2:105309/105‑1 (MQ=255)
                                                                    tGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                                                                               >  1:105309/1‑105 (MQ=255)
                                                                      ggCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                                                                               <  1:579222/103‑1 (MQ=255)
                                                                      ggCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                                                                               >  2:579222/1‑103 (MQ=255)
                                                                      ggCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCAt                                                                                <  2:1027285/150‑1 (MQ=255)
                                                                      ggCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCAt                                                                                >  2:680839/1‑150 (MQ=255)
                                                                       gCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                                                                               <  1:1135200/150‑1 (MQ=255)
                                                                        tttttGGAGCTGCGTTGCAGCCCATTCGCGGGTAACCCCAGTCGTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATgg                                                                              <  1:616286/149‑1 (MQ=39)
                                                                        cTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATgg                                                                              <  2:505860/150‑1 (MQ=255)
                                                                        cTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATgg                                                                              <  1:904234/150‑1 (MQ=255)
                                                                         ttttGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATggg                                                                             >  1:830052/1‑150 (MQ=255)
                                                                               ggCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAggg                                                                       >  1:312932/1‑150 (MQ=255)
                                                                                     tttGCAGCCCATTAGATTATAAACCCGCTCCTACGAAAGCCGCGTATGTGTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATTAGAGAGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAccgc                                                                 >  1:947769/1‑150 (MQ=38)
                                                                                          aGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcc                                                            <  2:520801/150‑1 (MQ=255)
                                                                                           gcgcATTCGCGGCTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGAGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGACCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcca                                                           >  2:550370/4‑150 (MQ=255)
                                                                                           gCCCATTCGCGGGTAACCCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcca                                                           <  1:1049899/150‑1 (MQ=255)
                                                                                           gCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcca                                                           >  1:507827/1‑150 (MQ=255)
                                                                                            cccATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTGTGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTccac                                                          >  2:697409/1‑150 (MQ=255)
                                                                                                   gcgGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                                                                                 <  2:1142262/104‑1 (MQ=255)
                                                                                                   gcgGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                                                                                 >  1:1142262/1‑104 (MQ=255)
                                                                                                       gTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCt                                               >  2:1059901/1‑150 (MQ=255)
                                                                                                       gTAAACCCGCTCCTACGAAAGCCGCGTATGATTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCt                                               >  2:536108/1‑150 (MQ=255)
                                                                                                            cccGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTcgccg                                          <  2:610179/150‑1 (MQ=255)
                                                                                                             ccGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTcgccgc                                         >  1:516623/1‑150 (MQ=255)
                                                                                                                    tACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTgg                                  <  2:613655/150‑1 (MQ=255)
                                                                                                                    tACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTgg                                  >  2:821507/1‑150 (MQ=255)
                                                                                                                     aCGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGc                                 >  1:1017640/1‑150 (MQ=255)
                                                                                                                       gAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGcac                               >  1:728903/1‑150 (MQ=255)
                                                                                                                                     ttttGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGc                 >  2:17991/1‑150 (MQ=255)
                                                                                                                                       ttGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCga               <  2:312932/150‑1 (MQ=255)
                                                                                                                                       ttGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCga               <  1:821507/150‑1 (MQ=255)
                                                                                                                                         gTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCgaga             <  1:284858/150‑1 (MQ=255)
                                                                                                                                         ccattcgcgGGTAAACCCGCGCATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCaa                    <  2:438102/138‑1 (MQ=255)
                                                                                                                                         ccattcgcgGGTAAACCCGCGCATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCaa                    >  1:438102/6‑143 (MQ=255)
                                                                                                                                            ttcgcgGGTAAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGc          >  1:954318/3‑150 (MQ=255)
                                                                                                                                            ggCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGc          >  2:389392/1‑150 (MQ=255)
                                                                                                                                                  ggTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCaa    <  2:1017640/150‑1 (MQ=255)
                                                                                                                                                   gTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGGGGCCTAGTGGCACACCATGGGCAATGCGAGACGCCACCAAg   >  1:968310/1‑150 (MQ=255)
                                                                                                                                                   gTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCgaga             <  2:249535/140‑1 (MQ=255)
                                                                                                                                                   gTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCgaga             >  1:249535/1‑140 (MQ=255)
                                                                                                                                                    tttACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTg                                   <  2:1160833/117‑1 (MQ=255)
                                                                                                                                                    tttACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTg                                   >  1:1160833/1‑117 (MQ=255)
                                                                                                                                                    tttACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAGc  <  1:17991/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAGC  >  NC_002947/4851587‑4851884

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: