Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,315,817 G→A 5.6% A89V (GCC→GTC)  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,8170GA 5.6% 812.5 / 19.9 250A89V (GCC→GTC) phnNribose 1,5‑bisphosphate phosphokinase
Reads supporting (aligned to +/- strand):  ref base G (125/110);  new base A (7/7);  total (133/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 252 total aligned reads displayed.

CACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGCATCGCG  >  NC_000913/4315668‑4315965
                                                                                                                                                     |                                                                                                                                                    
   tGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGt                                                                                                                                                   <  2:12963/150‑1 (MQ=255)
    ggCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTg                                                                                                                                                  <  1:558416/150‑1 (MQ=255)
      ggTTTTCACGGCCACGGGTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGct                                                                                                                                                <  2:1931808/149‑2 (MQ=255)
        ttttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGCCGGATTCGACGTGctgg                                                                                                                                              <  1:1871669/150‑4 (MQ=255)
         tttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGcc                                                                                                                                             <  2:939676/150‑1 (MQ=255)
         tttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGcc                                                                                                                                             <  2:271330/150‑1 (MQ=255)
          ttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGc                                                                                                                                                 >  1:1826955/1‑145 (MQ=255)
          ttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGc                                                                                                                                                 <  2:1826955/145‑1 (MQ=255)
          ttCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCa                                                                                                                                            <  1:261243/150‑1 (MQ=255)
           tCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCag                                                                                                                                           >  2:722385/1‑150 (MQ=255)
            cACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCaga                                                                                                                                          <  1:792934/150‑1 (MQ=255)
            cACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCaga                                                                                                                                          <  2:1154771/150‑1 (MQ=255)
            cACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCaga                                                                                                                                          <  2:222318/150‑1 (MQ=255)
             aCGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCagag                                                                                                                                         <  1:1812353/150‑1 (MQ=255)
              cGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAgaga                                                                                                                                        >  1:992784/1‑150 (MQ=255)
                gCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGCCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATc                                                                                                                                      <  1:745551/150‑1 (MQ=255)
                 ccACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCa                                                                                                                                     >  1:510908/1‑150 (MQ=255)
                    cGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAAtc                                                                                                                                  >  1:638566/1‑150 (MQ=255)
                    cGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAAtc                                                                                                                                  >  1:1174244/1‑150 (MQ=255)
                       ttttCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcg                                                                                                                               <  2:1917031/150‑1 (MQ=255)
                        tttCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcga                                                                                                                              >  2:1380754/1‑150 (MQ=255)
                          tCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacg                                                                                                                            <  2:1711102/150‑1 (MQ=255)
                              ggCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccga                                                                                                                        >  1:411465/1‑150 (MQ=255)
                              ggCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccga                                                                                                                        >  2:625489/1‑150 (MQ=255)
                                 gcTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTcgacgccgacgc                                                                                                                     >  1:1529954/1‑150 (MQ=255)
                                   tGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCa                                                                                                                   <  1:1487605/150‑1 (MQ=255)
                                        gaggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCataata                                                                                                              <  1:1025247/150‑1 (MQ=255)
                                         aggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATAc                                                                                                             <  2:748259/150‑1 (MQ=255)
                                          ggaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACa                                                                                                            >  1:365481/1‑150 (MQ=255)
                                           gaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAg                                                                                                           <  2:239153/150‑1 (MQ=255)
                                           gaTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAg                                                                                                           <  2:413363/150‑1 (MQ=255)
                                            aTCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGa                                                                                                          >  1:98525/1‑150 (MQ=255)
                                                ccGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATAGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAg                                                                                                           >  1:166807/1‑145 (MQ=255)
                                                ccGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATAGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAg                                                                                                           <  2:166807/145‑1 (MQ=255)
                                                 cGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGtt                                                                                                     >  1:113497/1‑150 (MQ=255)
                                                 cGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGtt                                                                                                     <  1:1677283/150‑1 (MQ=255)
                                                   gCGAAACCTGTAAACAGCCGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAg                                                                                                   <  1:1363841/150‑1 (MQ=255)
                                                       aaCCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTg                                                                                               <  1:1380754/150‑1 (MQ=255)
                                                         ccTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGcc                                                                                             >  1:1997574/1‑150 (MQ=255)
                                                            gTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGc                                                                                          >  2:77156/1‑150 (MQ=255)
                                                               aaCAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCTTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTa                                                                                       <  1:156497/150‑1 (MQ=255)
                                                                 cAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAg                                                                                     <  2:551833/150‑1 (MQ=255)
                                                                  aGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAgg                                                                                    <  2:365481/150‑1 (MQ=255)
                                                                  aGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAgg                                                                                    <  2:1207020/150‑1 (MQ=255)
                                                                  aGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAgg                                                                                    <  2:1320084/150‑1 (MQ=255)
                                                                  aGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAgg                                                                                    <  1:1202775/150‑1 (MQ=255)
                                                                     cGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCa                                                                                 >  1:1333184/1‑150 (MQ=255)
                                                                       ggCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGAcc                                                                                                        >  2:1953179/1‑125 (MQ=255)
                                                                            cagcGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGAtt                                                                          >  2:1956770/1‑150 (MQ=255)
                                                                            cagcGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGAtt                                                                          >  1:685387/1‑150 (MQ=255)
                                                                                 ccGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGcc                                                                     >  2:1449072/1‑150 (MQ=255)
                                                                                 ccGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGca                                                                     >  1:1465995/1‑149 (MQ=255)
                                                                                       gATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGCTCAATCTCGACGCCGCCGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCcgcgc                                                               >  2:825950/1‑150 (MQ=255)
                                                                                         tAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCgg                                                             >  2:405970/1‑150 (MQ=255)
                                                                                          aGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCggg                                                            <  1:936682/150‑1 (MQ=255)
                                                                                          aGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCggg                                                            <  2:322396/150‑1 (MQ=255)
                                                                                           gcgcgcCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGt                                                           <  1:722385/150‑1 (MQ=255)
                                                                                            cgcgcccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTa                                                          >  1:971959/1‑150 (MQ=255)
                                                                                            cgcgcccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTa                                                          >  1:920035/1‑150 (MQ=255)
                                                                                                cccgcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTaaaaa                                                      <  2:1467857/150‑1 (MQ=255)
                                                                                                   gcgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACt                                                   <  2:81803/150‑1 (MQ=255)
                                                                                                    cgccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTc                                                  <  1:1220993/150‑1 (MQ=255)
                                                                                                      ccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCt                                                >  2:65189/1‑150 (MQ=255)
                                                                                                      ccTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCt                                                <  2:1633716/150‑1 (MQ=255)
                                                                                                        tGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgc                                              >  1:331587/1‑150 (MQ=255)
                                                                                                         gCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTCAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgct                                             <  2:98525/150‑1 (MQ=255)
                                                                                                          cGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctc                                            >  2:644011/1‑150 (MQ=255)
                                                                                                          cGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctc                                            >  2:591050/1‑150 (MQ=255)
                                                                                                             cAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTgctcgct                                         >  1:418405/1‑150 (MQ=255)
                                                                                                               gATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCa                                       >  1:1958269/1‑150 (MQ=255)
                                                                                                                 tGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAgg                                     >  2:964318/1‑150 (MQ=255)
                                                                                                                   ggCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGc                                   >  1:1734118/1‑150 (MQ=255)
                                                                                                                   ggCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGc                                   >  2:1159255/1‑150 (MQ=255)
                                                                                                                   ggCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCACGCTCAGGGc                                   <  1:1975409/150‑1 (MQ=255)
                                                                                                                    gcgcGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCg                                  <  2:1997574/150‑1 (MQ=255)
                                                                                                                     cgcgTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGCCCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGa                                 >  2:1343267/1‑150 (MQ=255)
                                                                                                                        gTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGatat                              >  1:625015/1‑150 (MQ=255)
                                                                                                                         tGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGAGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATg                             >  2:540002/1‑150 (MQ=255)
                                                                                                                          gAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATgg                            >  1:1927160/1‑150 (MQ=255)
                                                                                                                          gAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATgg                            <  2:446940/150‑1 (MQ=255)
                                                                                                                           aGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGt                           >  1:845191/1‑150 (MQ=255)
                                                                                                                              cGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGtttt                        <  2:1529954/150‑1 (MQ=255)
                                                                                                                               gTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTc                       <  1:688218/150‑1 (MQ=255)
                                                                                                                               gTTGACCAGCAAGTCTAATCCGGCGTCAAGCCACAGATCAAACTCGACTCAGACGCCATAATACAGACCGTTAGCGTGCCATCTTAAGGCCAACAGATTTTGACACGCGCGGGTAAAAAACTCCTGCTCGCTCAGGCCGATATGAtttta                       >  2:1987170/1‑149 (MQ=255)
                                                                                                                                tctggctgCACGCCGGATTCGACGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCa                      >  1:1864360/8‑150 (MQ=255)
                                                                                                                                  gACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACt                    <  1:649987/150‑1 (MQ=255)
                                                                                                                                    ccAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTc                  >  2:1235831/1‑150 (MQ=255)
                                                                                                                                    ccAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTc                  >  2:689746/1‑150 (MQ=255)
                                                                                                                                      tgCACGCCGGATTCGACGTGCAGCCAGAGAACAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCg                >  2:1753461/2‑150 (MQ=255)
                                                                                                                                      aGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATg                             >  1:821129/1‑137 (MQ=255)
                                                                                                                                      aGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCg                >  1:1948875/1‑150 (MQ=255)
                                                                                                                                       gCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCgg               >  1:1011686/1‑150 (MQ=255)
                                                                                                                                        cACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGgc              <  1:757792/150‑1 (MQ=255)
                                                                                                                                          cGTCGAATCCTGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGgcgc            >  1:1792133/1‑150 (MQ=255)
                                                                                                                                          cGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGgcgc            >  1:1346217/1‑150 (MQ=255)
                                                                                                                                            tCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCTTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTg          <  2:1807067/150‑1 (MQ=255)
                                                                                                                                               aaTCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGCa       >  1:286538/1‑150 (MQ=255)
                                                                                                                                               aaTCCGGCGTGCAGACAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCCCGCGGGAAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTGTCACTTCCGGCGCTGGCa       >  2:1856690/1‑150 (MQ=255)
                                                                                                                                                 tCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGCATc     <  2:331587/150‑1 (MQ=255)
                                                                                                                                                  ccGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGCATcg    >  1:4729/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACTGGCGTTTTCACGGCCACGGTTTTCCAGGCGCTGGCGGAGGATCTCCGGCGAAACCTGTAAACAGACGGGCAGCAGCGCCGATTGATAGCGCGCCCGCGCCTGCGGCAGATGGGCGCGTGAGCCGTTGACCAGCACGTCGAATCCGGCGTGCAGCCAGAGATCAATCTCGACGCCGACGCCATAATACAGACCGTTAGCGTGCCAGCTTAAGGCCAACAGATTTTGCCCCGCGCGGGTAAAAAACTCCTGCTCGCTCAGGGCGATATGGTTTTCACTTCCGGCGCTGGCATCGCG  >  NC_000913/4315668‑4315965

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: