Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,403,867 C→T 5.9% A233A (GCG→GCA abgA ← p‑aminobenzoyl‑glutamate hydrolase, A subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,403,8670CT 5.9% 212.2 / 2.4 68A233A (GCG→GCAabgAp‑aminobenzoyl‑glutamate hydrolase, A subunit
Reads supporting (aligned to +/- strand):  ref base C (28/36);  new base T (2/2);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

CCTGCATAACGCCCACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCCGCGCGCCACGCGTACC  >  NC_000913/1403719‑1404013
                                                                                                                                                    |                                                                                                                                                  
cctgcATAACGCCCACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTgc                                                                                                                                                    <  2:197022/149‑1 (MQ=255)
  tgcATAACGCCCACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTgcgc                                                                                                                                                  >  1:307417/1‑149 (MQ=255)
   gcATAACGCCCACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGcgcg                                                                                                                                                 >  3:312454/1‑149 (MQ=255)
           cccACGGTTACGCTGGAAGCTCCTTCGCTGTGAGGGGCGATTGCATGCAGTGCAAGAGGGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTTAGCGGCGGTACCGGTGAAGTGCGCGTCAAAttt                                                                                                                                         <  5:58796/149‑1 (MQ=255)
             cACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGAAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCTTGCGCGACGGTACCGGTGAAGTGCGCGTCAAAttttg                                                                                                                                       >  3:115810/1‑147 (MQ=255)
              aCGTTTACTCCGGCAGCTCCATCGCGGCGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTtggt                                                                                                                                      <  1:318643/149‑1 (MQ=255)
               cGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAATAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTCCGCGTCAAATTtggtg                                                                                                                                     <  7:42763/149‑1 (MQ=255)
                       ctGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGtt                                                                                                                                  >  7:196702/1‑144 (MQ=255)
                       ctGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGtt                                                                                                                                  <  8:196702/144‑1 (MQ=255)
                           aaGCTCCTTCGCTGTGCGGGGCGTTTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTACTGGGTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATaaaa                                                                                                                         <  6:16780/149‑1 (MQ=255)
                           aaGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATaaaa                                                                                                                         <  7:317415/149‑1 (MQ=255)
                           aaGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATaaaa                                                                                                                         <  7:218134/149‑1 (MQ=255)
                                ccTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTAt                                                                                                                     >  2:279604/1‑148 (MQ=255)
                                 cTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTAt                                                                                                                     <  8:325523/147‑1 (MQ=255)
                                 cTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTAt                                                                                                                     <  4:312454/147‑1 (MQ=255)
                                 cTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCa                                                                                                                   <  5:212557/149‑1 (MQ=255)
                                      cTGTGCGGGGCGATTGCATGGAGTGCAAGAGTGGCTTGTGCTTCCGCCAACAAGGAATTGTGAACGTCTTCTGGTTTTGCGCCTGCGTGAGCGTCGGTACCTGTTAAGTTCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGcc                                                                                                              <  5:120684/149‑1 (MQ=255)
                                      cTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGcc                                                                                                              >  6:192970/1‑149 (MQ=255)
                                            ggggCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCTCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCacac                                                                                                        >  3:130772/1‑149 (MQ=255)
                                              ggCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCAcacca                                                                                                      <  7:66906/149‑1 (MQ=255)
                                              ggCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCAcacca                                                                                                      >  4:148466/1‑149 (MQ=255)
                                                     gcatgcaGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGcc                                                                                               >  7:64197/1‑149 (MQ=255)
                                                       atgcaGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGGGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCg                                                                                             <  4:194000/149‑1 (MQ=255)
                                                          caGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAg                                                                                          <  7:19661/149‑1 (MQ=255)
                                                          caGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAg                                                                                          <  1:279604/149‑1 (MQ=255)
                                                          caGTGCAAGAGTGGCTTCTGCGGCCGCCAACAAGCCATTGTGACCGTCTTCAGGAGTTGCGCCTGCGTGCGCGGCGGTACCGGGGAAGTGCGCGTCAAATTGGGTGGGTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAg                                                                                          <  2:94862/149‑1 (MQ=255)
                                                            gTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAgg                                                                                         >  1:293043/1‑148 (MQ=255)
                                                               cAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTAc                                                                                      <  3:148466/148‑1 (MQ=255)
                                                               cAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACg                                                                                     >  7:243213/1‑149 (MQ=255)
                                                                aaGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTCGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCCCGGTGCCAGCAGGTACGc                                                                                    >  6:112234/1‑149 (MQ=255)
                                                                    gTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGt                                                                                <  3:97909/149‑1 (MQ=255)
                                                                            tgctgcCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAAtgt                                                                        <  8:64197/149‑1 (MQ=255)
                                                                               tgcCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCa                                                                     >  3:84278/1‑149 (MQ=255)
                                                                                  cgccAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACgg                                                                  >  7:68756/1‑149 (MQ=255)
                                                                                   gccAACAACGCATTGTGCCCGACTTCTGGTTTTGCGCCTCCGTGAGCGGCGGTACAGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTCCCGCCCACCACCGTGCCCGCACTTACGCCAGTGCCAATGTGCAccg                                                                  >  6:260529/1‑146 (MQ=255)
                                                                                    ccAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCAcacc                                                                                                       <  8:193457/110‑1 (MQ=255)
                                                                                    ccAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCAcacc                                                                                                       >  7:193457/1‑110 (MQ=255)
                                                                                    ccAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCCGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCa                                                                <  2:45416/149‑1 (MQ=255)
                                                                                      aacaaGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCGAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAg                                                               >  5:82451/1‑148 (MQ=255)
                                                                                          aGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTaaaa                                                          >  8:170201/1‑149 (MQ=255)
                                                                                          aGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTaaaa                                                          >  8:165479/1‑149 (MQ=255)
                                                                                               ttGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATc                                                     <  2:183309/149‑1 (MQ=255)
                                                                                                 gtgACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCa                                                    <  5:260167/148‑1 (MQ=255)
                                                                                                  tgACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCaa                                                   >  1:66799/1‑148 (MQ=255)
                                                                                                  tgACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCaa                                                   <  5:113139/148‑1 (MQ=255)
                                                                                                  tgACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAAc                                                  <  2:66799/149‑1 (MQ=255)
                                                                                                   gACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAAca                                                 >  6:334556/1‑149 (MQ=255)
                                                                                                      cGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAAcatca                                              >  5:153738/1‑149 (MQ=255)
                                                                                                       gTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAAcatcat                                             >  1:177860/1‑149 (MQ=255)
                                                                                                        tcttctGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAgg                                                                                         >  4:267262/1‑104 (MQ=255)
                                                                                                        tcttctGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAgg                                                                                         <  3:267262/104‑1 (MQ=255)
                                                                                                          ttctGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCCTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGACAGTa                                                             <  2:81398/130‑1 (MQ=255)
                                                                                                          ttctGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCCTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGACAGTa                                                             >  1:81398/1‑130 (MQ=255)
                                                                                                            ctGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCTGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTAc                                         <  4:115810/148‑1 (MQ=255)
                                                                                                            ctGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTAc                                         <  2:177860/148‑1 (MQ=255)
                                                                                                               gTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGaca                                     <  8:243213/149‑1 (MQ=255)
                                                                                                                ttttGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGTGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGacac                                    >  8:140818/1‑149 (MQ=255)
                                                                                                                ttttGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGTGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGacac                                    >  2:303305/1‑149 (MQ=255)
                                                                                                                  ttGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACAcc                                   <  1:342676/148‑1 (MQ=255)
                                                                                                                  ttGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGGTGCCATATAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACAcc                                   <  7:23040/148‑1 (MQ=255)
                                                                                                                   tGCGCCTTCGTTAGCGGCGGTACAGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCt                                  <  4:130772/148‑1 (MQ=255)
                                                                                                                   tGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCt                                  <  2:307417/148‑1 (MQ=255)
                                                                                                                           cGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGa                           <  5:192970/147‑1 (MQ=255)
                                                                                                                           cGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGAc                          >  1:50136/1‑148 (MQ=255)
                                                                                                                           cGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCAAACCACGGTGCCCTCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGa                           <  4:84278/147‑1 (MQ=255)
                                                                                                                                    cggTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCcgc                >  6:159848/1‑149 (MQ=255)
                                                                                                                                      gTACCGGTGAAGTGTGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCcgcgc              <  7:140818/149‑1 (MQ=255)
                                                                                                                                      gTACCGGTGAAGTGTGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCcgcgc              <  1:303305/149‑1 (MQ=255)
                                                                                                                                            gTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCCGCGCGCCAcg        >  7:96690/1‑149 (MQ=255)
                                                                                                                                            gTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCCGCGCGCCAcg        >  1:135787/1‑149 (MQ=255)
                                                                                                                                             tGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCCGCGCGCCAcgc       >  7:168674/1‑149 (MQ=255)
                                                                                                                                                  ggcgcGTCAAATTTGGTGGGTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCCGCGCGCCACGCGTAcc  <  5:131300/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCTGCATAACGCCCACGTTTACTCTGGAAGCTCCTTCGCTGTGCGGGGCGATTGCATGCAGTGCAAGAGTGGCTTGTGCTGCCGCCAACAAGGCATTGTGACCGTCTTCTGGTTTTGCGCCTGCGTGAGCGGCGGTACCGGTGAAGTGCGCGTCAAATTTGGTGGTTGCCATAAAATTATCACTGCCGCACACCACGGTGCCCGCAGGTACGCCAGTGCCAATGTGCACGGCAGTAAAATAATCAACATCATCTACGACACCTGCATCGACCATCGCCCGCGCGCCACGCGTACC  >  NC_000913/1403719‑1404013

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: