Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,147,260 C→T 5.8% G41G (GGC→GGT frwD → putative enzyme IIB component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,147,2600CT 5.8% 2419.6 / 79.0 737G41G (GGC→GGTfrwDputative enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base C (343/348);  new base T (21/22);  total (365/372)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 771 total aligned reads displayed.

ACTTCAGGTTACCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGC  >  NC_000913/4147112‑4147408
                                                                                                                                                    |                                                                                                                                                    
aCTTCAGGTTACCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGgc                                                                                                                                                      >  1:5749848/1‑149 (MQ=255)
  ttCTGGTTACCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACCTATATGGCGGCGGAACGGCTGGAAAAGTGGTGCCTGTTAGAGAAGTGGGGAGGCAGCATTGAAACTCAGGgcgc                                                                                                                                                    >  4:1136571/1‑149 (MQ=255)
  ttCAGGTTACCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGgcgc                                                                                                                                                    >  3:4004940/1‑149 (MQ=255)
  ttCAGGTTACCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGgcgc                                                                                                                                                    >  3:6706850/1‑149 (MQ=255)
  ttCAGGTTACCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGgcgc                                                                                                                                                    >  4:1643722/1‑149 (MQ=255)
      ggTTACCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTg                                                                                                                                                >  1:976310/1‑149 (MQ=255)
      ggTTACCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTg                                                                                                                                                <  1:4652263/149‑1 (MQ=255)
       gTTACCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTgg                                                                                                                                               >  3:4206659/1‑149 (MQ=255)
       gTTACCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTAGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTgg                                                                                                                                               <  3:4386642/149‑1 (MQ=255)
       gTTACCGTGGGAGGTTAAAATGGCATACCTAGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTgg                                                                                                                                               <  1:5074326/149‑1 (MQ=255)
          aCCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGaa                                                                                                                                            <  4:606166/149‑1 (MQ=255)
           ccGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAAc                                                                                                                                           <  2:6952737/149‑1 (MQ=255)
            cGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACg                                                                                                                                          >  1:4308290/1‑149 (MQ=255)
             gTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACgg                                                                                                                                         <  1:1986166/149‑1 (MQ=255)
              tGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGga                                                                                                                                        <  3:3797663/149‑1 (MQ=255)
              tGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGga                                                                                                                                        >  4:3647563/1‑149 (MQ=255)
              tGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGga                                                                                                                                        <  2:3618892/149‑1 (MQ=255)
                 gAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGaa                                                                                                                                     <  3:2175683/149‑1 (MQ=255)
                 gAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGaa                                                                                                                                     <  3:6932202/149‑1 (MQ=255)
                   gttgaaaaTGGCTGACCTGGTGGCAGAAACCGCCTGCGTAAGCTGCGTGACGCATACTTAAATGGCGGCGGAACAGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGTGCGCTGGGAACGGAGAATc                                                                                                                                   <  2:1586054/145‑1 (MQ=255)
                   ggTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATc                                                                                                                                   <  1:3629422/149‑1 (MQ=255)
                   ggTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATc                                                                                                                                   >  4:5798190/1‑149 (MQ=255)
                   ggTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATc                                                                                                                                   <  4:6706850/149‑1 (MQ=255)
                    gTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCg                                                                                                                                  >  4:7087172/1‑149 (MQ=255)
                     ttAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGt                                                                                                                                 <  4:4265318/149‑1 (MQ=255)
                       aaaaTGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGttt                                                                                                                               <  2:5518183/149‑1 (MQ=255)
                         aaTGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAg                                                                                                                             >  1:100220/1‑149 (MQ=255)
                          aTGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGc                                                                                                                            <  2:6027285/149‑1 (MQ=255)
                               aTACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATg                                                                                                                       >  1:338262/1‑149 (MQ=255)
                               aTACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATg                                                                                                                       >  2:450971/1‑149 (MQ=255)
                                tACCTGGTGGCAGAAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGa                                                                                                                      <  4:2385163/149‑1 (MQ=255)
                                  ccTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAgg                                                                                                                    <  4:4864120/149‑1 (MQ=255)
                                     ggtggCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGata                                                                                                                 <  3:3443314/149‑1 (MQ=255)
                                         gCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATccg                                                                                                             >  2:7137625/1‑149 (MQ=255)
                                               aCCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGTGCGCTGGGAACGGAGAATCGTTTAGCTGATGAGGATATCCGCCTGGc                                                                                                       <  2:6523169/149‑1 (MQ=255)
                                                  gCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGa                                                                                                    <  4:2906913/149‑1 (MQ=255)
                                                  gCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGa                                                                                                    >  1:1337452/1‑149 (MQ=255)
                                                   ccTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGAt                                                                                                   <  4:4515162/149‑1 (MQ=255)
                                                     tGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGTGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGt                                                                                                 <  3:2416629/149‑1 (MQ=255)
                                                      gCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGtt                                                                                                <  2:623295/149‑1 (MQ=255)
                                                        gTAAGCGGCGTGGCGCATACTTATATGGTGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGTGCGCTGGGAACGGAGAATCGtt                                                                                                                                <  4:2575321/115‑1 (MQ=255)
                                                        gTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGc                                                                                              >  4:5550962/1‑149 (MQ=255)
                                                          aaGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGctc                                                                                            <  1:6291892/149‑1 (MQ=255)
                                                           aGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGTGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGctct                                                                                           <  4:3190325/149‑1 (MQ=255)
                                                             cggcgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGt                                                                                         <  2:3248882/149‑1 (MQ=255)
                                                              ggcgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGct                                                                                             >  1:6115549/1‑144 (MQ=255)
                                                              ggcgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGtt                                                                                        <  3:2668385/149‑1 (MQ=255)
                                                              ggcgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGtt                                                                                        <  4:4050757/149‑1 (MQ=255)
                                                                  tGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGAtt                                                                                    <  3:3028374/149‑1 (MQ=255)
                                                                   ggCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTa                                                                                   <  2:5826941/149‑1 (MQ=255)
                                                                    gcgcATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTAc                                                                                  <  1:3918910/149‑1 (MQ=255)
                                                                           cTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGTGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGTGGATGTTGCTCTGTTGATTACGGATATc                                                                           <  3:549016/149‑1 (MQ=255)
                                                                            ttATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCg                                                                          >  4:6280597/1‑149 (MQ=255)
                                                                             tataTGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGa                                                                         <  3:4724298/149‑1 (MQ=255)
                                                                              atatGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAg                                                                        <  4:1388924/149‑1 (MQ=255)
                                                                                  ggcggcGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGc                                                                                                                            >  3:1951271/1‑93 (MQ=255)
                                                                                   gcggcgGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGc                                                                   <  2:4268198/149‑1 (MQ=255)
                                                                                        ggAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTg                                                              >  1:761804/1‑149 (MQ=255)
                                                                                            cGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGa                                                          <  1:1849188/149‑1 (MQ=255)
                                                                                            cGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGa                                                          <  2:3366932/149‑1 (MQ=255)
                                                                                                tGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGTGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGa                                                          >  4:5031177/1‑145 (MQ=255)
                                                                                                tGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGTGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGa                                                          <  3:5031177/145‑1 (MQ=255)
                                                                                                tGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGc                                                             <  2:1914401/142‑1 (MQ=255)
                                                                                                tGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGa                                                      >  3:1372333/1‑149 (MQ=255)
                                                                                                 ggAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGAt                                                     >  2:2487728/1‑149 (MQ=255)
                                                                                                  gAAAAGTTGTGCCTGTTAGAGAAGTGGGGATTCAGCATTGAAAGTCAGGGCGCGCTGGGAAc                                                                                                                                           >  2:2278802/1‑62 (MQ=255)
                                                                                                    aaaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGc                                                             >  3:108334/1‑138 (MQ=255)
                                                                                                    aaaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGAtt                                                    >  1:6539969/1‑147 (MQ=255)
                                                                                                    aaaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTg                                                  <  1:3085712/149‑1 (MQ=255)
                                                                                                       gTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGAt                                                     <  4:210295/143‑1 (MQ=255)
                                                                                                        ttGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGTGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACa                                              >  2:48851/1‑149 (MQ=255)
                                                                                                          gtgCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACAtt                                            >  2:5295802/1‑149 (MQ=255)
                                                                                                             ccTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGcc                                         <  1:3004862/149‑1 (MQ=255)
                                                                                                               tGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGc                                       >  2:5174412/1‑149 (MQ=255)
                                                                                                               tGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACAGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGc                                       >  3:1775891/1‑149 (MQ=255)
                                                                                                                 ttAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGt                                                                                                 <  1:28494/89‑1 (MQ=255)
                                                                                                                 ttAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTa                                     >  3:4603448/1‑149 (MQ=255)
                                                                                                                  tAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAt                                    <  1:6950876/149‑1 (MQ=255)
                                                                                                                    gagaAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCGGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAtgt                                  >  2:3922630/1‑149 (MQ=255)
                                                                                                                        aGTGGGGAGTCAGCATTGAAACTCAGGGTGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCaa                              <  1:5325851/149‑1 (MQ=255)
                                                                                                                        aGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCaa                              <  2:3065444/149‑1 (MQ=255)
                                                                                                                        aGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCaa                              <  3:6078092/149‑1 (MQ=255)
                                                                                                                          tGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATg                            >  3:2770237/1‑149 (MQ=255)
                                                                                                                           ggggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGata                                                                             >  1:5865963/1‑99 (MQ=255)
                                                                                                                            gggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAATGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgca                          <  4:4655067/149‑1 (MQ=255)
                                                                                                                               aGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgcagca                       <  1:2381157/149‑1 (MQ=255)
                                                                                                                                gTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCAt                      <  1:151609/149‑1 (MQ=255)
                                                                                                                                 tCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATc                     >  4:3891737/1‑149 (MQ=255)
                                                                                                                                    gCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTAc                  <  1:7125384/149‑1 (MQ=255)
                                                                                                                                    gCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTAc                  <  4:5404909/149‑1 (MQ=255)
                                                                                                                                     cATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACg                 >  3:5670473/1‑149 (MQ=255)
                                                                                                                                      aTTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGc                >  2:4740213/1‑149 (MQ=255)
                                                                                                                                       ttGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCa               >  2:1337359/1‑149 (MQ=255)
                                                                                                                                       ttGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCa               >  1:4789535/1‑149 (MQ=255)
                                                                                                                                       ttGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCa               >  3:967834/1‑149 (MQ=255)
                                                                                                                                          aaaCTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTc            <  4:1939139/149‑1 (MQ=255)
                                                                                                                                              tCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGc        >  3:6413886/1‑149 (MQ=255)
                                                                                                                                              tCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGc        <  2:292378/149‑1 (MQ=255)
                                                                                                                                                 gggCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTg     >  3:1620163/1‑149 (MQ=255)
                                                                                                                                                    cgcgcTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGc  <  2:781423/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTTCAGGTTACCGTGGGAGGTTAAAATGGCATACCTGGTGGCAGTAACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGC  >  NC_000913/4147112‑4147408

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