Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,937,044 Δ1 bp 9.1% coding (1461/1776 nt) fucI → L‑fucose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,937,0440G. 9.1% 290.4 / 22.1 88coding (1461/1776 nt)fucIL‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base G (34/46);  new base . (3/5);  total (37/51)
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.91e-01

CACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGG  >  NC_000913/2936895‑2937193
                                                                                                                                                     |                                                                                                                                                      
cACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGAtgtg                                                                                                                                                        <  2:887327/150‑1 (MQ=255)
cACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCAACATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGAtgtg                                                                                                                                                        >  2:1241318/1‑150 (MQ=255)
cACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGGTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGAtgtg                                                                                                                                                        <  2:782185/150‑1 (MQ=255)
  cGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCa                                                                                                                                                      <  1:241977/150‑1 (MQ=255)
       aCTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACa                                                                                                                                                >  2:358847/1‑150 (MQ=255)
          tCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATcc                                                                                                                                             >  2:438663/1‑150 (MQ=255)
              tGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTcaa                                                                                                                                         >  1:725347/1‑150 (MQ=255)
              tGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTCATGCACATCCttcg                                                                                                                                         <  2:374707/150‑4 (MQ=255)
               ggcggcGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTcaac                                                                                                                                        <  2:1243949/150‑1 (MQ=255)
               ggcggcGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTCATGCACATCCttcgg                                                                                                                                        <  2:716098/150‑5 (MQ=255)
                gcggcgGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTTTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTCATGCACATCCttcggc                                                                                                                                       <  1:480399/150‑6 (MQ=255)
                gcggcgGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTCATGCACATCCttcggc                                                                                                                                       <  1:475062/150‑6 (MQ=255)
                 cggcggTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTcaacaa                                                                                                                                      <  2:394435/150‑1 (MQ=255)
                      gTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCa                                                                                                                                 <  2:880645/150‑1 (MQ=255)
                         aCTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCa                                                                                                                              >  1:1194396/1‑150 (MQ=255)
                          ctctTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCaa                                                                                                                             >  1:701552/1‑150 (MQ=255)
                               cccGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCaa                                                                                                                        >  1:1296236/1‑150 (MQ=255)
                                 cGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAAcc                                                                                                                      <  1:358847/150‑1 (MQ=255)
                                     tCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGc                                                                                                                  <  2:1250807/150‑1 (MQ=255)
                                      ccTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGcc                                                                                                                 >  1:265233/1‑150 (MQ=255)
                                      ccTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGcc                                                                                                                 >  1:308569/1‑150 (MQ=255)
                                      ccTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGcc                                                                                                                 <  2:338412/150‑1 (MQ=255)
                                          aCCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAacc                                                                                                             <  1:1470292/150‑1 (MQ=255)
                                           ccGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTTCTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAacca                                                                                                            >  2:444371/1‑150 (MQ=255)
                                           ccGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAacca                                                                                                            <  1:573762/150‑1 (MQ=255)
                                             gAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAaccacc                                                                                                          <  2:945064/150‑1 (MQ=255)
                                             gAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAaccacc                                                                                                          >  2:138495/1‑150 (MQ=255)
                                              aaGGCGGCGTCCCCTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGAGGAATTGAAGAAGGATGTGCATGACATCATCAACAAACGCACCAACTCAACCTGGCCAAACAACt                                                                                                         >  2:287441/1‑150 (MQ=255)
                                                ggcggcGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTgg                                                                                                       <  1:1012269/150‑1 (MQ=255)
                                                  cggcgTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTcaac                                                                                                                                        <  1:482011/115‑1 (MQ=255)
                                                  cggcgTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGtt                                                                                                     <  1:636524/150‑1 (MQ=255)
                                                  cggcgTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATAGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTcaac                                                                                                                                        >  2:482011/1‑115 (MQ=255)
                                                   ggcggcCCGTTCACCATGACTCGTGGCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGttt                                                                                                    <  2:1351247/150‑1 (MQ=255)
                                                        cccGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAAc                                                                                                                                    >  2:163284/1‑113 (MQ=255)
                                                        cccGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAAc                                                                                                                                    <  1:163284/113‑1 (MQ=255)
                                                          cGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACcgc                                                                                             >  2:695633/1‑150 (MQ=255)
                                                            ttCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCgcgc                                                                                           <  2:943522/150‑1 (MQ=255)
                                                             tCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCCCCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGcc                                                                                          <  1:312556/150‑1 (MQ=255)
                                                              cACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCt                                                                                         <  1:1002339/150‑1 (MQ=255)
                                                               aCCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTc                                                                                        <  2:98205/150‑1 (MQ=255)
                                                               aCCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTc                                                                                        >  2:1019806/1‑150 (MQ=255)
                                                                ccATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCa                                                                                       >  2:919381/1‑150 (MQ=255)
                                                                    gACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCgg                                                                                   <  2:701552/150‑1 (MQ=255)
                                                                    gACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCgg                                                                                   <  2:866315/150‑1 (MQ=255)
                                                                      cTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTa                                                                                 <  1:787823/150‑1 (MQ=255)
                                                                      cTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTa                                                                                 <  2:283124/150‑1 (MQ=255)
                                                                        cGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTaaa                                                                               >  1:1353072/1‑150 (MQ=255)
                                                                        cGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTaaa                                                                               >  2:1116330/1‑150 (MQ=255)
                                                                         gtgtCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAg                                                                              <  1:695633/150‑1 (MQ=255)
                                                                          tgtCAACATCATCAAAGGCCTGGGACCGGTACTGGAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGTATGTGCATGAGATGTTCAATAAACGCATCATCTAAAACTGTCCAAGCACCTGTTTTTCACTGATCCTCTCCTGTAATgg                                                                             <  2:572491/150‑1 (MQ=255)
                                                                           gtCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAggg                                                                            <  2:534851/150‑1 (MQ=255)
                                                                           gtCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAggg                                                                            <  2:344138/150‑1 (MQ=255)
                                                                           gtCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAggg                                                                            <  1:763829/150‑1 (MQ=255)
                                                                            tCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGc                                                                           <  2:1020661/150‑1 (MQ=255)
                                                                             cAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGcc                                                                          <  1:723697/150‑1 (MQ=255)
                                                                              aaCATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGc                                                                           >  2:663610/1‑148 (MQ=255)
                                                                              aaCATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGc                                                                           <  1:663610/148‑1 (MQ=255)
                                                                               aCATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGt                                                                        <  1:1301555/150‑1 (MQ=255)
                                                                                catcatCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACGAACTCAACCTGGCCATCCACCTGGTTTGGACCGCGCCTCACCGGTAAAGGGCCGtt                                                                       <  2:847293/150‑1 (MQ=255)
                                                                                catcatCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGtt                                                                       >  1:1054326/1‑150 (MQ=255)
                                                                                   catcaAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTAc                                                                    >  2:168584/1‑150 (MQ=255)
                                                                                     tcaAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACgg                                                                  >  1:346132/1‑150 (MQ=255)
                                                                                       aaaGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGAt                                                                >  1:624839/1‑150 (MQ=255)
                                                                                          ggCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGAtgtg                                                             >  2:78531/1‑150 (MQ=255)
                                                                                          ggCCTGGGACCGGTACTGAAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTCGCACCGCGCCTCACCGGTAAATGGCCGTTTACGGAtgtg                                                             <  2:588590/150‑1 (MQ=255)
                                                                                                 gACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCgg                                                      <  2:1194396/150‑1 (MQ=255)
                                                                                                  aCCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGt                                                     <  1:168584/150‑1 (MQ=255)
                                                                                                     ggTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAAt                                                  >  2:472687/1‑150 (MQ=255)
                                                                                                       tACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATgg                                                <  2:448929/150‑1 (MQ=255)
                                                                                                       tACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATgg                                                <  1:919381/150‑1 (MQ=255)
                                                                                                        aCTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGc                                               <  1:28599/150‑1 (MQ=255)
                                                                                                           gCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAACCGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGaa                                            >  1:1452325/1‑150 (MQ=255)
                                                                                                                tCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTggg                                       >  2:882051/1‑150 (MQ=255)
                                                                                                                    ggAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGgcgc                                   <  1:43530/150‑1 (MQ=255)
                                                                                                                     gAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCt                                  <  2:624839/150‑1 (MQ=255)
                                                                                                                       aGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTaa                                >  1:1001153/1‑150 (MQ=255)
                                                                                                                        ggCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCAGCGCGCCTGACCGGTAAAGGGCCGTTTTCGGATGTGTACTCGGTAATGGCGAACTTGGGCGCTAAc                               <  2:725347/150‑1 (MQ=255)
                                                                                                                         gCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAAcc                              >  1:209468/1‑150 (MQ=255)
                                                                                                                             gAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATgg                          <  2:308569/150‑1 (MQ=255)
                                                                                                                               gCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACg                                                                   >  2:666362/1‑107 (MQ=255)
                                                                                                                               gCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACg                                                                   <  1:666362/107‑1 (MQ=255)
                                                                                                                               gCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATgggg                        <  1:381090/150‑1 (MQ=255)
                                                                                                                                  tGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTc                     >  1:1269574/1‑150 (MQ=255)
                                                                                                                                      aTTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGAc                 >  2:1340978/1‑150 (MQ=255)
                                                                                                                                          gttgaGGATGTCATGCACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATc             >  1:850185/5‑150 (MQ=255)
                                                                                                                                          ccGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATc             >  1:1372340/1‑150 (MQ=255)
                                                                                                                                           ttgaGGATGTCATGCACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCg            >  1:449661/4‑150 (MQ=255)
                                                                                                                                            tgaGGATGTCATGCACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCgg           >  1:851546/3‑150 (MQ=255)
                                                                                                                                            gAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCgg           <  2:1225145/150‑1 (MQ=255)
                                                                                                                                             gaGGATGTCATGCACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGc          <  2:449661/149‑1 (MQ=255)
                                                                                                                                                gATGTGCATGACATCCTCAACAAACGTACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCAc       >  1:432207/1‑150 (MQ=255)
                                                                                                                                                     gCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTgg  >  1:1215848/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CACGAATACTTCCGTGGCGGCGGTTACTCTTCCCGCTTCCTTACCGAAGGCGGCGTCCCGTTCACCATGACTCGTGTCAACATCATCAAAGGCCTGGGACCGGTACTGCAAATCGCGGAAGGCTGGAGCGTGGAATTGCCGAAGGATGTGCATGACATCCTCAACAAACGCACCAACTCAACCTGGCCAACCACCTGGTTTGCACCGCGCCTCACCGGTAAAGGGCCGTTTACGGATGTGTACTCGGTAATGGCGAACTGGGGCGCTAACCATGGGGTTCTGACCATCGGCCACGTTGG  >  NC_000913/2936895‑2937193

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: