Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,805,351 G→C 7.2% R200R (CGC→CGG)  waaP ← lipopolysaccharide core heptose (I) kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,805,3510GC 7.2% 915.4 / 44.0 272R200R (CGC→CGG) waaPlipopolysaccharide core heptose (I) kinase
Reads supporting (aligned to +/- strand):  ref base G (132/119);  new base C (2/18);  total (134/138)
Fisher's exact test for biased strand distribution p-value = 2.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 274 total aligned reads displayed.

CATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTA  >  NC_000913/3805207‑3805498
                                                                                                                                                |                                                                                                                                                   
 aTGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGgcgc                                                                                                                                               >  1:224555/1‑150 (MQ=255)
  tGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGcgcg                                                                                                                                              <  1:1010160/150‑1 (MQ=255)
   gTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCgg                                                                                                                                             <  1:1690887/150‑1 (MQ=255)
   gTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCgg                                                                                                                                             >  2:1021933/1‑150 (MQ=255)
    ttCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGt                                                                                                                                            <  2:312545/150‑1 (MQ=255)
     tCAATAAATTCTGTTCAAGAGACTATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCGCTGGCGCGATa                                                                                                                                           <  1:1737819/150‑1 (MQ=255)
      cAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTAc                                                                                                                                          >  2:503864/1‑150 (MQ=255)
      cAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTAc                                                                                                                                          >  2:1153836/1‑150 (MQ=255)
       aaTAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCGCTGGCGCGATaaa                                                                                                                                         <  2:714953/150‑3 (MQ=255)
       aaTAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCGCTGGCGCGATaaa                                                                                                                                         <  1:495692/150‑3 (MQ=255)
        aTAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCGCTGGCGCGATaaag                                                                                                                                        <  2:1349558/150‑4 (MQ=255)
        aTAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAAGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCGCTGGCGCGATaaag                                                                                                                                        <  2:478359/150‑4 (MQ=255)
         tAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCGCTGGCGCGATaaaga                                                                                                                                       <  1:334714/150‑5 (MQ=255)
         tAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCGCTGGCGCGATaaaga                                                                                                                                       <  1:1764409/150‑5 (MQ=255)
         tAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACttt                                                                                                                                       <  1:1018398/150‑1 (MQ=255)
           aaTTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCGCTGGCGCGATaaagacc                                                                                                                                     <  2:396149/150‑7 (MQ=255)
             ttCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCa                                                                                                                                   <  1:92088/150‑1 (MQ=255)
             ttCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCa                                                                                                                                   >  2:1743959/1‑150 (MQ=255)
                 gTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTa                                                                                                                               >  1:1050884/1‑150 (MQ=255)
                    cAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCt                                                                                                                            >  1:451021/1‑150 (MQ=255)
                     aaGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCtg                                                                                                                           >  1:1974050/1‑150 (MQ=255)
                       gagcCAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGTCGAGGTACTTTTGCACGTATCtgtg                                                                                                                         >  1:1076530/1‑150 (MQ=255)
                       gagaCAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGAT‑‑GAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGcc                                                                                                                       <  1:1338391/150‑1 (MQ=255)
                         gaCAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGcc                                                                                                                       >  1:1241914/1‑150 (MQ=255)
                           cAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCg                                                                                                                     >  1:483980/1‑150 (MQ=255)
                           cAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCg                                                                                                                     >  1:1278402/1‑150 (MQ=255)
                            aaTATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGa                                                                                                                    >  2:1929443/1‑150 (MQ=255)
                            aaTATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGa                                                                                                                    <  1:1873897/150‑1 (MQ=255)
                             atatTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGAt                                                                                                                   <  1:1793777/150‑1 (MQ=255)
                              tatTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATg                                                                                                                  >  1:731425/1‑150 (MQ=255)
                                ttttGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATACa                                                                                                                <  2:336871/150‑1 (MQ=255)
                                    gCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATc                                                                                                            <  1:1674921/150‑1 (MQ=255)
                                     cGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCg                                                                                                           >  1:367407/1‑150 (MQ=255)
                                      gTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGa                                                                                                          >  2:1185620/1‑150 (MQ=255)
                                        aaaGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATa                                                                                                        <  2:1241914/150‑1 (MQ=255)
                                           ggCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACt                                                                                                     <  1:906367/150‑1 (MQ=255)
                                                cccAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAAt                                                                                                >  2:137727/1‑150 (MQ=255)
                                                cccAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCAAGTATCTGTGCCCGATGCAGATCGATAACTGAAAt                                                                                                >  2:803768/1‑150 (MQ=255)
                                                 ccAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAAtt                                                                                               >  2:61340/1‑150 (MQ=255)
                                                     aaTAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTa                                                                                           <  1:747987/150‑1 (MQ=255)
                                                     aaTAAACCTTCAAAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTa                                                                                           <  2:1080311/150‑1 (MQ=255)
                                                       tAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAAt                                                                                         <  1:1481663/150‑1 (MQ=255)
                                                           ccTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCAt                                                                                     >  1:1993365/1‑150 (MQ=255)
                                                             ttCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCt                                                                                   >  1:1134258/1‑150 (MQ=255)
                                                             ttCATAAATCGCCAGAAATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATAGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGACCGATCCAGATCGATAACTGAAAATTTTAATTAATca                                                                                   >  2:920527/1‑149 (MQ=255)
                                                              tCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTAATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCtt                                                                                  >  2:345104/1‑150 (MQ=255)
                                                               cATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTc                                                                                 <  1:973427/150‑1 (MQ=255)
                                                                aTAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTcc                                                                                <  2:1850379/150‑1 (MQ=255)
                                                                 tAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTccc                                                                               >  2:1425389/1‑150 (MQ=255)
                                                                  aaaTCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCg                                                                              >  2:1577287/1‑150 (MQ=255)
                                                                   aaTCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCgg                                                                             >  2:1354548/1‑150 (MQ=255)
                                                                       gCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAg                                                                         <  2:1098677/150‑1 (MQ=255)
                                                                       gCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAg                                                                         >  2:1384116/1‑150 (MQ=255)
                                                                       gCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAg                                                                         >  1:861745/1‑150 (MQ=255)
                                                                        ccAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGt                                                                        <  1:1185620/150‑1 (MQ=255)
                                                                        ccAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGt                                                                        <  2:911017/150‑1 (MQ=255)
                                                                               tctcTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGc                                                                 >  2:625958/1‑150 (MQ=255)
                                                                               tctcTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGc                                                                 >  2:1463971/1‑150 (MQ=255)
                                                                               tctcTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGc                                                                 <  1:1463971/150‑1 (MQ=255)
                                                                               tctcTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGc                                                                 >  1:1397801/1‑150 (MQ=255)
                                                                                ctctTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCa                                                                >  1:442278/1‑150 (MQ=255)
                                                                                  ctTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAGCCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTGCTTTTGCACGTATCTGTGCCCGATGCAGATCGATACCTGAAATTTTTAATTCATCTTCCCGGGCAGTAAATGGCgaa                                                              >  2:203634/1‑150 (MQ=255)
                                                                                   tttGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAAt                                                             >  2:1737518/1‑150 (MQ=255)
                                                                                   tttGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAAt                                                             <  1:135067/150‑1 (MQ=255)
                                                                                      gCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCTCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGaa                                                          <  1:622979/150‑1 (MQ=255)
                                                                                       cGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAg                                                         <  2:861745/150‑1 (MQ=255)
                                                                                          cAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTa                                                                       >  2:1447403/1‑133 (MQ=255)
                                                                                           aGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCaaa                                                     <  1:1929443/150‑1 (MQ=255)
                                                                                            ggCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCAAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCaaaa                                                    >  2:941604/1‑150 (MQ=255)
                                                                                              ccAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAg                                                  >  2:456493/1‑150 (MQ=255)
                                                                                               cAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGt                                                 >  2:543010/1‑150 (MQ=255)
                                                                                               cAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGt                                                 <  1:668836/150‑1 (MQ=255)
                                                                                                  tatTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGc                                              >  1:465201/1‑150 (MQ=255)
                                                                                                  tatTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGc                                              <  1:164887/150‑1 (MQ=255)
                                                                                                  tatTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGc                                              <  2:1616261/150‑1 (MQ=255)
                                                                                                   atTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGGCAGTAAATGGCAAATGAAGCAAAAAGTGGCa                                             >  2:224849/1‑150 (MQ=255)
                                                                                                    ttCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCaa                                            >  2:1814895/1‑150 (MQ=255)
                                                                                                       aTTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATg                                         <  2:569370/150‑1 (MQ=255)
                                                                                                            tgaGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACa                                    >  2:218450/1‑150 (MQ=255)
                                                                                                                  atatAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAAt                                                             <  2:1904316/119‑1 (MQ=255)
                                                                                                                  atatAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAAt                                                             >  1:1904316/1‑119 (MQ=255)
                                                                                                                    ataAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAAt                                                             <  1:1453378/117‑1 (MQ=255)
                                                                                                                    ataAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGt                            <  1:594386/150‑1 (MQ=255)
                                                                                                                     taAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTg                           >  2:1483900/1‑150 (MQ=255)
                                                                                                                       aaCCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGAt                         <  1:1028434/150‑1 (MQ=255)
                                                                                                                       aaCCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGAt                         >  2:1689447/1‑150 (MQ=255)
                                                                                                                          cAATGAGGTCTTTATCGCGCGAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCTCGGTGATTTa                      <  1:282104/150‑1 (MQ=255)
                                                                                                                          cAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTa                      <  1:803768/150‑1 (MQ=255)
                                                                                                                           aaTCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTAt                     >  1:640058/1‑150 (MQ=255)
                                                                                                                               aGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCt                 >  1:1593377/1‑150 (MQ=255)
                                                                                                                                gggatttATCGCGGCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATTCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTTGCAAATGTAACAGTCGCGGTGATTTATCCCTg                <  2:281605/146‑1 (MQ=255)
                                                                                                                                 gTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGc               <  1:1624514/150‑1 (MQ=255)
                                                                                                                                 gTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGc               <  2:1050884/150‑1 (MQ=255)
                                                                                                                                 gTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGc               <  2:1278402/150‑1 (MQ=255)
                                                                                                                                 gTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGc               >  1:1199323/1‑150 (MQ=255)
                                                                                                                                   cTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAg             >  2:1544771/1‑150 (MQ=255)
                                                                                                                                    tttATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGt            <  1:543010/150‑1 (MQ=255)
                                                                                                                                        tCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGc        <  1:1088566/150‑1 (MQ=255)
                                                                                                                                          gcgcCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCAt      >  1:685092/1‑150 (MQ=255)
                                                                                                                                             ccAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCAtttt   >  1:1426953/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGCACGTATCTGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTA  >  NC_000913/3805207‑3805498

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