Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,623 T→A 6.3% E42V (GAA→GTA) ‡ 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,6230TA 6.3% 431.4 / 14.3 143E42V (GAA→GTA) ‡hslUATPase component of the HslVU protease
Reads supporting (aligned to +/- strand):  ref base T (60/74);  new base A (6/3);  total (66/77)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

Only 100 of 146 total aligned reads displayed.

GTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAGCTT  >  NC_000913/4121475‑4121772
                                                                                                                                                    |                                                                                                                                                     
gTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTc                                                                                                                                                      >  1:467519/1‑150 (MQ=255)
gTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTc                                                                                                                                                      <  1:2153130/150‑1 (MQ=255)
 tAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                                                                                                       >  1:1020331/1‑148 (MQ=255)
 tAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                                                                                                       <  2:1020331/148‑1 (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)
    cccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                                                                                                                  <  2:859980/150‑1 (MQ=255)
     ccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGa                                                                                                                                                 <  2:1854657/150‑1 (MQ=255)
       aCTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGc                                                                                                                                               <  2:777968/150‑1 (MQ=255)
        cTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCt                                                                                                                                              <  2:2124854/150‑1 (MQ=255)
         ttCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTg                                                                                                                                             <  1:2013857/150‑1 (MQ=255)
          tCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGc                                                                                                                                            >  1:2130812/1‑150 (MQ=255)
           cGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCa                                                                                                                                           >  1:563137/1‑150 (MQ=255)
           cGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCAACGAAGAGCTGcg                                                                                                                                           <  1:2436316/150‑2 (MQ=255)
             gTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCAACGAAGAGCTGcgcc                                                                                                                                         <  2:1370429/150‑4 (MQ=255)
               gAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGTGCAGCTCAACGAAGAGCTGcgccat                                                                                                                                       <  1:2247882/150‑6 (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)
                      gTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAg                                                                                                                                <  2:387906/150‑1 (MQ=255)
                         gCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGCGCGATTTCAGGTTTACCGACACCGGTCGGGCCGATCATCAGGATATTGTTCGGGGTCACTTCATGGCGCAGGTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                                                                                                             >  2:1334652/1‑150 (MQ=255)
                               aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGc                                                                                                                       <  1:2433515/150‑1 (MQ=255)
                               aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGc                                                                                                                       <  2:575477/150‑1 (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                                                                                                                  <  2:703051/150‑1 (MQ=255)
                                       gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                                                                                                               >  1:73625/1‑150 (MQ=255)
                                         aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                                                                                                                   <  2:2043054/144‑1 (MQ=255)
                                         aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                                                                                                                   >  1:2043054/1‑144 (MQ=255)
                                          cGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGc                                                                                                            >  2:2212045/1‑150 (MQ=255)
                                             cgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcac                                                                                                         <  1:1451087/150‑1 (MQ=255)
                                              gcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcaca                                                                                                        <  1:1005393/150‑1 (MQ=255)
                                               cATTCGCCAGCTTAGCCAGACGACGGGTGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                                                                                                      <  2:1977989/135‑1 (MQ=255)
                                               cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                                                                                                       >  2:1404837/1‑150 (MQ=255)
                                               cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                                                                                                       >  1:2203884/1‑150 (MQ=255)
                                               cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                                                                                                       >  1:2182737/1‑150 (MQ=255)
                                               cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGGCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                                                                                                       <  2:520928/150‑1 (MQ=255)
                                                   cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACg                                                                                                   <  1:642236/150‑1 (MQ=255)
                                                   cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAGTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAAcc                                                                                                   >  2:1139695/1‑149 (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:713983/1‑150 (MQ=255)
                                                                      cgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCg                                                                                >  1:134972/1‑150 (MQ=255)
                                                                          cGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCgatga                                                                            <  1:1958424/150‑1 (MQ=255)
                                                                                   ttttACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTAt                                                                   <  1:958156/150‑1 (MQ=255)
                                                                                    tttACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATc                                                                  >  2:1720702/1‑150 (MQ=255)
                                                                                     ttACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATcc                                                                 <  2:932150/150‑1 (MQ=255)
                                                                                      tACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCa                                                                <  1:1950575/150‑1 (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:1245244/1‑150 (MQ=255)
                                                                                                         gATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgcg                                             <  1:599093/150‑1 (MQ=255)
                                                                                                         gATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgcg                                             >  1:1038536/1‑150 (MQ=255)
                                                                                                          atcatcAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGt                                            <  1:64053/150‑1 (MQ=255)
                                                                                                           tcatcaGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTg                                           <  1:1782466/150‑1 (MQ=255)
                                                                                                              tcaGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTgggg                                        <  1:1612221/150‑1 (MQ=255)
                                                                                                                aGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTc                                      <  2:111202/150‑1 (MQ=255)
                                                                                                                   atatTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTCAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCAtt                                   >  1:240017/1‑150 (MQ=255)
                                                                                                                    tatTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCAttt                                  <  1:87444/150‑1 (MQ=255)
                                                                                                                    tatTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCAttt                                  >  1:596892/1‑150 (MQ=255)
                                                                                                                     atTTTTCGGGGTCACTTCATTGGGCAGCTCTTCGTTGAGCTGCATGCGACGGCAGCGGTTACGCAGAGCAATGGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTc                                 <  2:878327/150‑1 (MQ=255)
                                                                                                                     atTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTc                                 <  2:1368814/150‑1 (MQ=255)
                                                                                                                     atTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTc                                 <  2:897738/150‑1 (MQ=255)
                                                                                                                       ttttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAg                               >  2:1605268/1‑150 (MQ=255)
                                                                                                                       ttttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAg                               <  1:1215539/150‑1 (MQ=255)
                                                                                                                       ttttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAg                               <  1:400447/150‑1 (MQ=255)
                                                                                                                        tttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGa                              <  1:1842786/150‑1 (MQ=255)
                                                                                                                              ggTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATggg                        >  1:87089/1‑150 (MQ=255)
                                                                                                                              ggTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATggg                        >  2:1288304/1‑150 (MQ=255)
                                                                                                                                 cACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCATGGGGTCATTTCAGACATGGgaga                     >  2:570792/1‑150 (MQ=255)
                                                                                                                                  aCTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGAt                    >  1:1002628/1‑150 (MQ=255)
                                                                                                                                     tCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCt                 >  2:2158460/1‑150 (MQ=255)
                                                                                                                                      cATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgcg                                             <  1:722043/121‑1 (MQ=255)
                                                                                                                                      cATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgcg                                             >  2:722043/1‑121 (MQ=255)
                                                                                                                                        tGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTAc              >  2:937497/1‑150 (MQ=255)
                                                                                                                                         catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg             >  1:1142691/4‑150 (MQ=255)
                                                                                                                                         catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg             >  2:242/4‑150 (MQ=255)
                                                                                                                                         catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg             >  2:439525/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)
                                                                                                                                               gCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGt       >  2:1708006/1‑150 (MQ=255)
                                                                                                                                               gCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGt       <  2:2218251/150‑1 (MQ=255)
                                                                                                                                                ctctTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGtt      <  2:1002628/150‑2 (MQ=255)
                                                                                                                                                ctctTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTa      <  2:318542/150‑1 (MQ=255)
                                                                                                                                                 tctTCGTTCAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAg     <  2:240017/150‑1 (MQ=255)
                                                                                                                                                    tCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGTGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAGCtt  >  2:1920078/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAGCTT  >  NC_000913/4121475‑4121772

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: