Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,315,811 A→G 5.8% F91S (TTC→TCC)  phnN ← ribose 1,5‑bisphosphate phosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,315,8110AG 5.8% 745.2 / 23.9 240F91S (TTC→TCC) phnNribose 1,5‑bisphosphate phosphokinase
Reads supporting (aligned to +/- strand):  ref base A (120/106);  new base G (7/7);  total (127/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

Only 100 of 246 total aligned reads displayed.

AATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGCA  >  NC_000913/4315663‑4315960
                                                                                                                                                    |                                                                                                                                                     
 aTTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAAt                                                                                                                                                     <  1:778106/150‑1 (MQ=255)
       cTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAACCCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCTGCACGCCGGATTCGACg                                                                                                                                               <  1:686256/150‑1 (MQ=255)
        tGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGt                                                                                                                                              <  2:12963/150‑1 (MQ=255)
            gTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGctg                                                                                                                                          <  1:135476/150‑3 (MQ=255)
             ttttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGctgg                                                                                                                                         <  1:1871669/150‑4 (MQ=255)
             ttttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGctgg                                                                                                                                         <  1:1490972/150‑4 (MQ=255)
              tttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGcc                                                                                                                                        <  2:939676/150‑1 (MQ=255)
              tttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGcc                                                                                                                                        >  1:98465/1‑150 (MQ=255)
               ttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGc                                                                                                                                            <  2:1826955/145‑1 (MQ=255)
               ttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGc                                                                                                                                            >  1:1826955/1‑145 (MQ=255)
               ttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCa                                                                                                                                       >  2:1975409/1‑150 (MQ=255)
               ttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCa                                                                                                                                       >  1:1807067/1‑150 (MQ=255)
                tCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCag                                                                                                                                      >  2:722385/1‑150 (MQ=255)
                 cACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCaga                                                                                                                                     <  2:1154771/150‑1 (MQ=255)
                 cACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCaga                                                                                                                                     >  2:287328/1‑150 (MQ=255)
                  aCGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCagag                                                                                                                                    >  1:1169615/1‑150 (MQ=255)
                     gCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATc                                                                                                                                 >  1:26593/1‑150 (MQ=255)
                     gCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATc                                                                                                                                 <  1:1343950/150‑1 (MQ=255)
                      ccACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCa                                                                                                                                >  1:510908/1‑150 (MQ=255)
                         cGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAAtc                                                                                                                             <  1:634973/150‑1 (MQ=255)
                         cGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAAtc                                                                                                                             <  1:641840/150‑1 (MQ=255)
                            ttttCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcg                                                                                                                          <  1:66784/150‑1 (MQ=255)
                               tCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacg                                                                                                                       <  2:1711102/150‑1 (MQ=255)
                                ccAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgc                                                                                                                      >  2:1068293/1‑150 (MQ=255)
                                 cAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgcc                                                                                                                     >  1:339605/1‑150 (MQ=255)
                                  aGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccg                                                                                                                    >  1:287253/1‑150 (MQ=255)
                                      gcTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccgacgc                                                                                                                >  1:1529954/1‑150 (MQ=255)
                                       cTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCgacgccgacgcc                                                                                                               <  1:1367130/150‑1 (MQ=255)
                                         ggCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCat                                                                                                             >  1:81803/1‑150 (MQ=255)
                                          gCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAACCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCata                                                                                                            >  1:412058/1‑150 (MQ=255)
                                             gaggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCataata                                                                                                         <  1:1025247/150‑1 (MQ=255)
                                              aggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATAc                                                                                                        <  2:748259/150‑1 (MQ=255)
                                              aggaGCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATAc                                                                                                        <  2:827143/150‑1 (MQ=255)
                                                 aTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGa                                                                                                     <  2:1239206/150‑1 (MQ=255)
                                                 aTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGa                                                                                                     >  1:98525/1‑150 (MQ=255)
                                                   ctcCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGAcc                                                                                                   >  1:1925559/1‑150 (MQ=255)
                                                     ccGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGt                                                                                                 >  1:1093775/1‑150 (MQ=255)
                                                     ccGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATAGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAg                                                                                                      <  2:166807/145‑1 (MQ=255)
                                                      cGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGtt                                                                                                <  1:1677283/150‑1 (MQ=255)
                                                       ggCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTa                                                                                               <  2:218579/150‑1 (MQ=255)
                                                       ggCGAAACCACTAAACAGACGGGCAGCAGCCCCGATTGATAGCGCGCCCGCGCCTGCCCCAGATGGGCGCGTGACCCCTTGCCCAGCACGTCGAATCCGCCGTGCAGCCAGAGATCAATCTCGCCGCCGACGCCAAAATACAGACCGATa                                                                                               >  2:582273/1‑150 (MQ=255)
                                                         cGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGc                                                                                             <  2:1925559/150‑1 (MQ=255)
                                                         cGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGc                                                                                             <  1:72675/150‑1 (MQ=255)
                                                          gAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCg                                                                                            <  1:906301/150‑1 (MQ=255)
                                                          gAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCg                                                                                            <  1:1978042/150‑1 (MQ=255)
                                                            aaCCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTg                                                                                          <  2:1956742/150‑1 (MQ=255)
                                                              ccTGTAAACAGACGGGCAGCAGCGCCGATTGATATCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGcc                                                                                        >  2:1993466/1‑150 (MQ=255)
                                                              ccTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGcc                                                                                        >  1:1487120/1‑150 (MQ=255)
                                                              ccTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGcc                                                                                        <  2:1174244/150‑1 (MQ=255)
                                                               cTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGctg                                                                                                                                          >  1:1855024/1‑97 (MQ=255)
                                                                tGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAg                                                                                      >  1:752121/1‑150 (MQ=255)
                                                                  tAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCt                                                                                    <  1:210689/150‑1 (MQ=255)
                                                                       aGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAgg                                                                               <  2:541941/150‑1 (MQ=255)
                                                                       aGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAgg                                                                               <  2:1207020/150‑1 (MQ=255)
                                                                        gACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGc                                                                              <  1:923417/150‑1 (MQ=255)
                                                                            ggCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAAc                                                                          <  1:1962443/150‑1 (MQ=255)
                                                                             gcagcaGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGCCCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACa                                                                         >  2:235379/1‑150 (MQ=255)
                                                                                gcagcGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGAt                                                                      >  1:446940/1‑150 (MQ=255)
                                                                                   gcgcCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGAtttt                                                                   >  1:125244/1‑150 (MQ=255)
                                                                                   gcgcCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGAtttt                                                                   >  1:1610127/1‑150 (MQ=255)
                                                                                      ccGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGcc                                                                >  2:1449072/1‑150 (MQ=255)
                                                                                      ccGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGcc                                                                >  2:68929/1‑150 (MQ=255)
                                                                                       cGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGccc                                                               <  2:632069/150‑1 (MQ=255)
                                                                                            gATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcgc                                                          >  2:797906/1‑150 (MQ=255)
                                                                                             aTAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGCCCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcgcg                                                         >  1:121648/1‑150 (MQ=255)
                                                                                               aGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCggg                                                       <  2:322396/150‑1 (MQ=255)
                                                                                                     cccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaaa                                                 <  2:1467857/150‑1 (MQ=255)
                                                                                                     cccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaaa                                                 >  1:55406/1‑150 (MQ=255)
                                                                                                        gcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACt                                              <  2:81803/150‑1 (MQ=255)
                                                                                                           ccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCt                                           <  2:1633716/150‑1 (MQ=255)
                                                                                                              gCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgct                                        <  2:1333184/150‑1 (MQ=255)
                                                                                                              gCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTCAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgct                                        <  2:98525/150‑1 (MQ=255)
                                                                                                                          cgcgTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGCCCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGa                            >  2:1343267/1‑150 (MQ=255)
                                                                                                                          cgcgTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGa                            >  2:649987/1‑150 (MQ=255)
                                                                                                                           gcgTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGat                           >  2:1935794/1‑150 (MQ=255)
                                                                                                                           gcgTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGat                           <  1:1955966/150‑1 (MQ=255)
                                                                                                                           gcgTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGat                           >  1:1174771/1‑150 (MQ=255)
                                                                                                                             gTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGatat                         >  1:625015/1‑150 (MQ=255)
                                                                                                                              tGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCACCGCGCGGGTAAAAAACTCCTTCTCGCTCAGGGCGATATg                        >  2:613002/1‑150 (MQ=255)
                                                                                                                                   cGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGtttt                   <  2:1529954/150‑1 (MQ=255)
                                                                                                                                    gTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGTGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTc                  <  2:510908/150‑1 (MQ=255)
                                                                                                                                    gTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTc                  <  2:971959/150‑1 (MQ=255)
                                                                                                                                     ttGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCa                 <  1:1665826/150‑1 (MQ=255)
                                                                                                                                         gctgCACGCCGGATTCGACGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTc             >  1:1721397/4‑150 (MQ=255)
                                                                                                                                          cAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTcc            <  1:964318/150‑1 (MQ=255)
                                                                                                                                           tgCACGCCGGATTCGACGTGCAGCCAGAGAACAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCg           >  2:1753461/2‑150 (MQ=255)
                                                                                                                                           atcACGTCGAATCACGCGAGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCg           >  2:1424456/3‑150 (MQ=255)
                                                                                                                                           aGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATg                        >  1:821129/1‑137 (MQ=255)
                                                                                                                                           aGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCg           >  2:440496/1‑150 (MQ=255)
                                                                                                                                           aGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCg           >  2:350538/1‑150 (MQ=255)
                                                                                                                                           aGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCg           >  1:375856/1‑150 (MQ=255)
                                                                                                                                           aGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCg           >  1:1948875/1‑150 (MQ=255)
                                                                                                                                            gCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCgg          >  1:139125/1‑150 (MQ=255)
                                                                                                                                            gCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCgg          >  1:1011686/1‑150 (MQ=255)
                                                                                                                                               cGTCGAATCCTGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGgcgc       >  1:1792133/1‑150 (MQ=255)
                                                                                                                                               cGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGgcgc       >  1:698745/1‑150 (MQ=255)
                                                                                                                                                 tCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCTTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTg     <  2:1807067/150‑1 (MQ=255)
                                                                                                                                                    aaTCCGGCGTGCAGCCAGAGATCAATCTCTACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGCa  <  2:375856/150‑1 (MQ=255)
                                                                                                                                                    aaTCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGCa  <  2:412058/150‑1 (MQ=255)
                                                                                                                                                    aaTCCGGCGTGCAGACAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCCCGCGGGAAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTGTCACTTCCGGCGCTGGCa  >  2:1856690/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AATTTCACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGCA  >  NC_000913/4315663‑4315960

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: ‑