Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,478,921 Δ1 bp intergenic (+32/+113) rraB → / ← yjgM protein inhibitor of RNase E/GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,478,9200C.99.3% 710.7 / ‑1.5 145intergenic (+31/+114)rraB/yjgMprotein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (67/77);  total (68/77)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

Only 100 of 148 total aligned reads displayed.

TGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATC  >  NC_000913/4478771‑4479066
                                                                                                                                                     |                                                                                                                                                  
 ggCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTc                                                                                                                                                    <  1:1626340/149‑1 (MQ=255)
    agagAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCaa                                                                                                                                              <  2:346904/151‑1 (MQ=255)
         aaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTgg                                                                                                                                         >  1:1073357/1‑151 (MQ=255)
            ttGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTg                                                                                                                                       >  1:151330/1‑150 (MQ=255)
              gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCCCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGaa                                                                                                                                     <  2:1395185/150‑1 (MQ=255)
              gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAAt                                                                                                                                    >  3:112691/1‑151 (MQ=255)
               aCGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAAtt                                                                                                                                   <  3:496208/151‑1 (MQ=255)
                cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAAttt                                                                                                                                  >  1:757995/1‑151 (MQ=255)
                cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAAttt                                                                                                                                  <  3:176812/151‑1 (MQ=255)
                  ttGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGt                                                                                                                                >  2:730590/1‑151 (MQ=255)
                    gAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAg                                                                                                                              <  2:1183624/151‑1 (MQ=255)
                       tACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGccc                                                                                                                           >  2:1211959/1‑151 (MQ=255)
                         cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCgg                                                                                                                         >  1:776741/1‑151 (MQ=255)
                         cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCgg                                                                                                                         >  1:860226/1‑151 (MQ=255)
                            cgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGAt                                                                                                                       <  2:1335315/150‑1 (MQ=255)
                            cgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGAt                                                                                                                       >  2:1075711/1‑150 (MQ=255)
                            cgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATa                                                                                                                      <  2:78965/151‑1 (MQ=255)
                             ggATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATa                                                                                                                      >  3:404616/1‑150 (MQ=255)
                               aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAgg                                                                                                                   <  2:1284785/151‑1 (MQ=255)
                                 ggggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCa                                                                                                                 <  1:1111264/151‑1 (MQ=255)
                                  gggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCAt                                                                                                                <  1:1076302/151‑1 (MQ=255)
                                     cACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCa                                                                                                             <  1:1565888/151‑1 (MQ=255)
                                      acttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCAt                                                                                                            >  2:1321851/1‑151 (MQ=255)
                                      acttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCAt                                                                                                            <  2:322444/151‑1 (MQ=255)
                                         tacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTc                                                                                                         >  1:823657/1‑151 (MQ=255)
                                           cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGt                                                                                                          <  2:1421480/148‑1 (MQ=255)
                                           cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcac                                                                                                       >  2:40316/1‑151 (MQ=255)
                                             ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcaca                                                                                                      <  4:404596/150‑1 (MQ=255)
                                             ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcac                                                                                                       <  2:762246/149‑1 (MQ=255)
                                             ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcac                                                                                                       <  2:1718627/149‑1 (MQ=255)
                                             ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACAt                                                                                                     >  3:520395/1‑151 (MQ=255)
                                              tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATc                                                                                                    >  2:564841/1‑151 (MQ=255)
                                                aaGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCg                                                                                                  >  1:785972/1‑151 (MQ=255)
                                                  gATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGc                                                                                                >  2:1520838/1‑151 (MQ=255)
                                                     ccTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAACGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCAt                                                                                              >  1:1447381/1‑150 (MQ=255)
                                                        aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATaaa                                                                                           <  2:1433799/150‑1 (MQ=255)
                                                        aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATaaa                                                                                           <  1:268956/150‑1 (MQ=255)
                                                        aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATaa                                                                                            <  1:1328158/149‑1 (MQ=255)
                                                         aCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACa                                                                                         <  1:515102/151‑1 (MQ=255)
                                                             cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTg                                                                                                                                 >  1:1159437/1‑107 (MQ=255)
                                                             cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTg                                                                                                                                 <  2:1159572/107‑1 (MQ=255)
                                                             cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAac                                                                                     <  4:593066/151‑1 (MQ=255)
                                                              gAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAacg                                                                                    >  4:113300/1‑151 (MQ=255)
                                                                aGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAacgta                                                                                  <  2:1716676/151‑1 (MQ=255)
                                                                    ggCGACGATGAAGATTTTGTCGATGAAGACGTTGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCTGCATAAACAAAacgtacgt                                                                               <  1:961215/150‑1 (MQ=255)
                                                                      cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACgttgt                                                                            >  2:941228/1‑151 (MQ=255)
                                                                      cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGACCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACgttg                                                                             >  1:816051/1‑151 (MQ=255)
                                                                      ccacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAGCCTGGTGAAGTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACgttg                                                                             <  1:1710160/148‑1 (MQ=255)
                                                                       gacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACgttgtt                                                                           <  2:190384/151‑1 (MQ=255)
                                                                         cgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTa                                                                          >  1:1383347/1‑150 (MQ=255)
                                                                         cgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTa                                                                          >  2:1385933/1‑150 (MQ=255)
                                                                             gAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGtt                                                                     <  1:1047088/151‑1 (MQ=255)
                                                                              aaGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTc                                                                    >  2:1374965/1‑151 (MQ=255)
                                                                               aGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCg                                                                   <  2:672379/151‑1 (MQ=255)
                                                                                  ttttGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGaa                                                                 <  3:315669/150‑1 (MQ=255)
                                                                                   tttGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAAt                                                                <  3:349712/150‑1 (MQ=255)
                                                                                    ttGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATg                                                               <  4:249224/150‑1 (MQ=255)
                                                                                      gTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGccc                                                            >  3:475176/1‑151 (MQ=255)
                                                                                        cGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcg                                                           <  1:1034970/150‑1 (MQ=255)
                                                                                         gATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcgcg                                                         >  2:1146917/1‑151 (MQ=255)
                                                                                         gATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcgcg                                                         >  1:1414352/1‑151 (MQ=255)
                                                                                         gATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcgcg                                                         >  1:1716283/1‑151 (MQ=255)
                                                                                         gATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcgcg                                                         >  2:775968/1‑151 (MQ=255)
                                                                                            gAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGt                                                      <  2:1528674/151‑1 (MQ=255)
                                                                                             aaGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGtt                                                     <  4:243799/151‑1 (MQ=255)
                                                                                             aaGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGtt                                                     >  1:996315/1‑151 (MQ=255)
                                                                                               gACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTg                                                    >  2:1093498/1‑150 (MQ=255)
                                                                                               gACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTg                                                    >  3:175644/1‑150 (MQ=255)
                                                                                                  gATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                <  2:573410/151‑1 (MQ=255)
                                                                                                  gATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAac                                                 >  1:750892/1‑150 (MQ=255)
                                                                                                   aTGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAacac                                               <  2:1654777/151‑1 (MQ=255)
                                                                                                      aCGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgcg                                            <  2:1322783/151‑1 (MQ=255)
                                                                                                       cGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCgg                                           <  1:603227/151‑1 (MQ=255)
                                                                                                         gAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGc                                         <  1:1022758/151‑1 (MQ=255)
                                                                                                           gTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                       >  2:1101693/1‑151 (MQ=255)
                                                                                                              cGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                    >  1:924047/1‑151 (MQ=255)
                                                                                                               accACTAATTAAGCGTAGCGCATCAGGACCAGGTAGTTCTGCAACCTGGTGAATTCGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                    >  1:1165166/2‑150 (MQ=255)
                                                                                                                ccACTAAGTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                     <  4:584829/148‑1 (MQ=255)
                                                                                                                 cACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttt                                 <  2:758100/151‑1 (MQ=255)
                                                                                                                  aCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttt                                 >  4:308394/1‑150 (MQ=255)
                                                                                                                  aCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttt                                 >  1:119934/1‑150 (MQ=255)
                                                                                                                  aCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                <  3:204806/151‑1 (MQ=255)
                                                                                                                   cTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAGCCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                <  2:1022547/150‑1 (MQ=255)
                                                                                                                     aaTTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTag                               >  3:368975/1‑149 (MQ=255)
                                                                                                                      aTTAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGt                            <  3:191346/151‑1 (MQ=255)
                                                                                                                       ttAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                           >  2:654500/1‑151 (MQ=255)
                                                                                                                       ttAAGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                           >  1:588113/1‑151 (MQ=255)
                                                                                                                          aGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgc                        <  1:1211816/151‑1 (MQ=255)
                                                                                                                          aGCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgc                        >  1:511646/1‑151 (MQ=255)
                                                                                                                           gCGTAGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgcg                       >  3:304590/1‑151 (MQ=255)
                                                                                                                               aGCGCATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAAc                   >  2:1732466/1‑151 (MQ=255)
                                                                                                                                 cgcATCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACAt                 <  4:404728/151‑1 (MQ=255)
                                                                                                                                    aTCAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGACCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGCTGAACACGCGAGACTTTTTTTAGAGTTCGCGCAACATCCg              >  1:114205/1‑151 (MQ=255)
                                                                                                                                      cAGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGc             <  2:396827/150‑1 (MQ=255)
                                                                                                                                       aGGCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc           >  1:1730750/1‑151 (MQ=255)
                                                                                                                                         gCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACtt         <  1:1305696/151‑1 (MQ=255)
                                                                                                                                         gCCCAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACtt         <  2:1696560/151‑1 (MQ=255)
                                                                                                                                           ccAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca       <  1:652146/151‑1 (MQ=255)
                                                                                                                                            cAGGTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac      >  2:790839/1‑151 (MQ=255)
                                                                                                                                               gTAGTTCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAAt   >  2:828285/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATC  >  NC_000913/4478771‑4479066

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: