Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 494,034 A→T 10.9% intergenic (+11/‑42) 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,0340AT10.9% 225.5 / 12.4 79intergenic (+11/‑42)dnaX/ybaBDNA polymerase III subunit tau/putative nucleoid‑associated protein YbaB
Reads supporting (aligned to +/- strand):  ref base A (31/39);  new base T (2/7);  total (33/46)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCG  >  NC_000913/493888‑494181
                                                                                                                                                  |                                                                                                                                                   
gcTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCt                                                                                                                                                  >  1:873725/1‑150 (MQ=255)
gcTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCt                                                                                                                                                  >  2:609207/1‑150 (MQ=255)
 cTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCtt                                                                                                                                                 <  2:412158/150‑1 (MQ=255)
   ggAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTAc                                                                                                                                               >  1:431175/1‑150 (MQ=255)
      gTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCTGTGCTTACgat                                                                                                                                            <  2:1064770/150‑3 (MQ=255)
       tGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGTGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCTGTGCTTACGATc                                                                                                                                           <  1:261812/150‑1 (MQ=255)
       tGGCGTCAGGCGATAAACGAAGAAAAACTTCCGCAGGCGCGCGAGTCCATTATTGCGGATATAAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAATAAAGTATCCGCCCCATTTGATCGTAAGCTGTGCTTACGATc                                                                                                                                           <  1:16855/150‑1 (MQ=255)
          cGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCTGTGCTTACGATcaaa                                                                                                                                        <  1:337200/150‑4 (MQ=255)
          cGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTc                                                                                                                                        >  2:304852/1‑150 (MQ=255)
            tCAGGCGATATACGAAGAAAAACTTGCGCAGGGGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGACGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCAt                                                                                                                                      <  2:360794/150‑1 (MQ=255)
            tCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCAt                                                                                                                                      >  1:804551/1‑150 (MQ=255)
              aGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCTGTGCTTACGATcaaatggg                                                                                                                                    <  2:428116/150‑8 (MQ=255)
              aGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATcc                                                                                                                                    <  1:1001733/150‑1 (MQ=255)
               cgCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCACATTTGATCGTAAGCACAGCTTACGTTCGTCATCCt                                                                                                                                   >  1:84481/2‑150 (MQ=255)
                    tataCGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACg                                                                                                                              >  2:206515/1‑150 (MQ=255)
                    tataCGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACg                                                                                                                              <  1:609207/150‑1 (MQ=255)
                       aCGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGa                                                                                                                           <  2:940854/150‑1 (MQ=255)
                        cGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGAt                                                                                                                          >  1:83835/1‑150 (MQ=255)
                        cGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGAt                                                                                                                          >  2:911153/1‑150 (MQ=255)
                          aagaaAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTg                                                                                                                        <  2:236308/150‑1 (MQ=255)
                              aaaaCTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTgagag                                                                                                                    >  1:124373/1‑150 (MQ=255)
                              aaaaCTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTgagag                                                                                                                    <  2:873725/150‑1 (MQ=255)
                                  cTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGaa                                                                                                                <  1:798871/150‑1 (MQ=255)
                                      cgcAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCt                                                                                                            >  2:1148855/1‑150 (MQ=255)
                                       gcAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTa                                                                                                           <  2:1420537/150‑1 (MQ=255)
                                       gcAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTa                                                                                                           >  2:667136/1‑150 (MQ=255)
                                             gcgcgAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTg                                                                                                     <  2:580240/150‑1 (MQ=255)
                                              cgcgAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTgg                                                                                                    <  2:804551/150‑1 (MQ=255)
                                              cgcgAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTgg                                                                                                    >  1:333474/1‑150 (MQ=255)
                                               gcgAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGt                                                                                                   <  2:1162399/150‑1 (MQ=255)
                                               gcgAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGt                                                                                                   <  1:157151/150‑1 (MQ=255)
                                                cgAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTa                                                                                                  <  2:102035/150‑1 (MQ=255)
                                                    tCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAgg                                                                                              <  1:642750/150‑1 (MQ=255)
                                                      cATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCg                                                                                            <  2:83835/150‑1 (MQ=255)
                                                      cATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCg                                                                                            >  1:233132/1‑150 (MQ=255)
                                                             gCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTggg                                                                                     <  2:919362/150‑1 (MQ=255)
                                                             gCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTggg                                                                                     <  2:341122/150‑1 (MQ=255)
                                                             gCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTggg                                                                                     >  2:78268/1‑150 (MQ=255)
                                                              cGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGt                                                                                    <  1:1082132/150‑1 (MQ=255)
                                                               ggATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTa                                                                                   <  1:1320955/150‑1 (MQ=255)
                                                                  taataataTTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAAcc                                                                                <  2:24255/150‑1 (MQ=255)
                                                                    ataataTTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAAGCtg                                                                              >  2:1395526/1‑150 (MQ=255)
                                                                    ataataTTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCtg                                                                              >  2:959089/1‑150 (MQ=255)
                                                                    ataataTTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCtg                                                                              >  1:1343041/1‑150 (MQ=255)
                                                                     taataTTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCtga                                                                             <  1:304852/150‑1 (MQ=255)
                                                                      aataTTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCtgat                                                                            <  1:1225274/150‑1 (MQ=255)
                                                                        tatTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTTAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCtgatga                                                                          >  1:706945/1‑150 (MQ=255)
                                                                        tatTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCtgatga                                                                          >  1:390831/1‑150 (MQ=255)
                                                                           tCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGaagc                                                                       >  2:210821/1‑150 (MQ=255)
                                                                            cAGACCCTGCGTCGGTTCTTCGATGCGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGaagcaa                                                                     >  1:788850/1‑150 (MQ=255)
                                                                               aCCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGaagcaagc                                                                   <  2:130922/150‑1 (MQ=255)
                                                                               aCCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGaagcaagc                                                                   <  2:390831/150‑1 (MQ=255)
                                                                                 ccTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATAGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACATATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGccc                                                                 >  2:270751/1‑150 (MQ=255)
                                                                                     cGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCcagca                                                             <  1:959089/150‑1 (MQ=255)
                                                                                       tCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGa                                                           <  1:911153/150‑1 (MQ=255)
                                                                                                 gATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGaaaa                                                 >  1:1119392/1‑150 (MQ=255)
                                                                                                  aTGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGaaaaa                                                <  1:414978/150‑1 (MQ=255)
                                                                                                     cGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatg                                             <  1:78268/150‑1 (MQ=255)
                                                                                                        aGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcag                                          <  2:1317258/150‑1 (MQ=255)
                                                                                                        aGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcag                                          <  2:233132/150‑1 (MQ=255)
                                                                                                            ggATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaaaa                                      >  2:559858/1‑150 (MQ=255)
                                                                                                            ggATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaaaa                                      >  1:149136/1‑150 (MQ=255)
                                                                                                             gATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaaaat                                     <  1:270516/150‑1 (MQ=255)
                                                                                                                gaagaaAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaattcagaaaatgca                                  <  2:272587/150‑1 (MQ=255)
                                                                                                                gaagaaAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaaaatgca                                  <  2:834266/150‑1 (MQ=255)
                                                                                                                    aaaGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcag                                          <  2:225016/138‑1 (MQ=255)
                                                                                                                    aaaGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcag                                          >  1:225016/1‑138 (MQ=255)
                                                                                                                     aaGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAag                             >  1:158198/1‑150 (MQ=255)
                                                                                                                             gCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaaaatgcag                                 >  2:1142920/1‑138 (MQ=255)
                                                                                                                             gCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaaaatgcag                                 <  1:1142920/138‑1 (MQ=255)
                                                                                                                              ccccATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCgcgc                    <  2:1343041/150‑1 (MQ=255)
                                                                                                                              ccccATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCgcgc                    <  1:840125/150‑1 (MQ=255)
                                                                                                                              ccccATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCgcgc                    >  1:666618/1‑150 (MQ=255)
                                                                                                                                 cATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGc                 <  2:788850/150‑1 (MQ=255)
                                                                                                                                     cgaacGTAAGCTGTGCTTACGATCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGa             >  1:956626/5‑150 (MQ=255)
                                                                                                                                     cgaacGTAAGCTGTGCTTACGATCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGa             >  1:447613/5‑150 (MQ=255)
                                                                                                                                       aacGTAAGCTGTGCTTACGATCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAg           <  2:447613/148‑1 (MQ=255)
                                                                                                                                            atgCACTGCTTACGTTCGTCATCCTTAACGTGATTTAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATTCAGAAAATGCAGGAAGAGATCGCGCTGCTGGAAGTCAcc      <  2:84481/148‑1 (MQ=255)
                                                                                                                                             aGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCg     >  2:1223943/1‑150 (MQ=255)
                                                                                                                                               cacaGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGc   <  1:206515/150‑1 (MQ=255)
                                                                                                                                                acaGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCg  <  2:48104/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCG  >  NC_000913/493888‑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: