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

CGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCG  >  NC_000913/493884‑494181
                                                                                                                                                    |                                                                                                                                                     
cgccgcTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGcac                                                                                                                                                      <  2:1145202/150‑1 (MQ=255)
    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                                                                                                                              <  1:609207/150‑1 (MQ=255)
                        tataCGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACg                                                                                                                              >  2:206515/1‑150 (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                                                                                                                    <  2:873725/150‑1 (MQ=255)
                                  aaaaCTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTgagag                                                                                                                    >  1:124373/1‑150 (MQ=255)
                                      cTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGaa                                                                                                                <  1:798871/150‑1 (MQ=255)
                                          cgcAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCt                                                                                                            >  2:1148855/1‑150 (MQ=255)
                                           gcAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTa                                                                                                           >  2:667136/1‑150 (MQ=255)
                                           gcAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTa                                                                                                           <  2:1420537/150‑1 (MQ=255)
                                                 gcgcgAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTg                                                                                                     <  2:580240/150‑1 (MQ=255)
                                                  cgcgAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTgg                                                                                                    >  1:333474/1‑150 (MQ=255)
                                                  cgcgAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTgg                                                                                                    <  2:804551/150‑1 (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:78268/1‑150 (MQ=255)
                                                                 gCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTggg                                                                                     <  2:919362/150‑1 (MQ=255)
                                                                 gCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTggg                                                                                     <  2:341122/150‑1 (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:233132/150‑1 (MQ=255)
                                                                                                            aGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcag                                          <  2:1317258/150‑1 (MQ=255)
                                                                                                                ggATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaaaa                                      >  1:149136/1‑150 (MQ=255)
                                                                                                                ggATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaaaa                                      >  2:559858/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                                 <  1:1142920/138‑1 (MQ=255)
                                                                                                                                 gCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAaaaatgcagaaaatgcag                                 >  2:1142920/1‑138 (MQ=255)
                                                                                                                                  ccccATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCgcgc                    <  2:1343041/150‑1 (MQ=255)
                                                                                                                                  ccccATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCgcgc                    >  1:666618/1‑150 (MQ=255)
                                                                                                                                  ccccATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCgcgc                    <  1:840125/150‑1 (MQ=255)
                                                                                                                                     cATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGc                 <  2:788850/150‑1 (MQ=255)
                                                                                                                                         cgaacGTAAGCTGTGCTTACGATCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGa             >  1:447613/5‑150 (MQ=255)
                                                                                                                                         cgaacGTAAGCTGTGCTTACGATCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGa             >  1:956626/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)
                                                                                                                                                    |                                                                                                                                                     
CGCCGCTGGAGTGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCG  >  NC_000913/493884‑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: