Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,966,303 G→T 7.8% A711E (GCA→GAA)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,966,3030GT 7.8% 397.3 / 17.5 115A711E (GCA→GAA) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base G (58/48);  new base T (5/4);  total (63/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 115 total aligned reads displayed.

ACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGA  >  NC_000913/2966159‑2966447
                                                                                                                                                |                                                                                                                                                
aCCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTc                                                                                                                                              <  3:242952/149‑1 (MQ=255)
  cGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCgg                                                                                                                                            <  2:166906/149‑1 (MQ=255)
  cGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCgg                                                                                                                                            >  2:345723/1‑149 (MQ=255)
   gTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGc                                                                                                                                           >  2:16286/1‑149 (MQ=255)
   gTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGACCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGc                                                                                                                                           >  2:801155/1‑149 (MQ=255)
    ttAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTTCTTCGGCa                                                                                                                                          >  1:35383/1‑149 (MQ=255)
           atatGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCCCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCaa                                                                                                                                   <  2:392228/149‑1 (MQ=255)
            tatGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAAt                                                                                                                                  <  4:21469/149‑1 (MQ=255)
              tGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTTCTTCGGCATAATCAATgc                                                                                                                                <  1:825754/149‑1 (MQ=255)
              tGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgc                                                                                                                                <  3:97633/149‑1 (MQ=255)
              tGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgc                                                                                                                                <  4:618492/149‑1 (MQ=255)
               gTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcg                                                                                                                               >  1:193644/1‑149 (MQ=255)
                tataTGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcgc                                                                                                                              >  1:275697/1‑149 (MQ=255)
                tataTGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgcgc                                                                                                                              >  2:425087/1‑149 (MQ=255)
                 atatGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGTTGCGCAAGATTTTCTGCTTCGGCATAATCAATGcgcc                                                                                                                             >  1:102966/1‑149 (MQ=255)
                    tGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTTCTTCGGCATAATCAATGCGCCgca                                                                                                                          >  2:695179/1‑149 (MQ=255)
                      aTCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCa                                                                                                                                          >  2:401887/1‑131 (MQ=255)
                      aTCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCa                                                                                                                                          <  1:401887/131‑1 (MQ=255)
                      aTCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCgcagc                                                                                                                        >  1:723261/1‑149 (MQ=255)
                       tCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCgcagca                                                                                                                       >  3:753688/1‑149 (MQ=255)
                         cgcgCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAgg                                                                                                                     >  4:543335/1‑149 (MQ=255)
                         cgcgCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAgg                                                                                                                     >  3:512057/1‑149 (MQ=255)
                          gcgcTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCaa                                                                                                                                   <  2:673512/134‑1 (MQ=255)
                          gcgcTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGt                                                                                                                    <  2:723261/149‑1 (MQ=255)
                            gcTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTAt                                                                                                                  >  3:194828/1‑149 (MQ=255)
                              tataACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTAttt                                                                                                                >  4:176096/1‑149 (MQ=255)
                               ataACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTAtttt                                                                                                               <  4:327424/149‑1 (MQ=255)
                                taACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATCTTg                                                                                                              <  4:788126/149‑1 (MQ=255)
                                     cTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCTAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCcgc                                                                                                         >  2:243648/1‑149 (MQ=255)
                                      tCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCcgcg                                                                                                        >  1:643612/1‑149 (MQ=255)
                                      tCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCcgcg                                                                                                        >  3:182199/1‑149 (MQ=255)
                                      tCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCcgcg                                                                                                        >  3:244559/1‑149 (MQ=255)
                                      tCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCcgcg                                                                                                        >  2:349353/1‑149 (MQ=255)
                                         gcgAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTa                                                                                                     >  3:193724/1‑149 (MQ=255)
                                         gcgAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTa                                                                                                     <  4:743415/149‑1 (MQ=255)
                                            aaTCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAg                                                                                                  <  1:16286/149‑1 (MQ=255)
                                             aTCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGa                                                                                                 >  3:850482/1‑149 (MQ=255)
                                               cAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAAc                                                                                               <  1:296452/149‑1 (MQ=255)
                                                 gcccgcccATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCaa                                                                                                                                   <  3:409074/111‑1 (MQ=255)
                                                   ccgcccATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTcc                                                                                           <  4:753688/149‑1 (MQ=255)
                                                   ccgcccATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTcc                                                                                           >  2:558159/1‑149 (MQ=255)
                                                    cgcccATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCg                                                                                          <  3:333571/149‑1 (MQ=255)
                                                      cccATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGtt                                                                                        <  4:824121/149‑1 (MQ=255)
                                                       acATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTc                                                                                       <  3:80787/148‑1 (MQ=255)
                                                         aTGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCAt                                                                                     >  3:156140/1‑149 (MQ=255)
                                                          tGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATa                                                                                    <  4:193724/149‑1 (MQ=255)
                                                          tGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATa                                                                                    <  4:704788/149‑1 (MQ=255)
                                                                     tCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCg                                                                         >  4:274221/1‑149 (MQ=255)
                                                                      ccATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGa                                                                        >  1:6321/1‑149 (MQ=255)
                                                                      ccATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGa                                                                        >  1:15847/1‑149 (MQ=255)
                                                                      ccATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGGTGTGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGa                                                                        <  2:643612/149‑1 (MQ=255)
                                                                      ccATAAAGGCTGCAACCTAATGGCGAACTTCGGTCGCCAGATGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGa                                                                        >  1:273549/1‑149 (MQ=255)
                                                                       cATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATgc                                                                                                                                >  4:219893/1‑92 (MQ=255)
                                                                         tAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATAc                                                                     <  3:274221/149‑1 (MQ=255)
                                                                           aaGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGTTTCGCCTTAATCATTGCGCCGCAGCAGTTATTTTGCCCGGGCTACAGAACGCCCGTTCATAGACAGATGGCGATAccc                                                                   >  4:35793/1‑149 (MQ=255)
                                                                           aaGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATAccc                                                                   >  1:126822/1‑149 (MQ=255)
                                                                            aGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCa                                                                  >  1:486230/1‑149 (MQ=255)
                                                                             ggCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTTCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATCGCGATACCCaa                                                                 <  1:695179/149‑1 (MQ=255)
                                                                             ggCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCaa                                                                 >  3:222777/1‑149 (MQ=255)
                                                                             ggCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCaa                                                                 >  3:814117/1‑149 (MQ=255)
                                                                             ggCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCaa                                                                 <  1:801155/149‑1 (MQ=255)
                                                                                tGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGcc                                                              >  1:186019/1‑149 (MQ=255)
                                                                                   aaCCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCaa                                                           <  4:194828/149‑1 (MQ=255)
                                                                                    aCCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAAt                                                          <  3:816553/149‑1 (MQ=255)
                                                                                    aCCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAAt                                                          <  3:543335/149‑1 (MQ=255)
                                                                                      cTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGa                                                        <  4:901742/149‑1 (MQ=255)
                                                                                       tGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAg                                                       <  3:176096/149‑1 (MQ=255)
                                                                                       tGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAg                                                       <  1:550426/149‑1 (MQ=255)
                                                                                           ggCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGAtt                                                   <  1:682008/149‑1 (MQ=255)
                                                                                           ggCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGAtt                                                   <  4:395750/149‑1 (MQ=255)
                                                                                           ggCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGAtt                                                   <  2:193644/149‑1 (MQ=255)
                                                                                           ggCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGAtt                                                   >  4:850622/1‑149 (MQ=255)
                                                                                            gCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTTCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTg                                                  >  2:597730/1‑149 (MQ=255)
                                                                                              gAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTTCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGAt                                                                       >  3:295574/1‑126 (MQ=255)
                                                                                              gAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTTCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGAt                                                                       <  4:295574/126‑1 (MQ=255)
                                                                                                aCTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTTCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCAc                                              <  1:597730/149‑1 (MQ=255)
                                                                                                 cTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACg                                             >  1:118676/1‑149 (MQ=255)
                                                                                                  ttCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGc                                            >  3:802020/1‑149 (MQ=255)
                                                                                                   tCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCATTGAGGATTGCCACGca                                           >  4:88030/1‑149 (MQ=255)
                                                                                                   tCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGca                                           >  2:577254/1‑149 (MQ=255)
                                                                                                        cGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATgg                                      <  4:359212/149‑1 (MQ=255)
                                                                                                          ccAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGa                                    >  1:744513/1‑149 (MQ=255)
                                                                                                          ccAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCAGAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGa                                    <  2:275697/149‑1 (MQ=255)
                                                                                                             gTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGa                                    <  4:378739/146‑1 (MQ=255)
                                                                                                                gcgcTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATcgccc                              <  3:185831/149‑1 (MQ=255)
                                                                                                                    ttCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCa                          >  3:667358/1‑149 (MQ=255)
                                                                                                                    ttCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCa                          <  4:850482/149‑1 (MQ=255)
                                                                                                                     tCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATGGACAGATGGTGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCAt                         >  4:912768/1‑149 (MQ=255)
                                                                                                                      ccAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCAtt                        >  3:659553/1‑149 (MQ=255)
                                                                                                                       cAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCAttt                       >  1:769755/1‑149 (MQ=255)
                                                                                                                        aGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTc                      >  4:8991/1‑149 (MQ=255)
                                                                                                                              cGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCAcc                   <  2:293126/146‑1 (MQ=255)
                                                                                                                              cGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCAcc                   >  1:293126/1‑146 (MQ=255)
                                                                                                                                 tGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACaa             >  1:654392/1‑149 (MQ=255)
                                                                                                                                    gcAAGATTTTCTTCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACg          >  3:662666/1‑149 (MQ=255)
                                                                                                                                    gcAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACg          >  3:472427/1‑149 (MQ=255)
                                                                                                                                      aaGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGga        <  2:15847/149‑1 (MQ=255)
                                                                                                                                      aaGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACACATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACgaa        >  1:31639/1‑147 (MQ=255)
                                                                                                                                         aTTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGAt     <  4:182199/149‑1 (MQ=255)
                                                                                                                                            ttCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGa  >  2:299251/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ACCGTTAAAAGATATGTATATGATCCGCGCTATAACCCTCCGCGAATCAGCCCGCCCATGCCGCGACGCTCCATAAAGGCTGCAACCTGATGGCGAACTTCGGTCGCCAGTTGCGCTTCCAGACTACGCTGCGCAAGATTTTCTGCTTCGGCATAATCAATGCGCCGCAGCAGGTATTTTGCCCGCGCTACAGAACGTCCGTTCATAGACAGATGGCGATACCCAAGCCCAATGAGGATTGCCACGCACATGGGATCGCCCGCCATTTCACCGCACAAACGGAGATCGA  >  NC_000913/2966159‑2966447

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