Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,906,607 C→T 5.6% intergenic (+24/+36) ygcG → / ← eno TPM domain protein, putative phosphatase/enolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,906,6070CT 5.6% 345.2 / 6.7 107intergenic (+24/+36)ygcG/enoTPM domain protein, putative phosphatase/enolase
Reads supporting (aligned to +/- strand):  ref base C (42/59);  new base T (4/2);  total (46/61)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 111 total aligned reads displayed.

AGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATA  >  NC_000913/2906459‑2906754
                                                                                                                                                    |                                                                                                                                                   
aGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGccc                                                                                                                                                     <  2:353142/149‑1 (MQ=255)
    aaGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGGGGTGGCGGCGCATCCGGCCGCTGGGAACTCAACCATAAAAAAAGCCCGCCCGGAg                                                                                                                                                 >  2:322980/1‑149 (MQ=255)
    aaGAGCATCGGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGAGGTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCAGCTGGTAACTCAACCATAAATAATGCCAGTCCGGAg                                                                                                                                                 <  2:518276/149‑1 (MQ=255)
      gagCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGc                                                                                                                                               <  2:981079/149‑1 (MQ=255)
       agCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCt                                                                                                                                              <  2:20714/149‑1 (MQ=255)
       agCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCt                                                                                                                                              <  2:38886/149‑1 (MQ=255)
        gCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTg                                                                                                                                             <  1:1023606/149‑1 (MQ=255)
        gCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTg                                                                                                                                             <  2:851915/149‑1 (MQ=255)
           tCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCa                                                                                                                                          >  1:210504/1‑149 (MQ=255)
             tGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGGTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAtt                                                                                                                                        >  1:592111/1‑149 (MQ=255)
              ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAttt                                                                                                                                       <  3:362947/149‑1 (MQ=255)
              ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAttt                                                                                                                                       >  3:271460/1‑149 (MQ=255)
               gTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCAtttt                                                                                                                                      <  2:1048492/149‑1 (MQ=255)
                tAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGGTCTTTTGGAGGTGGCGGCGCATCCGGCCGCTGGGAACTCAACCATAAAAAAAGCCAGCCCGGAGGCGGGCAttttt                                                                                                                                     >  4:339629/1‑149 (MQ=255)
                  gAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTaa                                                                                                                                   <  4:325658/149‑1 (MQ=255)
                   aaGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTaaa                                                                                                                                  <  4:600441/149‑1 (MQ=255)
                   aaGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTaaa                                                                                                                                  <  4:928340/149‑1 (MQ=255)
                     gTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATc                                                                                                                                >  3:56439/1‑149 (MQ=255)
                        aTCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGa                                                                                                                             >  3:350998/1‑149 (MQ=255)
                           aTTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCgg                                                                                                                                                   >  3:681070/1‑124 (MQ=255)
                           aTTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaa                                                                                                                          >  1:457015/1‑149 (MQ=255)
                            ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCCTCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaaa                                                                                                                         >  4:54155/1‑149 (MQ=255)
                            ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaaa                                                                                                                         <  2:330935/149‑1 (MQ=255)
                            ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCAGCCGCTGGGAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaaa                                                                                                                         >  1:121428/1‑149 (MQ=255)
                            ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGGGGCGGCGCGTCCGGCCGCTGGGAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaaa                                                                                                                         >  3:396544/1‑149 (MQ=255)
                            ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGCGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaaa                                                                                                                         <  2:886933/149‑1 (MQ=255)
                            ttCAGACAACAGCGGTTCCGCTGGAGGCCCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGGATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaaa                                                                                                                         <  4:772340/149‑1 (MQ=255)
                             tCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAag                                                                                                                        <  1:553800/149‑1 (MQ=255)
                              cAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAagt                                                                                                                       >  1:964453/1‑149 (MQ=255)
                              cAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGGAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAagt                                                                                                                       >  3:199537/1‑149 (MQ=255)
                                gACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAagtca                                                                                                                     <  4:271460/149‑1 (MQ=255)
                                 acaacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAagtcag                                                                                                                    >  3:754698/1‑149 (MQ=255)
                                   aacaGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGAGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCCGCCCGGCGGATGGCATTTTTAAATCAGATAAagtcagcc                                                                                                                  >  1:36112/1‑147 (MQ=255)
                                       gCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTAt                                                                                                              <  2:536847/149‑1 (MQ=255)
                                       gCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTAt                                                                                                              <  2:648840/149‑1 (MQ=255)
                                          gTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGcc                                                                                                           >  3:301867/1‑149 (MQ=255)
                                           ttCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTTTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCt                                                                                                          >  3:375311/1‑149 (MQ=255)
                                           ttCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCt                                                                                                          <  3:664732/149‑1 (MQ=255)
                                           ttCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGGTCTTTTGGTGGGGGCGGCGCACCCGGCCGCTGGTAACTCCACCCTAAAAAATGCCAGACCGGGGGGTTGCATTTTTAAATTCGATAAAGTGACTCTTAAGCCt                                                                                                          >  3:764121/1‑149 (MQ=255)
                                              cGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGc                                                                                                       <  3:339629/149‑1 (MQ=255)
                                                    aGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGa                                                                                                 <  3:846604/149‑1 (MQ=255)
                                                    aGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGa                                                                                                 <  3:214161/149‑1 (MQ=255)
                                                     ggCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAt                                                                                                <  3:795987/149‑1 (MQ=255)
                                                      gCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGAtc                                                                                               <  2:210504/149‑1 (MQ=255)
                                                          aGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGc                                                                                                                                           >  4:347395/1‑101 (MQ=255)
                                                          aGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGc                                                                                                                                           <  3:347395/101‑1 (MQ=255)
                                                           gATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCttt                                                                                          <  1:182045/149‑1 (MQ=255)
                                                            aTAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTa                                                                                         <  1:46784/149‑1 (MQ=255)
                                                                    ggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGtt                                                                                 <  2:964453/149‑1 (MQ=255)
                                                                    ggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGtt                                                                                 <  1:700729/149‑1 (MQ=255)
                                                                    ggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGtt                                                                                 <  4:199537/149‑1 (MQ=255)
                                                                    ggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGtt                                                                                 >  4:861341/1‑149 (MQ=255)
                                                                      aTTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGGGGCGCATCCGGCCGCTGGTAACTCAACCATATAAAATGCCAACCCGGAGGCTAGCATTTTTAAATCCGATAAAGTCAGTCGTATGCCTGGCCTTTGATCTCTTTACGACCGgtgt                                                                               <  3:153382/149‑1 (MQ=255)
                                                                        ttAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTAc                                                                             <  2:819033/149‑1 (MQ=255)
                                                                            cggcggGGGCGGTTCTACTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAAACAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTg                                                                         <  3:259434/149‑1 (MQ=255)
                                                                              gcggGGGCGGTTCTTCTGGAGGTGGCGGCGCATTCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGTGGCTAGCATTTTTAAATCAGCTAACGTCAGTCTTATGCCTGGCCGTTGATCTCGTTGCGACCGTTGTCCGGTGCt                                                                       <  2:402254/149‑1 (MQ=255)
                                                                              gcggGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCt                                                                       <  4:398983/149‑1 (MQ=255)
                                                                                gggggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCttt                                                                     <  1:175832/149‑1 (MQ=255)
                                                                                 ggggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCttat                                                                    <  4:754698/149‑3 (MQ=255)
                                                                                   ggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTc                                                                  <  1:909426/149‑1 (MQ=255)
                                                                                      tgTTCTTATGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGcc                                                               <  4:764121/148‑1 (MQ=255)
                                                                                      ggTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGcc                                                               >  2:307973/1‑149 (MQ=255)
                                                                                       gTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGccc                                                              >  4:353986/1‑149 (MQ=255)
                                                                                        ttcttcTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCa                                                             >  1:382523/1‑149 (MQ=255)
                                                                                           ttctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCagag                                                          >  3:461899/1‑149 (MQ=255)
                                                                                               ggAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAACAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttc                                                      <  4:751645/149‑1 (MQ=255)
                                                                                               ggAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttc                                                      <  4:91589/149‑1 (MQ=255)
                                                                                                gAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttct                                                     <  3:175246/149‑1 (MQ=255)
                                                                                                 aGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGcttctt                                                    <  2:614187/149‑1 (MQ=255)
                                                                                                  ggTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCttcttc                                                   <  4:170423/149‑1 (MQ=255)
                                                                                                    tGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGa                                                 <  2:1034492/149‑1 (MQ=255)
                                                                                                    tGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGa                                                 <  3:353986/149‑1 (MQ=255)
                                                                                                       cggcgCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATAc                                              <  3:283354/149‑1 (MQ=255)
                                                                                                              tCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAg                                       >  3:461319/1‑149 (MQ=255)
                                                                                                                 ggCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTg                                    <  2:592111/149‑1 (MQ=255)
                                                                                                                      cTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGt                               <  1:402976/149‑1 (MQ=255)
                                                                                                                       tGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTa                              <  1:259809/149‑1 (MQ=255)
                                                                                                                       tGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTa                              >  2:391556/1‑149 (MQ=255)
                                                                                                                         gTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTAtt                            >  3:543153/1‑149 (MQ=255)
                                                                                                                           aaCTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTa                          <  3:583682/149‑1 (MQ=255)
                                                                                                                            gagttaCCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAg                         <  4:499953/144‑1 (MQ=255)
                                                                                                                              tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTc                                                                  >  4:478171/1‑106 (MQ=255)
                                                                                                                              tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCa                       >  4:745489/1‑149 (MQ=255)
                                                                                                                              tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCa                       >  1:528369/1‑149 (MQ=255)
                                                                                                                              tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCa                       >  4:105979/1‑149 (MQ=255)
                                                                                                                              tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCa                       >  3:277002/1‑149 (MQ=255)
                                                                                                                                 tgatttaaaaaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGtt                                 >  4:270839/7‑136 (MQ=38)
                                                                                                                                 tgatttaaaaaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGtt                                 <  3:270839/130‑1 (MQ=38)
                                                                                                                                  gatttaaaaaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAacac                   >  1:623398/6‑149 (MQ=39)
                                                                                                                                   cATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACg                  <  4:629133/149‑1 (MQ=255)
                                                                                                                                    aTAAAAAATGCCAGCCCAGAGGCTGGCATTTTTAAATCAGAAACAGTCAATCTTATGCATGGCCTTTGAACTCTTTACGACCGTTGTACGGTCCTTTTTCGCCCAGAGCTTCTTCGATGCGAATCAGCTGGTTGTATATAGCAACACgg                 >  4:648715/1‑149 (MQ=255)
                                                                                                                                     tAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGt                >  4:539565/1‑149 (MQ=255)
                                                                                                                                       aaaaaTGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCATTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCCGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCa              >  4:301341/1‑149 (MQ=255)
                                                                                                                                        aaaaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAg             >  3:130820/1‑149 (MQ=255)
                                                                                                                                         aaaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGGCTTATGCCTGGCCTTTGATCTCTTTACGTCCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCATCTGGTTGTATTTAGCAACACGGTCAGa            >  2:784690/1‑149 (MQ=255)
                                                                                                                                           aTGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAAc          <  2:187606/149‑1 (MQ=255)
                                                                                                                                               cAGCCCGGTGGCTGGAATTTTTAACTCAGATAAAGTAATTCTTATTCCTGGCCTTTGATGTTTTTACTACCGTTGTACCGTGCTTTTTCGCCCAGAGCTTCTTCGTTACGAATCAGCTCGTTGTATTTAGCAACACGGTCAGAACGGCt      >  4:874228/1‑149 (MQ=255)
                                                                                                                                                  cccGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCAt   <  2:121428/149‑1 (MQ=255)
                                                                                                                                                  cccGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCAt   <  2:17977/149‑1 (MQ=255)
                                                                                                                                                   ccGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATa  <  1:244196/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATA  >  NC_000913/2906459‑2906754

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