Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 494,032 A→T 8.2% intergenic (+9/‑44) dnaX → / → ybaB DNA polymerase III subunit tau/putative nucleoid‑associated protein YbaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913494,0320AT 8.2% 490.7 / 23.5 168intergenic (+9/‑44)dnaX/ybaBDNA polymerase III subunit tau/putative nucleoid‑associated protein YbaB
Reads supporting (aligned to +/- strand):  ref base A (69/85);  new base T (3/11);  total (72/96)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

Only 100 of 172 total aligned reads displayed.

ACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCG  >  NC_000913/493883‑494181
                                                                                                                                                     |                                                                                                                                                     
   ccgcTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAg                                                                                                                                                    >  2:774366/1‑150 (MQ=255)
   ccgcTGGAATGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAg                                                                                                                                                    >  2:1213423/1‑150 (MQ=255)
         gAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACg                                                                                                                                              <  1:1970120/150‑1 (MQ=255)
          aGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGt                                                                                                                                             >  1:892207/1‑150 (MQ=255)
           gTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGAACGTAAGCTGTGCTTACgat                                                                                                                                            <  1:74271/150‑3 (MQ=255)
            tGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTc                                                                                                                                           >  2:1885583/1‑150 (MQ=255)
             ggCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCTCCATTTGATCGTAAGCTGTGCTTACGATca                                                                                                                                          <  2:1321945/150‑2 (MQ=255)
             ggCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCg                                                                                                                                          <  1:1971684/150‑1 (MQ=255)
             ggCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCg                                                                                                                                          >  1:1121486/1‑150 (MQ=255)
             ggCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGAACGTAAGCTGTGCTTACGATca                                                                                                                                          <  2:1262033/150‑2 (MQ=255)
              gCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGt                                                                                                                                         >  2:866101/1‑150 (MQ=255)
              gCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGt                                                                                                                                         <  1:573070/150‑1 (MQ=255)
              gCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGt                                                                                                                                         <  1:954992/150‑1 (MQ=255)
               cGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTc                                                                                                                                        <  1:600462/150‑1 (MQ=255)
                 tCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGAACGTAAGCTGTGCTTACGATcaaatg                                                                                                                                      <  1:519330/150‑6 (MQ=255)
                  cAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCTGTGCTTACGATcaaatgg                                                                                                                                     <  1:146818/150‑7 (MQ=255)
                  cAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCTGTGCTTACGATcaaatgg                                                                                                                                     <  1:452895/150‑7 (MQ=255)
                  cAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGAACGTAAGCTGTGCTTACGATcaaatgg                                                                                                                                     <  2:2077647/150‑7 (MQ=255)
                   aGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATcc                                                                                                                                    >  1:1648793/1‑150 (MQ=255)
                   aGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGAACGTAAGCTGTGCTTACGATcaaatggg                                                                                                                                    <  1:638491/150‑8 (MQ=255)
                    ggCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCt                                                                                                                                   >  2:2026978/1‑150 (MQ=255)
                    ggCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCt                                                                                                                                   >  1:648205/1‑150 (MQ=255)
                    ggCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCt                                                                                                                                   <  2:241865/150‑1 (MQ=255)
                       gATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTaa                                                                                                                                >  1:2003811/1‑150 (MQ=255)
                        atatACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAAc                                                                                                                               <  2:1814339/150‑1 (MQ=255)
                          ataCGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGATTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGt                                                                                                                             >  1:215333/1‑150 (MQ=255)
                           taCGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCa                                                                                                                                       <  2:1927116/139‑1 (MQ=255)
                            aCGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGa                                                                                                                           >  2:1039238/1‑150 (MQ=255)
                            aCGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGa                                                                                                                           >  2:2341430/1‑150 (MQ=255)
                             cGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGAt                                                                                                                          >  1:562333/1‑150 (MQ=255)
                               aagaaAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTg                                                                                                                        >  1:1387351/1‑150 (MQ=255)
                               aagaaAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTg                                                                                                                        >  1:1458774/1‑150 (MQ=255)
                                agaaAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTga                                                                                                                       >  1:1095271/1‑150 (MQ=255)
                                 gaaAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTgag                                                                                                                      <  1:2392363/150‑1 (MQ=255)
                                   aaaaCTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTgagag                                                                                                                    >  2:1884975/1‑150 (MQ=255)
                                     aaCTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTgagagag                                                                                                                  >  2:1308440/1‑150 (MQ=255)
                                         tGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAAc                                                                                                              <  2:274953/150‑1 (MQ=255)
                                          gcgcAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAAcc                                                                                                             >  1:1403025/1‑150 (MQ=255)
                                           cgcAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCt                                                                                                            <  2:517761/150‑1 (MQ=255)
                                                         tCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTCCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAgg                                                                                              <  1:361368/150‑1 (MQ=255)
                                                         tCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCAATTGATCGTAAGCACAGCTTACGTTCGTAATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAgg                                                                                              >  2:2479641/1‑150 (MQ=255)
                                                           cATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCg                                                                                            >  2:1495722/1‑150 (MQ=255)
                                                            attattGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTATGCACAGCTTACGTTCGTCATCCTTATCGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCgg                                                                                           <  2:1637799/150‑1 (MQ=255)
                                                                 tGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTgg                                                                                      <  1:1098462/150‑1 (MQ=255)
                                                                 tGCGGATAATAATATTCAGACCCTGCGTCGGGTCTTCGATGCGGAGCTGGAGGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTgg                                                                                      <  2:414218/150‑1 (MQ=255)
                                                                  gCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGTAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTggg                                                                                     <  2:1682803/150‑1 (MQ=255)
                                                                  gCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTCCGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTggg                                                                                     <  1:1056127/150‑1 (MQ=255)
                                                                   cGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGt                                                                                    <  1:1800709/150‑1 (MQ=255)
                                                                    ggATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTa                                                                                   <  2:607045/150‑1 (MQ=255)
                                                                       taataataTTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAAcc                                                                                >  2:1792818/1‑150 (MQ=255)
                                                                       taataataTTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCAGGATGAAGAAAGTATCCGCCCCATTTGATCCTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAAcc                                                                                >  2:2265429/1‑150 (MQ=255)
                                                                        aataattTTGAGTGCGTGGGTGGGTTCTTCGATGGGGAGCTGGTTTTAGAATGTTTCCGCCCCATTTGCTCGTAAGCACAGCTTACGTGAGTCATCGTTAACGTGATTGAGAGAGAAACCTTTGTTTGGTAAAGGCGGTCTGGGTAACCt                                                                               <  2:1733339/150‑1 (MQ=255)
                                                                         ataataTTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCtg                                                                              >  2:83744/1‑150 (MQ=255)
                                                                          taataTTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTTTCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAAcc                                                                                                             <  1:1790177/118‑1 (MQ=255)
                                                                          taataTTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTTTCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAAcc                                                                                                             >  2:1790177/1‑118 (MQ=255)
                                                                            atatTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCtgatg                                                                           <  1:1007287/150‑1 (MQ=255)
                                                                             tatTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCtgatga                                                                          <  1:1138867/150‑1 (MQ=255)
                                                                               ttCAGACGCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTTATGaag                                                                        <  1:1213423/150‑1 (MQ=255)
                                                                               ttCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGaag                                                                        >  2:1533335/1‑150 (MQ=255)
                                                                                   gACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGaagcaag                                                                    >  2:2181597/1‑150 (MQ=255)
                                                                                   gACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGaagcaag                                                                    >  1:1696079/1‑150 (MQ=255)
                                                                                      ccTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGccc                                                                 <  1:894953/150‑1 (MQ=255)
                                                                                           gTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCcagcag                                                            <  2:621214/150‑1 (MQ=255)
                                                                                           gTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCcagcag                                                            <  1:774366/150‑1 (MQ=255)
                                                                                             cGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGAt                                                          <  2:714742/150‑1 (MQ=255)
                                                                                              ggTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTTGGTAACCTGATGAAGCAAGCCCAGCAGATg                                                         <  1:269844/150‑1 (MQ=255)
                                                                                              ggTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATg                                                         >  2:2359463/1‑150 (MQ=255)
                                                                                                ttcttcGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCTTTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCa                                                       <  1:2290330/150‑1 (MQ=255)
                                                                                                      gATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGagaa                                                 <  2:138180/150‑1 (MQ=255)
                                                                                                      gATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGaaaa                                                 <  2:1864569/150‑1 (MQ=255)
                                                                                                         gCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaat                                              <  2:1998977/150‑1 (MQ=255)
                                                                                                          cGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatg                                             >  2:105351/1‑150 (MQ=255)
                                                                                                          cGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatg                                             <  1:987861/150‑1 (MQ=255)
                                                                                                              gCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcaga                                         >  2:642593/1‑150 (MQ=255)
                                                                                                              gCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcaga                                         <  1:2353089/150‑1 (MQ=255)
                                                                                                               cTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaa                                        <  2:155673/150‑1 (MQ=255)
                                                                                                               cTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaa                                        <  2:1176001/150‑1 (MQ=255)
                                                                                                                     gaagaaAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCTTTCTGGGTAACCTGATGATGGAAGCCCAGCAGATTCAAGAAaaaatgcagaaaattca                                  <  2:573773/150‑1 (MQ=255)
                                                                                                                     gaagaaAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaaaatgca                                  <  1:2178755/150‑1 (MQ=255)
                                                                                                                       agaaAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAgg                                >  2:909791/1‑150 (MQ=255)
                                                                                                                       agaaAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAgg                                >  1:1616990/1‑150 (MQ=255)
                                                                                                                       agaaAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAgg                                >  1:2414819/1‑150 (MQ=255)
                                                                                                                                 cGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgc                                            <  2:1586341/128‑1 (MQ=255)
                                                                                                                                 cGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgc                                            >  1:1586341/1‑128 (MQ=255)
                                                                                                                                 cGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATcgc                      >  1:1261702/1‑150 (MQ=255)
                                                                                                                                 cGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATcgc                      >  1:1814364/1‑150 (MQ=255)
                                                                                                                                  gCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATcgcg                     >  2:2229006/1‑150 (MQ=255)
                                                                                                                                   ccccATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCgcgc                    <  1:2193845/150‑1 (MQ=255)
                                                                                                                                     ccATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAg                  >  1:1038953/1‑150 (MQ=255)
                                                                                                                                     ccATTTGATCGTAAGCAAAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAg                  <  1:2309528/150‑1 (MQ=255)
                                                                                                                                            aacGTAAGCTGTGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAg           >  2:823193/3‑150 (MQ=255)
                                                                                                                                            aTCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAg           <  2:1121486/150‑1 (MQ=255)
                                                                                                                                            aTCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAACGCAGGAAGAGATCGCGCAGCTGGAAg           <  2:648205/150‑1 (MQ=255)
                                                                                                                                            aTCGTAAGCACAGCTTACGTTCGTCATACTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGCAAGAGATCGCGCAGCTGGAAg           >  1:721678/1‑150 (MQ=255)
                                                                                                                                              cGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTc         <  1:256802/150‑1 (MQ=255)
                                                                                                                                                tAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCAc       <  2:1223983/150‑1 (MQ=255)
                                                                                                                                                  aGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCg     <  1:1885583/150‑1 (MQ=255)
                                                                                                                                                    cacaGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGc   <  1:591598/150‑1 (MQ=255)
                                                                                                                                                    cacaGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGc   <  2:411600/150‑1 (MQ=255)
                                                                                                                                                     acaGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCg  <  2:562333/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCG  >  NC_000913/493883‑494181

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: