Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,627 T→A 5.5% N41Y (AAC→TAC) ‡ hslU ← molecular chaperone and ATPase component of HslUV protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,121,6270TA 5.5% 353.1 / 7.3 124N41Y (AAC→TAC) ‡hslUmolecular chaperone and ATPase component of HslUV protease
Reads supporting (aligned to +/- strand):  ref base T (58/58);  new base A (4/3);  total (62/62)
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.33e-01

Only 100 of 126 total aligned reads displayed.

GCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAGCTTAAT  >  NC_000913/4121478‑4121775
                                                                                                                                                     |                                                                                                                                                    
ccccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGtt                                                                                                                                                      <  1:1162553/149‑1 (MQ=255)
 cccACTTCTGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                                                                                                                     <  2:391254/150‑1 (MQ=255)
 cccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                                                                                                                     >  1:313890/1‑150 (MQ=255)
 cccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                                                                                                                     >  2:2061293/1‑150 (MQ=255)
 cccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                                                                                                                     <  2:1820218/150‑1 (MQ=255)
   cACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAg                                                                                                                                                   <  1:298951/150‑1 (MQ=255)
       tCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCAACGAAGAGCTGc                                                                                                                                               <  1:976477/150‑1 (MQ=255)
        cGGTGAAATAGGTCGCATCAACTTTGATGAACGGCCCATTCGCCAGCTTAGCCAGACGACGCGCGATTTCAGATATACCGACCCCGGTCGGGCCGATCATCAGGATATTATTCGTGGTCACTTAATGGCGCAGCTCTTCGTTGATCTGCa                                                                                                                                              >  2:1411902/1‑150 (MQ=255)
          gTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATg                                                                                                                                            >  2:129004/1‑150 (MQ=255)
          gTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATg                                                                                                                                            <  1:1495980/150‑1 (MQ=255)
            gAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGcc                                                                                                                                          >  2:978654/1‑149 (MQ=255)
            gAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                                                                                                          <  1:12576/150‑1 (MQ=255)
              aTTTGGTCGCTTCAACTTTGATGAACGGCGCATTCTCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGATCTTCGTTGAGCTGCATGCGAc                                                                                                                                        >  2:1935683/1‑150 (MQ=255)
              aTTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGTTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                                                                                                        >  1:127269/1‑150 (MQ=255)
              aTTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGGCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                                                                                                        <  2:1162397/150‑1 (MQ=255)
                   gTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAg                                                                                                                                   >  1:635814/1‑150 (MQ=255)
                   gTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAg                                                                                                                                   >  1:896506/1‑150 (MQ=255)
                    tCGCTTCAACATTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGCCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGc                                                                                                                                  <  1:565300/150‑1 (MQ=255)
                      gCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGTCGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                                                                                                                <  1:2208431/150‑1 (MQ=255)
                      gCTTCAACTTAGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGGCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                                                                                                                <  1:771287/150‑1 (MQ=255)
                           aaCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCCGCTCTTCGTTGAGCTGCCTGCGACGCCCGCGGTTACg                                                                                                                           >  2:391995/1‑150 (MQ=255)
                            aCTTTGATCAACGCCGAATAAGCCAGCTTAGCCAGACCACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGGCACTTCATGGGGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGc                                                                                                                          <  2:1073974/150‑1 (MQ=255)
                             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                                                                                                         >  2:1365529/1‑150 (MQ=255)
                              tttGATGAACGGCGCATGCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                                                                                                        <  1:950453/150‑1 (MQ=255)
                                tgatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                                                                                                                      >  2:1619494/1‑150 (MQ=255)
                                 gatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGc                                                                                                                     >  1:68288/1‑150 (MQ=255)
                                  atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCATTGAGCTGCATGCGACGCCAGCGGTTCCGCAGAGCa                                                                                                                    >  2:1837584/1‑150 (MQ=255)
                                   aagctggCGAATGCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCaa                                                                                                                   >  2:410009/6‑150 (MQ=255)
                                    gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                                                                                                                  >  1:506765/1‑150 (MQ=255)
                                      aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGCGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCg                                                                                                                >  1:1844227/1‑150 (MQ=255)
                                         gcgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCca                                                                                                             >  2:832550/1‑150 (MQ=255)
                                         gcgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCca                                                                                                             <  2:1005777/150‑1 (MQ=255)
                                           gcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcaca                                                                                                           <  1:2183232/150‑1 (MQ=255)
                                            cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                                                                                                          <  1:1890114/150‑1 (MQ=255)
                                               tCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAAc                                                                                                       >  2:1254922/1‑150 (MQ=255)
                                                cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACg                                                                                                      >  1:1507554/1‑150 (MQ=255)
                                                 gCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGc                                                                                                     >  1:1767897/1‑150 (MQ=255)
                                                   cAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCtt                                                                                                   >  2:730700/1‑150 (MQ=255)
                                                       ttAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCg                                                                                               <  1:1207302/150‑1 (MQ=255)
                                                          gCCAGACGACGGGCGATTTCAGTTTTACCGACTCAGGTCGGGCCGAACATCAGGATATATTTCGGGGTCACTTCATGGCGCAGCTCTTCGATGAGCTCCATGCGAAGCCAGCGGATACGCAGAGCAATCGCCACAGATCGCTTGGCGTTg                                                                                            >  1:375424/1‑150 (MQ=255)
                                                          gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg                                                                                            <  2:1165124/150‑1 (MQ=255)
                                                          gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCGCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGttt                                                                                            >  2:589692/1‑149 (MQ=255)
                                                          cccAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg                                                                                            <  2:1845174/149‑1 (MQ=255)
                                                            cAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTc                                                                                          >  1:891049/1‑150 (MQ=255)
                                                                 gacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGc                                                                                     >  2:1620874/1‑150 (MQ=255)
                                                                 gacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGc                                                                                     <  1:1872321/150‑1 (MQ=255)
                                                                  acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGcc                                                                                    <  1:827596/150‑1 (MQ=255)
                                                                  acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGcc                                                                                    <  2:732052/150‑1 (MQ=255)
                                                                    gggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCGCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCga                                                                                  <  2:974630/150‑1 (MQ=255)
                                                                     ggCGATTTCAGTTCTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGAGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCgat                                                                                 >  1:465273/1‑150 (MQ=255)
                                                                      gCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCgatg                                                                                <  2:1459631/150‑1 (MQ=255)
                                                                      gCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCgatg                                                                                <  2:1587942/150‑1 (MQ=255)
                                                                           ttCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGAtgtg                                                                           >  1:1463435/1‑150 (MQ=255)
                                                                               gTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTa                                                                       >  2:711017/1‑150 (MQ=255)
                                                                                 tttCCCGACACCGGGCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCGCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATc                                                                     <  1:1335222/150‑1 (MQ=255)
                                                                                   tACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCa                                                                   >  1:2189823/1‑150 (MQ=255)
                                                                                   tACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCa                                                                   <  2:1907242/150‑1 (MQ=255)
                                                                                    aCCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAg                                                                  >  1:836095/1‑150 (MQ=255)
                                                                                     ccGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGt                                                                 >  1:1307407/1‑150 (MQ=255)
                                                                                       gACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTc                                                               <  1:1246002/150‑1 (MQ=255)
                                                                                          acCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCCCTTCATGGCGCAGCTCTTCGTTGCGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGGTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCt                                                            >  2:1077385/1‑150 (MQ=255)
                                                                                          acCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCt                                                            <  2:595479/150‑1 (MQ=255)
                                                                                           ccGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTATAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTg                                                           <  2:639228/150‑1 (MQ=255)
                                                                                              gTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTGCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACg                                                        <  1:1541325/150‑1 (MQ=255)
                                                                                                  ggCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGTAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGAttt                                                    >  2:1391581/1‑150 (MQ=255)
                                                                                                  ggCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGAttt                                                    >  1:382869/1‑150 (MQ=255)
                                                                                                     cGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgc                                                 >  1:564492/1‑150 (MQ=255)
                                                                                                         catcaGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTgg                                             >  1:1688034/1‑150 (MQ=255)
                                                                                                         catcaGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTgg                                             >  1:324874/1‑150 (MQ=255)
                                                                                                          atcaGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTggg                                            >  2:1974900/1‑150 (MQ=255)
                                                                                                           tcaGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTgggg                                           <  2:583521/150‑1 (MQ=255)
                                                                                                                  atTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGATGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCAGACGATTTCGCCTGGGGTCATTTc                                    >  2:536980/1‑150 (MQ=255)
                                                                                                                  aattttCGGGGTCACATCAAGGCGCAGCTCTTCGTTTAGCTGCAGGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTc                                    <  2:381858/148‑1 (MQ=255)
                                                                                                                    ttttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAg                                  <  1:70415/150‑1 (MQ=255)
                                                                                                                    ttttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGAGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAg                                  >  1:888144/1‑150 (MQ=255)
                                                                                                                     tttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAGCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGa                                 <  1:1770818/150‑1 (MQ=255)
                                                                                                                      ttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGTCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGAc                                <  1:978654/150‑1 (MQ=255)
                                                                                                                       tCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACa                               <  1:871183/150‑1 (MQ=255)
                                                                                                                        cGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACAt                              >  2:957669/1‑150 (MQ=255)
                                                                                                                         ggggTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATg                             >  2:13799/1‑150 (MQ=255)
                                                                                                                         ggggTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATg                             >  1:82690/1‑150 (MQ=255)
                                                                                                                          gggTCACTTCATGGCGCAGCTCTTCGTTGCGCTGCATGCGACGCCAGCGGTAACGCAGAGCAATCGCCACAGAACGCTTGGCGTTCTCCTGGCCGATGATGTGCTGATCCAGTTCGCTGACGAATACGCGTGGGGTCAGTTCAGACATgg                            >  1:1060621/1‑150 (MQ=255)
                                                                                                                            ggcACTTCATGTCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCGACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGTTTTCGCGTGGGTTCATTTCAGACATGGga                          <  1:1411902/148‑1 (MQ=255)
                                                                                                                             tCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCCGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGgag                         <  1:410009/150‑1 (MQ=255)
                                                                                                                             tCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGgag                         <  1:1912033/150‑1 (MQ=255)
                                                                                                                             gcACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGgag                         <  2:1466927/149‑1 (MQ=255)
                                                                                                                                 ttCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATcc                     <  2:506765/150‑1 (MQ=255)
                                                                                                                                  tCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCt                    <  1:1432755/150‑1 (MQ=255)
                                                                                                                                  tCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCt                    >  2:526754/1‑150 (MQ=255)
                                                                                                                                   tcgcatgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGCCGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCtt                   >  1:662957/7‑150 (MQ=255)
                                                                                                                                   tcgcatgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCtt                   >  2:497992/7‑150 (MQ=255)
                                                                                                                                     gcatgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTAc                 >  1:1757477/5‑150 (MQ=255)
                                                                                                                                      catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg                >  1:2069537/4‑150 (MQ=255)
                                                                                                                                        cgcAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCt              >  2:108325/1‑150 (MQ=255)
                                                                                                                                            gCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGt          <  2:1463435/150‑1 (MQ=255)
                                                                                                                                               ctTCGTGGAGCTGCATGCGACGCCAGCGGGTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAGc       <  1:1139520/150‑1 (MQ=255)
                                                                                                                                                ttCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAGCt      <  2:1507554/150‑1 (MQ=255)
                                                                                                                                                 tCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAACCGCCACAGCACGCTTTGCGTTGTCATGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAGCtt     <  1:1101085/150‑1 (MQ=255)
                                                                                                                                                   gTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAGCTTaa   >  2:1616095/1‑150 (MQ=255)
                                                                                                                                                    ttGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAGCTTAAt  <  1:279130/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTAGCTTAAT  >  NC_000913/4121478‑4121775

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: