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)
    agagAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCaa                                                                                                                                              <  2:346904/151‑1 (MQ=255)
         aaTTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTgg                                                                                                                                         >  1:1073357/1‑151 (MQ=255)
            ttGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTg                                                                                                                                       >  1:151330/1‑150 (MQ=255)
              gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCCCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGaa                                                                                                                                     <  2:1395185/150‑1 (MQ=255)
              gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAAt                                                                                                                                    >  3:112691/1‑151 (MQ=255)
               aCGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAAtt                                                                                                                                   <  3:496208/151‑1 (MQ=255)
                cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAAttt                                                                                                                                  >  1:757995/1‑151 (MQ=255)
                cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAAttt                                                                                                                                  <  3:176812/151‑1 (MQ=255)
                  ttGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGt                                                                                                                                >  2:730590/1‑151 (MQ=255)
                    gAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAg                                                                                                                              <  2:1183624/151‑1 (MQ=255)
                       tACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGccc                                                                                                                           >  2:1211959/1‑151 (MQ=255)
                         cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCgg                                                                                                                         >  1:776741/1‑151 (MQ=255)
                         cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCgg                                                                                                                         >  1:860226/1‑151 (MQ=255)
                            cgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGAt                                                                                                                       <  2:1335315/150‑1 (MQ=255)
                            cgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGAt                                                                                                                       >  2:1075711/1‑150 (MQ=255)
                            cgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATa                                                                                                                      <  2:78965/151‑1 (MQ=255)
                             ggATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATa                                                                                                                      >  3:404616/1‑150 (MQ=255)
                               aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAgg                                                                                                                   <  2:1284785/151‑1 (MQ=255)
                                 ggggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCa                                                                                                                 <  1:1111264/151‑1 (MQ=255)
                                  gggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCAt                                                                                                                <  1:1076302/151‑1 (MQ=255)
                                     cACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCa                                                                                                             <  1:1565888/151‑1 (MQ=255)
                                      acttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCAt                                                                                                            >  2:1321851/1‑151 (MQ=255)
                                      acttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCAt                                                                                                            <  2:322444/151‑1 (MQ=255)
                                         tacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTc                                                                                                         >  1:823657/1‑151 (MQ=255)
                                           cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGt                                                                                                          <  2:1421480/148‑1 (MQ=255)
                                           cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcac                                                                                                       >  2:40316/1‑151 (MQ=255)
                                             ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcaca                                                                                                      <  4:404596/150‑1 (MQ=255)
                                             ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcac                                                                                                       <  2:762246/149‑1 (MQ=255)
                                             ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcac                                                                                                       <  2:1718627/149‑1 (MQ=255)
                                             ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACAt                                                                                                     >  3:520395/1‑151 (MQ=255)
                                              tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATc                                                                                                    >  2:564841/1‑151 (MQ=255)
                                                aaGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCg                                                                                                  >  1:785972/1‑151 (MQ=255)
                                                  gATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGc                                                                                                >  2:1520838/1‑151 (MQ=255)
                                                     ccTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAACGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCAt                                                                                              >  1:1447381/1‑150 (MQ=255)
                                                        aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATaaa                                                                                           <  2:1433799/150‑1 (MQ=255)
                                                        aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATaaa                                                                                           <  1:268956/150‑1 (MQ=255)
                                                        aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATaa                                                                                            <  1:1328158/149‑1 (MQ=255)
                                                         aCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACa                                                                                         <  1:515102/151‑1 (MQ=255)
                                                             cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTg                                                                                                                                 >  1:1159437/1‑107 (MQ=255)
                                                             cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTg                                                                                                                                 <  2:1159572/107‑1 (MQ=255)
                                                             cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAac                                                                                     <  4:593066/151‑1 (MQ=255)
                                                              gAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAacg                                                                                    >  4:113300/1‑151 (MQ=255)
                                                                aGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAacgta                                                                                  <  2:1716676/151‑1 (MQ=255)
                                                                    ggCGACGATGAAGATTTTGTCGATGAAGACGTTGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCTGCATAAACAAAacgtacgt                                                                               <  1:961215/150‑1 (MQ=255)
                                                                      cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACgttgt                                                                            >  2:941228/1‑151 (MQ=255)
                                                                      cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGACCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACgttg                                                                             >  1:816051/1‑151 (MQ=255)
                                                                      ccacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAGCCTGGTGAAGTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACgttg                                                                             <  1:1710160/148‑1 (MQ=255)
                                                                       gacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACgttgtt                                                                           <  2:190384/151‑1 (MQ=255)
                                                                         cgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTa                                                                          >  1:1383347/1‑150 (MQ=255)
                                                                         cgaTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTa                                                                          >  2:1385933/1‑150 (MQ=255)
                                                                             gAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGtt                                                                     <  1:1047088/151‑1 (MQ=255)
                                                                              aaGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTc                                                                    >  2:1374965/1‑151 (MQ=255)
                                                                               aGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCg                                                                   <  2:672379/151‑1 (MQ=255)
                                                                                  ttttGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGaa                                                                 <  3:315669/150‑1 (MQ=255)
                                                                                   tttGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAAt                                                                <  3:349712/150‑1 (MQ=255)
                                                                                    ttGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATg                                                               <  4:249224/150‑1 (MQ=255)
                                                                                      gTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGccc                                                            >  3:475176/1‑151 (MQ=255)
                                                                                        cGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcg                                                           <  1:1034970/150‑1 (MQ=255)
                                                                                         gATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcgcg                                                         >  2:1146917/1‑151 (MQ=255)
                                                                                         gATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcgcg                                                         >  1:1414352/1‑151 (MQ=255)
                                                                                         gATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcgcg                                                         >  1:1716283/1‑151 (MQ=255)
                                                                                         gATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcgcg                                                         >  2:775968/1‑151 (MQ=255)
                                                                                            gAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGt                                                      <  2:1528674/151‑1 (MQ=255)
                                                                                             aaGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGtt                                                     <  4:243799/151‑1 (MQ=255)
                                                                                             aaGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGtt                                                     >  1:996315/1‑151 (MQ=255)
                                                                                               gACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTg                                                    >  2:1093498/1‑150 (MQ=255)
                                                                                               gACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTg                                                    >  3:175644/1‑150 (MQ=255)
                                                                                                  gATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                <  2:573410/151‑1 (MQ=255)
                                                                                                  gATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAac                                                 >  1:750892/1‑150 (MQ=255)
                                                                                                   aTGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAacac                                               <  2:1654777/151‑1 (MQ=255)
                                                                                                      aCGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgcg                                            <  2:1322783/151‑1 (MQ=255)
                                                                                                       cGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCgg                                           <  1:603227/151‑1 (MQ=255)
                                                                                                         gAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGc                                         <  1:1022758/151‑1 (MQ=255)
                                                                                                           gTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAt                                       >  2:1101693/1‑151 (MQ=255)
                                                                                                              cGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                    >  1:924047/1‑151 (MQ=255)
                                                                                                               accACTAATTAAGCGTAGCGCATCAGGACCAGGTAGTT‑CTGCAACCTGGTGAATTCGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                    >  1:1165166/2‑150 (MQ=255)
                                                                                                                ccACTAAGTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                     <  4:584829/148‑1 (MQ=255)
                                                                                                                 cACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttt                                 <  2:758100/151‑1 (MQ=255)
                                                                                                                  aCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttt                                 >  4:308394/1‑150 (MQ=255)
                                                                                                                  aCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttttt                                 >  1:119934/1‑150 (MQ=255)
                                                                                                                  aCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                <  3:204806/151‑1 (MQ=255)
                                                                                                                   cTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAGCCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTa                                <  2:1022547/150‑1 (MQ=255)
                                                                                                                     aaTTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTag                               >  3:368975/1‑149 (MQ=255)
                                                                                                                      aTTAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGt                            <  3:191346/151‑1 (MQ=255)
                                                                                                                       ttAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                           >  2:654500/1‑151 (MQ=255)
                                                                                                                       ttAAGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                           >  1:588113/1‑151 (MQ=255)
                                                                                                                          aGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgc                        <  1:1211816/151‑1 (MQ=255)
                                                                                                                          aGCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgc                        >  1:511646/1‑151 (MQ=255)
                                                                                                                           gCGTAGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgcg                       >  3:304590/1‑151 (MQ=255)
                                                                                                                               aGCGCATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAAc                   >  2:1732466/1‑151 (MQ=255)
                                                                                                                                 cgcATCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACAt                 <  4:404728/151‑1 (MQ=255)
                                                                                                                                    aTCAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGACCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGCTGAACACGCGAGACTTTTTTTAGAGTTCGCGCAACATCCg              >  1:114205/1‑151 (MQ=255)
                                                                                                                                      cAGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGc             <  2:396827/150‑1 (MQ=255)
                                                                                                                                       aGGCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc           >  1:1730750/1‑151 (MQ=255)
                                                                                                                                         gCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACtt         <  1:1305696/151‑1 (MQ=255)
                                                                                                                                         gCCCAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACtt         <  2:1696560/151‑1 (MQ=255)
                                                                                                                                           ccAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca       <  1:652146/151‑1 (MQ=255)
                                                                                                                                            cAGGTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac      >  2:790839/1‑151 (MQ=255)
                                                                                                                                               gTAGTT‑CTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAAt   >  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: ‑