Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,580 G→A 9.5% intergenic (+25/+216) ybiT → / ← ybiU ABC‑F family putative regulatory ATPase/DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913857,5800GA 9.5% 769.7 / 51.7 255intergenic (+25/+216)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base G (121/107);  new base A (11/13);  total (134/121)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

Only 100 of 255 total aligned reads displayed.

TAAGTTAATGATGCAGAAGCCGAACATTCTGATCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCG  >  NC_000913/857312‑857841
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                     
                               aTCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGc                                                                                                                                                                                                                                  <  2:270431/275‑1 (MQ=255)
                                             gACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGAc                                                                                                                                                                                                                                                <  2:490189/247‑1 (MQ=255)
                                                                        tGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaaa                                                                                                                                                                                                                                       <  2:483070/229‑1 (MQ=255)
                                                                                                              cGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGGAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCTATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATTCGCTACGCTTATCCTGCCTCCGCAGTCCCTGCAATATATTTAATCTTCAGGACTTTGTCGGCCGGATAAGGCg                                                                                                                                                   >  4:486338/1‑275 (MQ=255)
                                                                                                                             cacaCGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCggg                                                                                                                                                                                                                                                                  <  6:1657558/149‑1 (MQ=255)
                                                                                                                              acacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgggg                                                                                                                                                                                                                                                                 >  7:221881/1‑149 (MQ=255)
                                                                                                                               cacGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgcgc                                                                                                                                                                                                                         >  4:542234/1‑188 (MQ=255)
                                                                                                                                acGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGt                                                                                                                                                                                                                                                                      >  8:484173/1‑142 (MQ=255)
                                                                                                                                acGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGt                                                                                                                                                                                                                                                                >  7:1137551/1‑148 (MQ=255)
                                                                                                                                 cGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGa                                                                                                                                                                                                                                                              <  6:830896/149‑1 (MQ=255)
                                                                                                                                  gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAc                                                                                                                                                                                                                                                             <  6:890601/149‑1 (MQ=255)
                                                                                                                                  gACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAc                                                                                                                                                                                                                                                             <  7:1423046/149‑1 (MQ=255)
                                                                                                                                        aggTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCCCCCCGGAACGCGTGATCGACTTTAGCGGGAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACaaa                                                                                                                         <  2:219583/273‑1 (MQ=255)
                                                                                                                                          gTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtg                                                                                                                                                                                                                                                     >  8:854794/1‑149 (MQ=255)
                                                                                                                                              gTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatga                                                                                                                                                                                                                                                 >  8:1320003/1‑149 (MQ=255)
                                                                                                                                               tAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                                                                                                                                                                                                >  6:310716/1‑133 (MQ=255)
                                                                                                                                               tAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGAc                                                                                                                                                                                                                                                >  5:1001610/1‑149 (MQ=255)
                                                                                                                                                aaGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACa                                                                                                                                                                                                                                               >  5:568098/1‑149 (MQ=255)
                                                                                                                                                   cTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAAt                                                                                                                                                                                                                                            <  5:335542/149‑1 (MQ=255)
                                                                                                                                                   cTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAAt                                                                                                                                                                                                                                            <  5:1391201/149‑1 (MQ=255)
                                                                                                                                                   cTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCa                                                                                                                                                                                           <  2:150511/198‑1 (MQ=255)
                                                                                                                                                     cccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACATGCGCAGTAAAGGGATCGAGTAAGTGCTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACCAATGCAAAATAGCCTGATATGCTACGCTTCTCAGGCCTACACATTACCTGCCATATATAGAATCTGCAGGAGTTTGTAGGCGGGATAAGGCCTTAAAGCCAAATCAAGAATTAACAAAACGCagtgttgat                                                                                                            >  4:219580/1‑267 (MQ=255)
                                                                                                                                                      ccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                                                                                                                                                                                                >  6:1087243/1‑126 (MQ=255)
                                                                                                                                                       cTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                                                                                                                                                                                                        >  6:169514/1‑149 (MQ=255)
                                                                                                                                                               ccGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                                                                                                                                                                                                >  7:389059/1‑149 (MQ=255)
                                                                                                                                                               ccGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                                                                                                                                                                                                >  5:1626616/1‑149 (MQ=255)
                                                                                                                                                               ccGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                                                                                                                                                                                                >  7:825209/1‑149 (MQ=255)
                                                                                                                                                               ccGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                                                                                                                                                                                                >  7:783655/1‑149 (MQ=255)
                                                                                                                                                                cGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGa                                                                                                                                                                                                                              >  4:124450/1‑150 (MQ=255)
                                                                                                                                                                 gCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGa                                                                                                                                                                                                                              >  8:368252/1‑149 (MQ=255)
                                                                                                                                                                 gCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGa                                                                                                                                                                                                                              <  6:804940/149‑1 (MQ=255)
                                                                                                                                                                  cATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGAt                                                                                                                                                                                                                             <  5:1054263/149‑1 (MQ=255)
                                                                                                                                                                   aTTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATg                                                                                                                                                                                                                            <  7:927686/149‑1 (MQ=255)
                                                                                                                                                                        ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTa                                                                                                                                                                                                                       >  8:204649/1‑149 (MQ=255)
                                                                                                                                                                        ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCggg                                                                                                                                                                                                                                                                  <  5:1350025/106‑1 (MQ=255)
                                                                                                                                                                        ggAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCggg                                                                                                                                                                                                                                                                  >  6:1350025/1‑106 (MQ=255)
                                                                                                                                                                         gAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTAc                                                                                                                                                                                                                      <  8:440518/149‑1 (MQ=255)
                                                                                                                                                                         gAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGGGATTATCCGTCACCCCGACAGTATCgg                                                                                                                                                                                                                                                                   >  5:1147140/1‑104 (MQ=255)
                                                                                                                                                                           aaTCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGCAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAAt                                                                                                                                                                                                                                            <  8:958062/125‑1 (MQ=255)
                                                                                                                                                                           aaTCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGc                                                                                                                                                                                                                    >  5:213713/1‑149 (MQ=255)
                                                                                                                                                                             tCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCtt                                                                                                                                                                                                                  >  8:136322/1‑149 (MQ=255)
                                                                                                                                                                               aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                >  8:773953/1‑149 (MQ=255)
                                                                                                                                                                               aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                <  6:373889/149‑1 (MQ=255)
                                                                                                                                                                                 cccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCa                                                                                                                                                                                                    >  3:48276/1‑159 (MQ=255)
                                                                                                                                                                                 cccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCa                                                                                                                                                                                                    <  1:48276/159‑1 (MQ=255)
                                                                                                                                                                                         gggttgtctaCTTTACTGGTAATTGCAGAGATTTCCTGCGTGGGAACGGGGCCGAGTAAGTGATTCACCGTCACCCCGATANGTTCGGGGTGACGGCGGCTGATGACAAATTCAAAATAGCCTGATGCGATACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACAt                                                                    <  3:234039/270‑1 (MQ=37)
                                                                                                                                                                                               tCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCt                                                                                                                                                                                                                                                      <  6:534423/95‑1 (MQ=255)
                                                                                                                                                                                                cGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCa                                                                                                                                                                                                                                         >  8:1277351/1‑107 (MQ=255)
                                                                                                                                                                                                 gACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCt                                                                                                                                                                                              >  8:1181939/1‑149 (MQ=255)
                                                                                                                                                                                                   cTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGc                                                                                                                                                                                            <  7:1679589/149‑1 (MQ=255)
                                                                                                                                                                                                    tttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCa                                                                                                                                                                                           >  7:611654/1‑149 (MQ=255)
                                                                                                                                                                                                     ttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCaa                                                                                                                                                                                          <  6:1438221/149‑1 (MQ=255)
                                                                                                                                                                                                      tAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAAGCACCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTGTTTTAATCCCACATCAACTAACTTcct                                                       >  1:234039/1‑278 (MQ=255)
                                                                                                                                                                                                        gCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa                                                                                                                                                                                  <  4:470215/154‑1 (MQ=255)
                                                                                                                                                                                                          ggTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                                                                                                                                                                                     <  7:136322/149‑1 (MQ=255)
                                                                                                                                                                                                                    aaaaTTACCTGCTGGTTAAAGGGATCGAGTAAGTGATTATCCGTAACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCa                                             <  4:635794/272‑1 (MQ=255)
                                                                                                                                                                                                                         tACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACt                                                                                                                                                                      <  8:739252/149‑1 (MQ=255)
                                                                                                                                                                                                                          aCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACAt                                        >  2:617530/1‑274 (MQ=255)
                                                                                                                                                                                                                            cTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTg                                                                                                                                                                   <  5:666070/149‑1 (MQ=255)
                                                                                                                                                                                                                             tGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGt                                                                                                                                                                  <  7:660577/149‑1 (MQ=255)
                                                                                                                                                                                                                              gCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTa                                                                                                                                                                 <  5:164756/149‑1 (MQ=255)
                                                                                                                                                                                                                               cGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                                                                                                                                                >  3:133406/1‑117 (MQ=255)
                                                                                                                                                                                                                               cGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                                                                                                                                                                                                <  1:133406/117‑1 (MQ=255)
                                                                                                                                                                                                                                 tagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAcc                                                                                                                                                                                               >  4:741193/1‑116 (MQ=255)
                                                                                                                                                                                                                                    taAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                                                                                                                                                                                           <  8:797871/101‑1 (MQ=255)
                                                                                                                                                                                                                                    taAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa                                                                                                                                                                                  <  6:1269407/126‑1 (MQ=255)
                                                                                                                                                                                                                                    taAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa                                                                                                                                                                                  >  5:1269407/1‑126 (MQ=255)
                                                                                                                                                                                                                                    taAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGt                                                                                                                                                                  <  7:1449769/142‑1 (MQ=255)
                                                                                                                                                                                                                                     aaaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCa                                                                                                                                                                                           <  2:741199/116‑1 (MQ=255)
                                                                                                                                                                                                                                      aaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                                                                                                                                       <  6:511660/103‑1 (MQ=255)
                                                                                                                                                                                                                                      aaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                                                                                                                                       >  5:511660/1‑103 (MQ=255)
                                                                                                                                                                                                                                       aGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGa                                                                                                                                                          <  8:577347/147‑1 (MQ=255)
                                                                                                                                                                                                                                       aGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATa                                                                                                                                                        >  8:145901/1‑149 (MQ=255)
                                                                                                                                                                                                                                           aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                                                                                                                                                    >  5:26542/1‑149 (MQ=255)
                                                                                                                                                                                                                                           aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                                                                                                                                                    >  5:967398/1‑149 (MQ=255)
                                                                                                                                                                                                                                            tCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAg                                                                                                            >  4:201119/1‑188 (MQ=255)
                                                                                                                                                                                                                                             cGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGt                                                                                                                                                  <  7:497596/149‑1 (MQ=255)
                                                                                                                                                                                                                                                gTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACttt                                                                                                                                                                    <  7:817573/128‑1 (MQ=255)
                                                                                                                                                                                                                                                gTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTaa                                                                                                                                               <  8:1555604/149‑1 (MQ=255)
                                                                                                                                                                                                                                                gTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTaa                                                                                                                                               >  5:842499/1‑149 (MQ=255)
                                                                                                                                                                                                                                                gTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTaa                                                                                                                                               <  7:1290116/149‑1 (MQ=255)
                                                                                                                                                                                                                                                gTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAAt                                                                                                        <  2:201122/188‑1 (MQ=255)
                                                                                                                                                                                                                                                 tAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATACGGCGTTAAc                                                                                                                                              >  5:740782/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     tGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgccg                                                                                                                                          >  8:821456/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          aTCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGttta       >  2:696398/1‑273 (MQ=255)
                                                                                                                                                                                                                                                           tCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAg                                                                                                                                    <  6:213713/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             cGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGc                                                                                                                                  <  8:1386737/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               tCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAt                                                                                                                                <  5:294606/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               tCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCTCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTCCCGTGCATGTTGAGGTCg  >  2:452234/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                cACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAtt                                                                                                                               <  7:1268637/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aCCCCGACAGTATCGGGGTGCCGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa                                                                                                                              >  6:68501/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 aCCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCa                                                                                                                                                                           >  8:1611589/1‑104 (MQ=255)
                                                                                                                                                                                                                                                                 aCCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCa                                                                                                                                                                           <  7:1611589/104‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aCCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                                                                                                                                                    <  8:1064962/127‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aCCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTa                                                                                                                              >  6:930917/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 aCCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATAAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGc                                                                                                                                                    >  7:1064962/1‑127 (MQ=255)
                                                                                                                                                                                                                                                                  ccccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTaa                                                                                                                             <  7:1320003/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cccGATACTGTCTGCGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCg                                                                                                                                                   >  7:125640/1‑126 (MQ=255)
                                                                                                                                                                                                                                                                    ccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACa                                                                                                                           >  8:1233599/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       aTACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgccgc                                                                                                                                         <  5:8034/132‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                     
TAAGTTAATGATGCAGAAGCCGAACATTCTGATCATGGACGAACCGACCAACCACCTGGATATGGAATCCATTGAGTCGCTGAACATGGCACTGGAACTGTATCAGGGCACGCTGATCTTTGTTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCG  >  NC_000913/857312‑857841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: