Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,622 T→A 6.4% E42D (GAA→GAT) ‡ hslU ← ATPase component of the HslVU protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,121,6220TA 6.4% 437.5 / 14.4 140E42D (GAA→GAT) ‡hslUATPase component of the HslVU protease
Reads supporting (aligned to +/- strand):  ref base T (60/71);  new base A (6/3);  total (66/74)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

Only 100 of 145 total aligned reads displayed.

GTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAG  >  NC_000913/4121475‑4121769
                                                                                                                                                   |                                                                                                                                                   
gTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTc                                                                                                                                                   >  1:467519/1‑150 (MQ=255)
gTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTc                                                                                                                                                   <  1:2153130/150‑1 (MQ=255)
 tAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                                                                                                    >  1:1020331/1‑148 (MQ=255)
 tAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCg                                                                                                                                                  <  2:627663/150‑1 (MQ=255)
 tAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCg                                                                                                                                                  <  2:1271908/150‑1 (MQ=255)
 tAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCg                                                                                                                                                  >  2:661385/1‑150 (MQ=255)
  aGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGt                                                                                                                                                 >  1:1614365/1‑150 (MQ=255)
   gCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGtt                                                                                                                                                <  1:711300/150‑1 (MQ=255)
    cccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                                                                                                               <  2:859980/150‑1 (MQ=255)
    cccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGAAGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                                                                                                               >  2:2340485/1‑150 (MQ=255)
       aCTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGc                                                                                                                                            <  2:777968/150‑1 (MQ=255)
       aCTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGc                                                                                                                                            <  1:2487826/150‑1 (MQ=255)
        cTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCt                                                                                                                                           <  2:2124854/150‑1 (MQ=255)
          tCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGc                                                                                                                                         >  1:2130812/1‑150 (MQ=255)
           cGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCAACGAAGAGCTGcg                                                                                                                                        <  1:2436316/150‑2 (MQ=255)
             gTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCAACGAAGAGCTGcgcc                                                                                                                                      <  2:1370429/150‑4 (MQ=255)
               gAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                                                                                                    >  1:166547/1‑150 (MQ=255)
                 aTTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                                                                                                  <  1:1511126/150‑1 (MQ=255)
                  tttGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACg                                                                                                                                 <  2:10143/150‑1 (MQ=255)
                     ggTCGCTTCAACTTTTATGAACGGCGCATTCGCCAGCTTAGCCAGACGACTGGCGATTTAAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGGTGCATGCGACGCCa                                                                                                                              <  2:511960/150‑1 (MQ=255)
                        cGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTTAGCTGCATGCGACGCCAGCg                                                                                                                           <  2:1526619/150‑1 (MQ=255)
                         gCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGCGCGATTTCAGGTTTACCGACACCGGTCGGGCCGATCATCAGGATATTGTTCGGGGTCACTTCATGGCGCAGGTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                                                                                                          >  2:1334652/1‑150 (MQ=255)
                            tCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTa                                                                                                                       <  1:143076/150‑1 (MQ=255)
                             cAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTAc                                                                                                                      <  2:101212/150‑1 (MQ=255)
                               aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGc                                                                                                                    <  2:575477/150‑1 (MQ=255)
                                cTTTGATGAACGGCGCATTCGCCAGCTTTGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                                                                                                   >  2:1640318/1‑150 (MQ=255)
                                  ttGATGAACGGCGCATTCGCCAGCTTAGCCAGCCGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCaga                                                                                                                 >  1:1171665/1‑150 (MQ=255)
                                   tgatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                                                                                                                >  1:111202/1‑150 (MQ=255)
                                   tgatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                                                                                                                <  2:457200/150‑1 (MQ=255)
                                    gatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGTGGTTACGCAGAGc                                                                                                               <  2:910566/150‑1 (MQ=255)
                                    gatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGc                                                                                                               >  1:1906651/1‑150 (MQ=255)
                                        aaCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATc                                                                                                           >  2:406459/1‑150 (MQ=255)
                                         aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                                                                                                                >  1:2043054/1‑144 (MQ=255)
                                         aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                                                                                                                <  2:2043054/144‑1 (MQ=255)
                                          cGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGc                                                                                                         >  2:2212045/1‑150 (MQ=255)
                                           ggCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                                                                                                        >  2:1080857/1‑150 (MQ=255)
                                              gcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcaca                                                                                                     <  1:1005393/150‑1 (MQ=255)
                                               cATTCGCCAGCTTAGCCAGACGACGGGTGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                                                                                                   <  2:1977989/135‑1 (MQ=255)
                                               cATTCGCCAGCTTAGCCAGACGACGGGTGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                                                                                                   >  1:1977989/1‑135 (MQ=255)
                                               cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                                                                                                    >  1:2182737/1‑150 (MQ=255)
                                               cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                                                                                                    >  1:2203884/1‑150 (MQ=255)
                                               cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGGCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                                                                                                    <  2:520928/150‑1 (MQ=255)
                                                   cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACg                                                                                                <  1:642236/150‑1 (MQ=255)
                                                     ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCt                                                                                              >  2:1842786/1‑150 (MQ=255)
                                                          ttAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCg                                                                                         <  2:1121867/150‑1 (MQ=255)
                                                            aGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGtt                                                                                       >  1:725698/1‑150 (MQ=255)
                                                             gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg                                                                                      >  2:1035195/1‑150 (MQ=255)
                                                                   cgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTgg                                                                                <  1:1378785/150‑1 (MQ=255)
                                                                     acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGcc                                                                              <  2:1117952/150‑1 (MQ=255)
                                                                     acgGGCGATTTCAGTTTTACAGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTAACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCATCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTACTGccc                                                                              >  2:2181765/1‑150 (MQ=255)
                                                                      cgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCg                                                                             >  1:134972/1‑150 (MQ=255)
                                                                        ggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCgat                                                                           <  1:952309/150‑1 (MQ=255)
                                                                          cGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCgatga                                                                         <  1:1958424/150‑1 (MQ=255)
                                                                                   ttttACCGACCCCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTAt                                                                >  2:1348489/1‑150 (MQ=255)
                                                                                    tttACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATc                                                               >  2:1720702/1‑150 (MQ=255)
                                                                                     ttACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATcc                                                              <  2:932150/150‑1 (MQ=255)
                                                                                        ccGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGt                                                           >  1:860988/1‑150 (MQ=255)
                                                                                        ccGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGt                                                           >  2:1446224/1‑150 (MQ=255)
                                                                                         cGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGtt                                                          <  1:2035067/150‑1 (MQ=255)
                                                                                          gACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTc                                                         >  1:1458779/1‑150 (MQ=255)
                                                                                            cacCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGc                                                       >  1:1584049/1‑150 (MQ=255)
                                                                                             acCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCt                                                      <  2:1402989/150‑1 (MQ=255)
                                                                                                ggTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGAc                                                   <  1:1544196/150‑1 (MQ=255)
                                                                                                ggTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGAc                                                   <  2:1614365/150‑1 (MQ=255)
                                                                                                       ccGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcg                                            <  2:725698/150‑1 (MQ=255)
                                                                                                         gATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgcg                                          >  1:1038536/1‑150 (MQ=255)
                                                                                                         gATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgcg                                          <  2:713983/150‑1 (MQ=255)
                                                                                                         gATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgcg                                          >  1:1245244/1‑150 (MQ=255)
                                                                                                         gATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgcg                                          <  1:599093/150‑1 (MQ=255)
                                                                                                          atcatcAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGt                                         <  1:64053/150‑1 (MQ=255)
                                                                                                           tcatcaGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTg                                        <  1:1782466/150‑1 (MQ=255)
                                                                                                              tcaGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTgggg                                     <  1:1612221/150‑1 (MQ=255)
                                                                                                                 ggATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCa                                  >  1:2208019/1‑150 (MQ=255)
                                                                                                                   atatTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTCAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCAtt                                >  1:240017/1‑150 (MQ=255)
                                                                                                                    tatTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCAttt                               >  1:596892/1‑150 (MQ=255)
                                                                                                                    tatTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCAttt                               <  1:87444/150‑1 (MQ=255)
                                                                                                                     atTTTTCGGGGTCACTTCATTGGGCAGCTCTTCGTTGAGCTGCATGCGACGGCAGCGGTTACGCAGAGCAATGGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTc                              <  2:878327/150‑1 (MQ=255)
                                                                                                                     atTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTc                              <  2:897738/150‑1 (MQ=255)
                                                                                                                       ttttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAg                            <  1:1215539/150‑1 (MQ=255)
                                                                                                                        tttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGa                           <  1:1842786/150‑1 (MQ=255)
                                                                                                                              ggTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATggg                     <  1:2409575/150‑1 (MQ=255)
                                                                                                                                 cACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCATGGGGTCATTTCAGACATGGgaga                  >  2:570792/1‑150 (MQ=255)
                                                                                                                                  aCTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGAt                 >  1:468421/1‑150 (MQ=255)
                                                                                                                                  aCTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGAt                 >  1:1002628/1‑150 (MQ=255)
                                                                                                                                  aCTTCATGCCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGAt                 >  2:1966974/1‑150 (MQ=255)
                                                                                                                                     tCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCt              >  2:2158460/1‑150 (MQ=255)
                                                                                                                                     tCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCt              <  2:1402355/150‑1 (MQ=255)
                                                                                                                                      cATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgcg                                          <  1:722043/121‑1 (MQ=255)
                                                                                                                                      cATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgcg                                          >  2:722043/1‑121 (MQ=255)
                                                                                                                                       aTGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTa            >  2:1282816/1‑150 (MQ=255)
                                                                                                                                        tGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTAc           >  2:937497/1‑150 (MQ=255)
                                                                                                                                         catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg          >  1:1142691/4‑150 (MQ=255)
                                                                                                                                         catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg          >  2:439525/4‑150 (MQ=255)
                                                                                                                                         catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg          >  2:303561/4‑150 (MQ=255)
                                                                                                                                         catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTCGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTAGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg          >  2:746683/4‑150 (MQ=255)
                                                                                                                                          gcgcAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGc         >  1:748371/1‑150 (MQ=255)
                                                                                                                                          atgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGc         >  1:174180/3‑150 (MQ=255)
                                                                                                                                                ctctTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGt  <  2:1002628/150‑2 (MQ=255)
                                                                                                                                                ctctTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTa   <  2:318542/150‑1 (MQ=255)
                                                                                                                                                 tctTCGTTCAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAg  <  2:240017/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAG  >  NC_000913/4121475‑4121769

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: