Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 494,042 T→A 8.6% intergenic (+19/‑34) 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,0420TA 8.6% 241.2 / 11.5 81intergenic (+19/‑34)dnaX/ybaBDNA polymerase III subunit tau/putative nucleoid‑associated protein YbaB
Reads supporting (aligned to +/- strand):  ref base T (32/42);  new base A (2/5);  total (34/47)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGC  >  NC_000913/493895‑494189
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                                                  >  2:911153/1‑150 (MQ=255)
                 cGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGAt                                                                                                                                  >  1:83835/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: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                                                                                                           <  1:157151/150‑1 (MQ=255)
                                        gcgAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGt                                                                                                           <  2:1162399/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:341122/150‑1 (MQ=255)
                                                      gCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTggg                                                                                             <  2:919362/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                                              >  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                            >  1:666618/1‑150 (MQ=255)
                                                                                                                       ccccATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCgcgc                            <  1:840125/150‑1 (MQ=255)
                                                                                                                       ccccATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCgcgc                            <  2:1343041/150‑1 (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)
                                                                                                                                               tACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGAAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTg    >  1:708894/1‑150 (MQ=255)
                                                                                                                                                aCGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTgg   >  2:1323375/1‑150 (MQ=255)
                                                                                                                                                aCGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTgg   <  1:642341/150‑1 (MQ=255)
                                                                                                                                                 cGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGgc  >  2:208122/1‑150 (MQ=255)
                                                                                                                                                 cGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCAAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGgc  >  1:1173107/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGCGTCAGGCGATATACGAAGAAAAACTTGCGCAGGCGCGCGAGTCCATTATTGCGGATAATAATATTCAGACCCTGCGTCGGTTCTTCGATGCGGAGCTGGATGAAGAAAGTATCCGCCCCATTTGATCGTAAGCACAGCTTACGTTCGTCATCCTTAACGTGATTGAGAGAGAAACCTATGTTTGGTAAAGGCGGTCTGGGTAACCTGATGAAGCAAGCCCAGCAGATGCAAGAAAAAATGCAGAAAATGCAGGAAGAGATCGCGCAGCTGGAAGTCACCGGCGAATCTGGC  >  NC_000913/493895‑494189

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: