Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,373,421 C→T A272T (GCG→ACG)  yjeJ ← protein YjeJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,373,4210CT98.5% 200.8 / NA 65A272T (GCG→ACG) yjeJprotein YjeJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/36);  minor base G (0/1);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAATTGTTCTTTGCCACTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTT  >  NC_000913/4373182‑4373667
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
gACAAACGTATAATGGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTACGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCAt                                                                                                                                                                                                                                              <  1:109232/250‑1 (MQ=255)
       ttattaTTGCCTGATGCTCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGCATGCTTGAGAATTTGTCGGCCGGCTAAGCCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGAAGTTATGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATccgg                                                                                                                                                                                                                                      >  2:175551/6‑249 (MQ=255)
        tataATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTgg                                                                                                                                                                                                                                     >  1:70928/1‑251 (MQ=255)
          taATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGt                                                                                                                                                                                                                                   >  1:128244/1‑251 (MQ=255)
              tGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAg                                                                                                                                                                                                                                <  1:475314/250‑1 (MQ=255)
               gccggATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCg                                                                                                                                                                                                                              <  2:298873/247‑1 (MQ=255)
               gCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGc                                                                                                                                                                                                                               <  1:255740/250‑1 (MQ=255)
                ccTGATGCGCTGCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCgg                                                                                                                                                                                                                             <  1:21422/251‑1 (MQ=255)
                    aTGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTgg                                                                                                                                                                                                                         >  1:451949/1‑251 (MQ=255)
                         cttcgcttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGCTGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGAGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGtt                                                                                                                                                                                                                    >  2:502512/1‑251 (MQ=255)
                         cttcgcttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGACCGGATAAGGCGTTTACGCAGCATCCGGCATGACCAACGATCACGTTATCAGTAATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGtt                                                                                                                                                                                                                    >  1:526680/1‑251 (MQ=255)
                             gcttATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTc                                                                                                                                                                                                                >  1:135306/1‑251 (MQ=255)
                                 aTCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTcgcc                                                                                                                                                                                                             <  1:218504/250‑1 (MQ=255)
                                       ccTACATAAGGCAGTGCAAGTGATTGAATGATTGAGAATTTGCAGGCCGGATAAGGCGTTTAAGCCGCATCAGGCATGAACAACGATAACGCTATCAGTCATTTCGGGCTGGATGGCCCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGATTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGAACGAATGTATCCTAGGTTAGCGGTTGGGCGTTATTCGCCGCTACgg                                                                                                                                                                                                      >  2:296450/1‑251 (MQ=255)
                                         tACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCAGTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTAAGGTg                                                                                                                                                                                                    >  2:432037/1‑251 (MQ=255)
                                                    ttGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGACCATAGAc                                                                                                                                                                                         >  1:63069/1‑251 (MQ=255)
                                                       cAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTg                                                                                                                                                                                      >  1:519105/1‑251 (MQ=255)
                                                        aaCTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGa                                                                                                                                                                                     <  1:354371/251‑1 (MQ=255)
                                                            gATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACAt                                                                                                                                                                                  <  1:249440/250‑1 (MQ=255)
                                                             attgaacgagtGAGAATTTGAAGACAGGATAAGGCGTTTACGAAGAATCCGGCCTGAACAACGATCACGTTATCAGTCAGTTAAGGATGAATGGCTCCAGCCTCGCATTACAAAAAATGAGATTATTTCGACTGTGGTGCAGTTTTGCTGGTGGCTTGAACACGTAAGAGATGCAGCGTACGTACGAGATGATCCTGGGTTAGCGGGTGGACGTTATTCGCCGCTACGGTGCGAACATAGACCtacaatgg                                                                                                                                                                                >  2:83999/11‑244 (MQ=255)
                                                                          gaATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGt                                                                                                                                                                   <  1:353270/251‑1 (MQ=255)
                                                                           aaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTa                                                                                                                                                                                                                                 <  1:131010/188‑1 (MQ=255)
                                                                           aaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTa                                                                                                                                                                                                                                 >  2:131010/1‑188 (MQ=255)
                                                                           aaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGtt                                                                                                                                                                  >  1:359139/1‑251 (MQ=255)
                                                                              ttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGtttt                                                                                                                                                                <  1:426003/250‑1 (MQ=255)
                                                                                       ggATAATGCGTTTACGCCGTATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTTCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGTGGGTGGGCGTTATTCGCCGCTACGGTGCGAGCATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTg                                                                                                                                                      <  2:233142/251‑1 (MQ=255)
                                                                                         aTAAGGCGTTTGCGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCCTTTCAGGCTGGATGTCTCCAGCCTCGCATTAGCAAAGATGAGATTATTCCGCCTGTGGGGCAGTTTGGCTGTTGGATTGAACACGTAAGGGATGCTGCGGGCGTAAGCATTGATCCTGGGTTATCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTgg                                                                                                                                                     <  2:369773/250‑1 (MQ=255)
                                                                                              gCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATACTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGa                                                                                                                                                >  2:87939/1‑250 (MQ=255)
                                                                                                     cgccgcATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACAGGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGc                                                                                                                                        <  1:495889/251‑1 (MQ=255)
                                                                                                      gccgtATCCGGCATGAACAACGATCACGTTATCAGTCATTTCAGGCTGGATGGTCCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTAGTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCa                                                                                                                                       <  2:69030/251‑1 (MQ=255)
                                                                                                          cATCCGGCTTGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATGTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCg                                                                                                                                   <  2:363836/251‑1 (MQ=255)
                                                                                                           aTCCGGCATGACCAACGTTCACGTTATCAGTCAGTTCAGGCTGGATGGTGCCAGCGTCGCAGTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTGAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCgg                                                                                                                                  <  2:58831/251‑1 (MQ=255)
                                                                                                                 cctgacccacgaataCGTTATCAGGCACTTCAGGCTGGATGGCTCCAGGCTCGCATTAGCAAAGCTGAGATTATTTCGCTCGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGGTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTg                                                                                                                            <  2:377972/237‑1 (MQ=255)
                                                                                                                     cccAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCg                                                                                                                        <  2:359139/250‑1 (MQ=255)
                                                                                                                       caacgaataCGTTATAAGTCATTTAAGACTGGATGGCTCTAGCCACGCATTAGCAAAAATGAGCTTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGCACCCGTAAGAGATGCAGCGTGCGGAAGATTTGATCCTGGGTTAGCGGTTGGGCGTTAGTCGCCGCTACGGCGCGACCATAGCCCTGACATGGCACCCCCGCGCGTTTGTTCCGTCGTGGTATGAAATCAAGATTTCGGCgggcagtcgat                                                                                                                      >  2:353631/9‑241 (MQ=255)
                                                                                                                           gATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCt                                                                                                                   >  1:95918/1‑250 (MQ=255)
                                                                                                                              cACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTGTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTAttt                                                                                                               <  1:348728/251‑1 (MQ=255)
                                                                                                                                       aGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATGATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCg                                                                                                      <  1:314986/251‑1 (MQ=255)
                                                                                                                                          aaTTTCAGGCGGGATGGCTAAAGCCTCGCATTAACAAAGATGAGATTATTTAGACTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAATTGATGCAGCGTGCGTAAGAATTGAGCCTGGGTTAGCGGTTGGGCTTTATTCGCCGCTACGGTGCGAACATAGACCTGACAGGGCACACCCGCGAGTTTTTTTAGTCGTGGGCTGAAATCAAGCATACGGCGAGTGATCGCTTCTAGTTCTTTACCCGGCg                                                                                                   <  2:328812/250‑1 (MQ=255)
                                                                                                                                              gcAGGCGGGATGGCTCCAGTCTAGAATTAGCAAAGAGGAGATTATTTCGCCTGTGGTGCAGTTTGGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACAGAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTg                                                                                                <  2:135306/249‑1 (MQ=255)
                                                                                                                                                gggcgggATGGCTCCAGCCTCGCTTTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCGGGTGGATGGAACACGTAAGTGAGGCAGCTTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCtt                                                                                             <  2:519105/246‑1 (MQ=255)
                                                                                                                                                 ggCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTc                                                                                            <  1:323129/251‑1 (MQ=255)
                                                                                                                                                             ccAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTc                                                                                  <  1:483950/249‑1 (MQ=255)
                                                                                                                                                               aGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTcacca                                                                              <  1:502512/251‑1 (MQ=255)
                                                                                                                                                               aGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTcacca                                                                              <  1:246610/251‑1 (MQ=255)
                                                                                                                                                                gCCTCGCATTAGCAAAGATGAGATTATTGCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTcaccac                                                                             <  1:441400/251‑1 (MQ=255)
                                                                                                                                                                    cGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGc                                                                                                                                                                                                                               >  2:444954/1‑101 (MQ=255)
                                                                                                                                                                    cGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGc                                                                                                                                                                                                                               <  1:444954/101‑1 (MQ=255)
                                                                                                                                                                      cATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGcgcc                                                                       <  1:175551/251‑1 (MQ=255)
                                                                                                                                                                      cATTAGCAAAGATGAGATTATTGCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGACCATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGcgcc                                                                       <  2:238884/251‑1 (MQ=255)
                                                                                                                                                                              aaGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAAtt                                                               >  2:475483/1‑251 (MQ=255)
                                                                                                                                                                                gattagaTTATTGCGCCTGTGGTGCAGTGTTGCTGTTGGTTTGAACACGTAATGTATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTAGTCGCCGCTACGGCGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGTGCTGAAATCAAGCATGCGGCGAGTGATCGCTTCTATTTCTTTTCCCGGCGTTGCTTCACGCGTTTTCGCCACTGCGCCGCTAAATTGt                                                             <  1:348412/247‑1 (MQ=255)
                                                                                                                                                                                    agaTTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTACATTGTTCtt                                                         >  2:156234/1‑251 (MQ=255)
                                                                                                                                                                                      attagtTCGCCTGTGGTGCAGTGTTGCTGGTGGATGGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATGCTGGGGTAGCGGTTGGGCGTTATTCGCCGCTACGGTTCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACTCGTGTTCTCCACTGCGCCGCTAAATTGTTCTTTg                                                       <  2:63069/251‑1 (MQ=255)
                                                                                                                                                                                            tCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAATTGTTCTTTGCCACtt                                                 >  1:497297/1‑251 (MQ=255)
                                                                                                                                                                                             cGCCTGTGGTGCAGTTATGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGGTCCTGGGTTAGCGGTTGGGAGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCAAACCCGCGAGTTTTTTAAGTAGTGGGCTGAAATAAAGCATTAGGCGGGTGATCGCTTCTATTTCTTTACCCGGCGATACTTCCCGCGTTATCACCACTGCGCCGCTAAATTGTTTTTTGCCCCttt                                                >  2:453537/1‑251 (MQ=255)
                                                                                                                                                                                              gCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAATTGTTCTTTGCCACTTTc                                               >  1:498183/1‑251 (MQ=255)
                                                                                                                                                                                                  gtggtgCGGTTTTGCTGGTGGATTGAACACGTAAGTGCGGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGCAATCAAGCATTCGGCGAGTGATCGCTTCTTTTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAATTTGTTCTTTGCCACTTTCATCt                                           >  2:88971/1‑251 (MQ=255)
                                                                                                                                                                                                               gcttgggtcttGAACACGGAAGTGATGTAGCGTGCGTAAGCATTGAGCCTGGTTTCGTGGTTGGGCGTTATGCGCCGATTCGGTTCGTAAATAGACCTGACAGGGCACGCACGCGATTTTTTGCTGTTGTGGGCTGAAATCAATCATTCTGCGATTTATCGCGTCTAGTTCTTTTCCCGGCGTTGCTTTACTCGTTTTCACCACTGCGCCGCCAAATTGTGCTTTGCCACTTTCATCTTTAAAGCGGTaaa                              <  2:48823/242‑1 (MQ=255)
                                                                                                                                                                                                                 tggtggATTGAACACGTAAGTGATGCAGCGTGCGTAAGCATTGTTCCGGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACACGGCACCCCCGCGAGCTTTTTCAGTCGTGGGCTGAAATCAATGATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAATTGTTCTTTGCCACTTTCATCTTTAAACCGGTAAAGc                            >  2:119627/1‑251 (MQ=255)
                                                                                                                                                                                                                             cacGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAATTGTTCTTTGCCACTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATa                 >  1:198584/1‑250 (MQ=255)
                                                                                                                                                                                                                              acGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTCTTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCCCCACTGCGCCGCTAAATTGTTCTTTGCCACTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAAt               >  2:259602/1‑251 (MQ=255)
                                                                                                                                                                                                                              acGTAAGTGATGCAGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAATTGTTCTTTGCCACTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAAt               <  1:87939/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aTGCAGCGTGCGTAAGCATTGATACTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGCCATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGCGTGATCGCTTCTATTTCTTTACCCGGCGATGCTTCACGCGTTTTCCCCACTGAGCCGATAAATTGTTCTTTGCCACTTTCATCTTTAAAGCGGTAAAGCACAGACAGATAATGGTCGAAc       >  2:44202/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aGCGTGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAATTGTTCTTTGCCACTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGtt  <  1:281742/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
GACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTACGGTGCGAACATAGACCTGACATGGCACACCCGCGAGTTTTTTCAGTCGTGGGCTGAAATCAAGCATTCGGCGAGTGATCGCTTCTATTTCTTTACCCGGCGTTGCTTCACGCGTTTTCACCACTGCGCCGCTAAATTGTTCTTTGCCACTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTT  >  NC_000913/4373182‑4373667

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: