Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,402,087 C→G 5.9% G390A (GGC→GCC)  abgB ← p‑aminobenzoyl‑glutamate hydrolase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,402,0870CG 5.9% 691.8 / 24.9 203G390A (GGC→GCC) abgBp‑aminobenzoyl‑glutamate hydrolase subunit B
Reads supporting (aligned to +/- strand):  ref base C (107/83);  new base G (10/2);  total (118/85)
Fisher's exact test for biased strand distribution p-value = 7.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 207 total aligned reads displayed.

TTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGACTGTTTTGCCGATCG  >  NC_000913/1401938‑1402236
                                                                                                                                                     |                                                                                                                                                     
 tgcggcCATAGTTTTCGCCGCCAGCAGCATTCGTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACTGCAAAACATTGGCTGAAACAATGGGCAACAGGCAGTTTCCAACTGACGTCGct                                                                                                                                                      <  2:795568/150‑2 (MQ=255)
  gctgcCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCg                                                                                                                                                     >  1:970960/1‑150 (MQ=255)
  gctgcCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCg                                                                                                                                                     >  1:1956899/1‑150 (MQ=255)
   ctgcCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGa                                                                                                                                                    <  1:841152/150‑1 (MQ=255)
     gcCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACa                                                                                                                                                  >  1:985631/1‑150 (MQ=255)
           gTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTc                                                                                                                                            <  1:1986492/150‑1 (MQ=255)
                cgccgcCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGc                                                                                                                                       <  1:1763238/150‑1 (MQ=255)
                  ccgccAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCg                                                                                                                                     <  2:1366658/150‑1 (MQ=255)
                  ccgccAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCg                                                                                                                                     >  2:585738/1‑150 (MQ=255)
                   cgccAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGc                                                                                                                                    <  1:157633/150‑1 (MQ=255)
                          aggaTTCCTTTATGAGCAATAGTTGTTCGCCGCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACTGCAAAACAGGGGCTGAAACACTGGGCAACTGGCAGTTTCCAACTGACGTCGCCGACATCAGGCGATGCCGCAAGCACg                                                                                                                             <  2:1300604/147‑1 (MQ=255)
                          agcaTTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACg                                                                                                                             >  1:443049/1‑150 (MQ=255)
                           gcaTTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGt                                                                                                                            >  2:1623460/1‑150 (MQ=255)
                           gcaTTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGt                                                                                                                            <  1:1623460/150‑1 (MQ=255)
                                ccTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCg                                                                                                                       >  2:1942201/1‑150 (MQ=255)
                                ccTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCg                                                                                                                       >  2:562925/1‑150 (MQ=255)
                                 cTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCgg                                                                                                                      <  2:1573814/150‑1 (MQ=255)
                                  tttATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACTGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGt                                                                                                                     >  2:488043/1‑150 (MQ=255)
                                   ttATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTg                                                                                                                    <  1:633347/150‑1 (MQ=255)
                                    tATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTgg                                                                                                                   >  1:1392604/1‑150 (MQ=255)
                                     aTGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGTGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTggc                                                                                                                  <  2:985631/150‑1 (MQ=255)
                                     aTGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTggc                                                                                                                  >  1:605263/1‑150 (MQ=255)
                                         gCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCCACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCa                                                                                                              >  2:1609231/1‑150 (MQ=255)
                                         gCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCa                                                                                                              >  2:1745387/1‑150 (MQ=255)
                                                 tGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGc                                                                                                      >  2:62194/1‑150 (MQ=255)
                                                 tGTTCGCCCCTGGCTAACCAGTTGCCACGTACGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGc                                                                                                      <  1:452487/150‑1 (MQ=255)
                                                  gTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCg                                                                                                     >  1:679888/1‑150 (MQ=255)
                                                   ttCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGa                                                                                                    >  2:1014175/1‑150 (MQ=255)
                                                      gCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCa                                                                                                              <  1:50779/137‑1 (MQ=255)
                                                       ccccTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTc                                                                                                >  2:1893035/1‑150 (MQ=255)
                                                            ggCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCa                                                                                                              <  1:1603166/131‑1 (MQ=255)
                                                              cTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGcc                                                                                         >  2:741750/1‑150 (MQ=255)
                                                                aaCCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGc                                                                                                                                    <  1:200893/105‑1 (MQ=255)
                                                                aaCCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGAc                                                                                                   >  1:348366/1‑138 (MQ=255)
                                                                  ccAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCACCTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTa                                                                                     >  2:1286525/1‑150 (MQ=255)
                                                                  ccAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAg                                                                                       <  2:920537/148‑1 (MQ=255)
                                                                        gCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGttt                                                                               >  2:1161842/1‑150 (MQ=255)
                                                                         ccACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTc                                                                              >  1:126268/1‑150 (MQ=255)
                                                                           aCGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCCGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCAc                                                                            >  2:236616/1‑150 (MQ=255)
                                                                               aTGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGcc                                                                                         >  2:828231/1‑133 (MQ=255)
                                                                               aTGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATg                                                                        <  2:679888/150‑1 (MQ=255)
                                                                               aTGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATg                                                                        >  2:68661/1‑150 (MQ=255)
                                                                               aTGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATg                                                                        <  2:443067/150‑1 (MQ=255)
                                                                                tGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGa                                                                       >  1:1034917/1‑150 (MQ=255)
                                                                                   aGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGt                                                                    <  1:1745387/150‑1 (MQ=255)
                                                                                   aGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGt                                                                    >  2:348065/1‑150 (MQ=255)
                                                                                    gCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTa                                                                   >  1:998628/1‑150 (MQ=255)
                                                                                      ggTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGCCTTCATTCGCCAGTACCGTTTCACGATGACGTAGt                                                                 >  1:818729/1‑150 (MQ=255)
                                                                                      ggTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGt                                                                 >  2:35137/1‑150 (MQ=255)
                                                                                       gtgtACCGACGGCAAAACAGGGGCTGAAACACTGGGCATCAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTg                                                                <  2:947302/150‑1 (MQ=255)
                                                                                         gtACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCa                                                              >  2:1277411/1‑150 (MQ=255)
                                                                                            ccGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCaaaa                                                           >  1:951449/1‑150 (MQ=255)
                                                                                            ccGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCaaaa                                                           >  1:757927/1‑150 (MQ=255)
                                                                                              gACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAAc                                                         >  2:184758/1‑150 (MQ=255)
                                                                                               aCGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAAcc                                                        <  1:741750/150‑1 (MQ=255)
                                                                                                cGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCt                                                       >  1:187080/1‑150 (MQ=255)
                                                                                                cGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCt                                                       >  1:53739/1‑150 (MQ=255)
                                                                                                 ggCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGTGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCtt                                                      >  2:20470/1‑150 (MQ=255)
                                                                                                 ggCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGAGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCtt                                                      >  1:1062065/1‑150 (MQ=255)
                                                                                                  gCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTg                                                     >  1:1809820/1‑150 (MQ=255)
                                                                                                  gCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTCATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTg                                                     >  2:1809218/1‑150 (MQ=255)
                                                                                                   cAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGc                                                    >  2:1470509/1‑150 (MQ=255)
                                                                                                    aaaaCAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGcc                                                   <  1:585738/150‑1 (MQ=255)
                                                                                                           gggCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCg                                            <  1:1803813/150‑1 (MQ=255)
                                                                                                            ggCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGc                                           >  1:1355941/1‑150 (MQ=255)
                                                                                                             gCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGcc                                          <  2:1187739/150‑1 (MQ=255)
                                                                                                             gCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGcc                                          >  1:913312/1‑150 (MQ=255)
                                                                                                              cTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCa                                         >  1:1169761/1‑150 (MQ=255)
                                                                                                                gAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCAcc                                       >  2:1108737/1‑150 (MQ=255)
                                                                                                                  aaCACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCgg                                     >  1:1220099/1‑150 (MQ=255)
                                                                                                                      cTGTGCAACAGTCAGTTTCGAACTGAGGTCGCCTACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTgc                                 <  2:421737/150‑1 (MQ=255)
                                                                                                                      cTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTgc                                 <  1:707714/150‑1 (MQ=255)
                                                                                                                        gggCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTgcgg                               >  1:376402/1‑150 (MQ=255)
                                                                                                                         ggCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCACATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTgcggc                              >  1:1432387/1‑150 (MQ=255)
                                                                                                                          gCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTgcggcg                             <  2:1993889/150‑1 (MQ=255)
                                                                                                                            aaCAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGat                           <  1:236616/150‑1 (MQ=255)
                                                                                                                              cAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGatat                         >  2:969832/1‑150 (MQ=255)
                                                                                                                              cAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGatat                         <  2:1787952/150‑1 (MQ=255)
                                                                                                                                ggCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGcc                                          >  2:1664072/1‑131 (MQ=255)
                                                                                                                                ggCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGAtatta                       >  1:1114745/1‑150 (MQ=255)
                                                                                                                                ggCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGAtatta                       <  1:924257/150‑1 (MQ=255)
                                                                                                                                ggCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGAtatta                       >  2:264470/1‑150 (MQ=255)
                                                                                                                                ggCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGAtatta                       >  2:1525041/1‑150 (MQ=255)
                                                                                                                                 gCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGAtattat                      <  1:327864/150‑1 (MQ=255)
                                                                                                                                    gTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCa                   <  1:1196395/150‑1 (MQ=255)
                                                                                                                                     tttCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGTGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAg                  <  1:20470/150‑1 (MQ=255)
                                                                                                                                     tttCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCaa                  >  2:1309231/1‑149 (MQ=255)
                                                                                                                                     tttCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAg                  <  2:320329/150‑1 (MQ=255)
                                                                                                                                     tttCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAg                  <  2:647876/150‑1 (MQ=255)
                                                                                                                                       tCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGAc                <  1:316589/150‑1 (MQ=255)
                                                                                                                                         cgaCTGATGTCGGCGACGTCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGACTg              >  1:791813/3‑150 (MQ=255)
                                                                                                                                          gaCTGATGTCGGCGACGTCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGAc                <  1:578124/146‑1 (MQ=255)
                                                                                                                                          gaCTGATGTCGGCGACGTCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGAc                >  2:578124/2‑147 (MQ=255)
                                                                                                                                          gaCTGATGTCGGCGACGTCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGACTGt             >  2:1079106/2‑150 (MQ=255)
                                                                                                                                            cTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGACTGttt           >  1:1640352/1‑150 (MQ=255)
                                                                                                                                                cGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGACTGTTTTGcc       <  2:1392604/150‑1 (MQ=255)
                                                                                                                                                cGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGACTGTTTTGcc       <  2:1961912/150‑1 (MQ=255)
                                                                                                                                                 gTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGACTGTTTTGCCg      <  2:1220099/150‑1 (MQ=255)
                                                                                                                                                  tCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGACTGTTTTGCCGa     >  1:499925/1‑150 (MQ=255)
                                                                                                                                                   cGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGACTGTTTTGCCGAt    >  1:115721/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATATGGAGCGACTTCATTCGCCAGTACCGTTTCACGATGACGTAGTGCAAAAACCTTGCCGTTTTCGCCACCGGTTGCGGCGATATTATTCAGACTGTTTTGCCGATCG  >  NC_000913/1401938‑1402236

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