Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,623 T→A 7.5% E42V (GAA→GTA) ‡ 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,6230TA 7.5% 335.6 / 13.0 129E42V (GAA→GTA) ‡hslUmolecular chaperone and ATPase component of HslUV protease
Reads supporting (aligned to +/- strand):  ref base T (61/57);  new base A (6/4);  total (67/62)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 129 total aligned reads displayed.

GCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTA  >  NC_000913/4121478‑4121768
                                                                                                                                                 |                                                                                                                                                 
gCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGtt                                                                                                                                               >  1:688474/1‑150 (MQ=255)
gCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGtt                                                                                                                                               >  2:408590/1‑150 (MQ=255)
 cccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                                                                                                              >  1:1581725/1‑150 (MQ=255)
 cccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGGCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                                                                                                              <  2:59858/150‑1 (MQ=255)
  ccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGa                                                                                                                                             <  2:436117/150‑1 (MQ=255)
     ctgcGGTGAATTTGGGCGCGGCTCCTTTGATCCACGGCGCCTTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGTTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGTTCTTCGTTGAGCt                                                                                                                                          <  1:67996/147‑1 (MQ=255)
      ttCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTg                                                                                                                                         <  1:102399/150‑1 (MQ=255)
        cGGTGAATTTGGTCGCTTAAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGCCACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTAGAGCTGcc                                                                                                                                       >  2:288602/1‑149 (MQ=255)
         ggTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCAt                                                                                                                                      >  1:1376565/1‑150 (MQ=255)
          gggAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCAACGAAGAGCTGcgcc                                                                                                                                     <  1:314865/148‑4 (MQ=255)
          gTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCAACGAAGAGCTGcgcc                                                                                                                                     <  1:1506060/150‑4 (MQ=255)
           tGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCGTGGCGCAGCTCTTCGTTGAGCTGCATGc                                                                                                                                    >  2:323954/1‑150 (MQ=255)
                 tGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGcc                                                                                                                              <  2:113002/150‑1 (MQ=255)
                   gTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAg                                                                                                                            <  1:1172948/150‑1 (MQ=255)
                     cGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCg                                                                                                                          <  2:707483/150‑1 (MQ=255)
                       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCATACGACTGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTAACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCCGCGGt                                                                                                                        >  1:1026089/1‑150 (MQ=255)
                       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                                                                                                        >  1:324913/1‑150 (MQ=255)
                       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                                                                                                        >  1:219790/1‑150 (MQ=255)
                          cAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTAc                                                                                                                     <  1:342537/150‑1 (MQ=255)
                           aaCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACg                                                                                                                    >  2:1151136/1‑150 (MQ=255)
                           aaCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACg                                                                                                                    <  2:1581725/150‑1 (MQ=255)
                           aaCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGAATTCAGTTTTACCGACACCGGTCGGGCCGATCATCCGGATATTTTTCGGGGTCTCTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGGTACg                                                                                                                    >  2:253812/1‑150 (MQ=255)
                            aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGc                                                                                                                   >  1:787747/1‑150 (MQ=255)
                             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                                                                                                  >  1:489401/1‑150 (MQ=255)
                               ttGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCCGCGCTTAGTTGAGCTGCATGCGACGCCAGCGGTTACGCaga                                                                                                                >  1:1033572/1‑150 (MQ=255)
                                tgatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                                                                                                               >  2:737983/1‑150 (MQ=255)
                                 gctgaACGGCGCAGGCGCCAGCGTAGCCAGACGCCGGGCGATTTCAGTTTGACCGACACCGGTCGGGCCGAGCATCAGGATATTTATCGGGGGCACTTCATGGCGCAGCTCTGCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGc                                                                                                              <  1:882608/148‑1 (MQ=255)
                                     aagggCGCATTCGGTAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGAAATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATc                                                                                                          <  2:594412/147‑1 (MQ=255)
                                      aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCg                                                                                                         <  2:688474/150‑1 (MQ=255)
                                      aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCg                                                                                                         <  1:494102/150‑1 (MQ=255)
                                         gcgcATTCGCCAGCTTAGCCAGACGACGGGGGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGAAATTATTCGGGGTCCCTTCATGGCGCCGCTCTTCGGTGAGCTGCATGCGACGCCAGCTGTAACGCAGAGCAATCGccc                                                                                                      >  1:332049/1‑149 (MQ=255)
                                          cgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcac                                                                                                     <  2:432401/150‑1 (MQ=255)
                                          cgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcac                                                                                                     >  2:715932/1‑150 (MQ=255)
                                            cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                                                                                                   <  1:551771/150‑1 (MQ=255)
                                                cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACg                                                                                               >  2:1327753/1‑150 (MQ=255)
                                                 gCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGc                                                                                              <  2:401505/150‑1 (MQ=255)
                                                  ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATAATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCt                                                                                             >  1:1507402/1‑150 (MQ=255)
                                                        tAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGt                                                                                       >  2:1145621/1‑150 (MQ=255)
                                                            cAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTc                                                                                   <  1:290205/150‑1 (MQ=255)
                                                             aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGAAACTTCATGCCCCAGCTCTTCGTTTAGCTGAATGCGACGCCACCGGTTCCGCATACCAAACGCCACAGAACGCTTGCCGTTGAcc                                                                                  >  1:1373213/1‑150 (MQ=255)
                                                             aGACGACGGGCGATTTCAGTGTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGTGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc                                                                                  <  2:671673/150‑1 (MQ=255)
                                                               acgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTg                                                                                <  2:1060627/150‑1 (MQ=255)
                                                               acgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGCCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTg                                                                                <  1:301679/150‑1 (MQ=255)
                                                                cgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATACTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTgg                                                                               >  2:1635907/1‑150 (MQ=255)
                                                                  ccgGGCGATTTCAGTTTTACCGACACCGGGCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGcc                                                                             <  1:1493418/149‑1 (MQ=255)
                                                                   cgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCGCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCg                                                                            <  1:7704/150‑1 (MQ=255)
                                                                     ggCGATTTCAGGTTTACCGACACCGGTCGGGCCAATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCgat                                                                          <  2:273901/150‑1 (MQ=255)
                                                                      gCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCgatg                                                                         >  2:988344/1‑150 (MQ=255)
                                                                      gCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGGTGTCCTGGCCgatg                                                                         >  1:1286679/1‑150 (MQ=255)
                                                                         aTTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGGTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGatgatg                                                                      >  1:72633/1‑150 (MQ=255)
                                                                          tttCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGAtgt                                                                     >  2:859653/1‑150 (MQ=255)
                                                                               tttttACCGAGACCGGTCGGGCCGATCATCAGTATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTa                                                                <  2:1378911/149‑1 (MQ=255)
                                                                                     ccGACACCGGTCGGGCCGATCATCAGTATATTTTTCGGGGTCACTTAATGGCGCAGATCTAAGTTAAGAAGCATGCGACGCCAGCGGATACGCAGAGCAATCGCCACACAACGCTTGGCGTTGTCCTGGCCGATGATCTTCTTATCCAGt                                                          >  2:1364159/1‑150 (MQ=255)
                                                                                     ccGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGt                                                          <  1:996339/150‑1 (MQ=255)
                                                                                        acacCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCg                                                       >  2:641558/1‑150 (MQ=255)
                                                                                          acCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCt                                                     <  2:1633477/150‑1 (MQ=255)
                                                                                            gggTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGa                                                   <  1:279616/149‑1 (MQ=255)
                                                                                             ggTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGAc                                                  <  1:641558/150‑1 (MQ=255)
                                                                                              gTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACg                                                 <  2:1600570/150‑1 (MQ=255)
                                                                                                 gggCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGAtt                                              >  1:1128727/1‑150 (MQ=255)
                                                                                                   gCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGGCGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGCTGTGCTTATCCAGTTCGCTGAAGAttta                                            >  1:1055579/1‑149 (MQ=255)
                                                                                                    ccGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcg                                           >  1:860045/1‑150 (MQ=255)
                                                                                                    ccGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcg                                           <  1:540590/150‑1 (MQ=255)
                                                                                                    ccGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcg                                           >  2:771204/1‑150 (MQ=255)
                                                                                                    ccGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCGTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcg                                           >  2:132327/1‑150 (MQ=255)
                                                                                                     cGATCATCAGGATATTTTTCGGGGTCACTTCATGGTGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTcgc                                          >  2:1043123/1‑150 (MQ=255)
                                                                                                         catcaGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTgg                                      <  2:736114/150‑1 (MQ=255)
                                                                                                              ggATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCa                                 >  2:1293835/1‑150 (MQ=255)
                                                                                                              ggATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTGAGCCAGTTCGCTGCCGATGTCGCGTGGGGTCa                                 >  2:1047190/1‑150 (MQ=255)
                                                                                                                 tatTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCAttt                              >  1:1483557/1‑150 (MQ=255)
                                                                                                                  atTTTTCGGGGTCTCTTCATGGCGCATGTCTTCGTTGAGTTGCATGCGACGCGAGCGGTTACGCAGAGCAATCGCCACAGAACTCTTGGAGTTGTCCTGGCCTATTATTTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTc                             <  2:992346/150‑1 (MQ=255)
                                                                                                                    ttttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGCACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAg                           >  2:827535/1‑150 (MQ=255)
                                                                                                                    ttttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAg                           >  1:444676/1‑150 (MQ=255)
                                                                                                                      ttCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTACGCGTGGGGTCATTTCAGAc                         <  1:1473337/150‑1 (MQ=255)
                                                                                                                        cGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGGTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGCCGATTTCGCGTGGGGTCATTTCAGACAt                       >  1:868912/1‑150 (MQ=255)
                                                                                                                          gggTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATgg                     <  2:324913/150‑1 (MQ=255)
                                                                                                                          gggTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATgg                     >  2:1698112/1‑150 (MQ=255)
                                                                                                                           ggTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATggg                    <  1:1364159/150‑1 (MQ=255)
                                                                                                                             tCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGgaga                 <  1:408590/150‑1 (MQ=255)
                                                                                                                             gcACTTCCTGGCGCAGCGCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCATTCGCCACAGAACGCTTGGCGTTGTCTTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGgag                  <  1:253812/149‑1 (MQ=255)
                                                                                                                                cTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACAATTTCGCGTGGGGTCATTTCAGACATGGGAGATc               >  1:657545/1‑150 (MQ=255)
                                                                                                                                cTTCATGGCGCAGCGCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGTCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATc               <  1:1008266/150‑1 (MQ=255)
                                                                                                                                 ttCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATcc              <  2:423239/150‑1 (MQ=255)
                                                                                                                                 ttCATGGCGCAGCTCTTCGTGGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATcc              <  2:1507402/150‑1 (MQ=255)
                                                                                                                                  tCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCt             >  2:1384266/1‑150 (MQ=255)
                                                                                                                                  tCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCt             >  1:34241/1‑150 (MQ=255)
                                                                                                                                   cATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCtt            >  1:415003/1‑150 (MQ=255)
                                                                                                                                   cATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCtt            >  1:742353/1‑150 (MQ=255)
                                                                                                                                     tGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTAc          >  2:576162/1‑150 (MQ=255)
                                                                                                                                      ggCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATAGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg         <  2:219790/150‑1 (MQ=255)
                                                                                                                                      catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTTGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg         >  2:962768/4‑150 (MQ=255)
                                                                                                                                      catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg         >  1:406477/4‑150 (MQ=255)
                                                                                                                                      catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg         >  2:509196/4‑150 (MQ=255)
                                                                                                                                      catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg         >  1:1679981/4‑150 (MQ=255)
                                                                                                                                      catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg         <  2:1302581/147‑1 (MQ=255)
                                                                                                                                      catgcAGCTCAACGAAGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACg         >  1:1302581/4‑150 (MQ=255)
                                                                                                                                       gcgcAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGc        <  1:1151136/150‑1 (MQ=255)
                                                                                                                                         gcAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCtt      <  1:621802/150‑1 (MQ=255)
                                                                                                                                         gcAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCtt      <  1:1406749/150‑1 (MQ=255)
                                                                                                                                             ctctTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTa  <  1:1327753/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCCTGGCCGATGATGTGCTTATCCAGTTCGCTGACGATTTCGCGTGGGGTCATTTCAGACATGGGAGATCCTTACGCTTTGTA  >  NC_000913/4121478‑4121768

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