Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,549,548 T→G 5.0% intergenic (+142/+98) yfeW → / ← yfeX penicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,549,5480TG 5.0% 314.1 / 6.4 100intergenic (+142/+98)yfeW/yfeXpenicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base T (53/42);  new base G (2/3);  total (55/45)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 104 total aligned reads displayed.

GAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCG  >  NC_000913/2549400‑2549694
                                                                                                                                                    |                                                                                                                                                  
gAAGTAATAATAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                                                                                                                    <  5:94208/149‑1 (MQ=255)
gAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGATGGTATGAGAGGGCGGACAGACGCAATAATGCCTGATGCGCTACGCTTATAAGGCCTACGTGGAGAAAACAATATATTGAGTTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                                                                                                                    >  2:263269/1‑149 (MQ=37)
      ataaAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                              >  6:284098/1‑149 (MQ=255)
      ataaAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                               <  8:219610/148‑1 (MQ=255)
       taaAAAGGCCAGTCGGTAATGTTACTGGCTGGTCAGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                             >  7:109418/1‑149 (MQ=255)
          aaaGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                          >  3:81671/1‑149 (MQ=255)
              gCCAGGCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                      <  2:328489/149‑1 (MQ=21)
                cagtcagtAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                       <  1:226613/146‑1 (MQ=25)
                cagtcagtAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCg                                                                                                                                    >  2:321382/1‑149 (MQ=255)
                cagtcagtAATGTTACTGCCTGGACTGAGAGAGCTGACAAACGCAAAACTGCATGATGCTCTACGCTTATCAGGCATAAGTGGAGCCCACAATATATTGAATTAGCACCATTTCGTAGGCCGGATAAGGCGTTCACGCCGCAAccgccg                                                                                                                                    >  6:117230/1‑149 (MQ=21)
                 agtaagtaATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCgt                                                                                                                                   <  4:219009/149‑1 (MQ=255)
                   tcagtAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCgtg                                                                                                                                  >  2:166182/1‑148 (MQ=255)
                   tcagtAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCgtg                                                                                                                                  >  3:133088/1‑148 (MQ=255)
                       tAATGTTACTGGCTGGTCTGAGATTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAAAATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAAGCCGCATCCGGCGTTTACa                                                                                                                              >  7:85068/1‑148 (MQ=255)
                       tAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACaa                                                                                                                             >  1:136194/1‑149 (MQ=255)
                       tAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCACATCCGGCGTGTACaa                                                                                                                             >  6:258218/1‑149 (MQ=255)
                                 ggctggTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTg                                                                                                                   >  8:257084/1‑149 (MQ=255)
                                  gctggTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACtt                                                                                                                    <  2:96851/147‑1 (MQ=255)
                                  gctggTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTgg                                                                                                                  <  6:341661/149‑1 (MQ=255)
                                   ctggTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGt                                                                                                                 <  5:284098/149‑1 (MQ=255)
                                         tGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTcaaca                                                                                                            >  1:81501/1‑148 (MQ=255)
                                           agagTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATc                                                                                                         >  8:61830/1‑149 (MQ=255)
                                             agTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCt                                                                                                        <  7:216128/148‑1 (MQ=255)
                                              gTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCACATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGa                                                                                                      <  5:258218/149‑1 (MQ=255)
                                               tGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAg                                                                                                     <  1:166182/149‑1 (MQ=255)
                                               tGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAg                                                                                                     >  8:350260/1‑149 (MQ=255)
                                                gCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAgg                                                                                                    >  2:303478/1‑149 (MQ=255)
                                                 cTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGc                                                                                                   >  8:47339/1‑149 (MQ=255)
                                                  tGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTg                                                                                                                   <  7:294969/132‑1 (MQ=255)
                                                  tGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTg                                                                                                                   >  8:294969/1‑132 (MQ=255)
                                                  tGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGc                                                                                                   >  2:21282/1‑148 (MQ=255)
                                                  tGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGc                                                                                                   >  2:14878/1‑148 (MQ=255)
                                                  tGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGc                                                                                                   >  8:306901/1‑148 (MQ=255)
                                                   gACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGt                                                                                                                 <  2:247870/133‑1 (MQ=255)
                                                   gACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGt                                                                                                                 >  1:247870/1‑133 (MQ=255)
                                                    aCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCca                                                                                                 <  1:191859/148‑1 (MQ=255)
                                                     cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAAGGCGCACTTGGTCAACAATCTGAGGCcagt                                                                                               <  8:109418/149‑1 (MQ=255)
                                                     cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCcagt                                                                                               <  1:303478/149‑1 (MQ=255)
                                                      aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCca                                                                                                 <  5:213537/146‑1 (MQ=255)
                                                      aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCca                                                                                                 >  6:213537/1‑146 (MQ=255)
                                                      aaaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACATGGTCAACAATCTGAGGCcagt                                                                                               >  6:321159/1‑148 (MQ=255)
                                                          gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGATCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTa                                                                                          <  2:7484/149‑1 (MQ=255)
                                                          gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACATTATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCcagtcagt                                                                                           >  5:401849/1‑148 (MQ=255)
                                                           cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTaa                                                                                         <  7:257084/149‑1 (MQ=255)
                                                                cttccTGATGAGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCTCTGTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTa                                                                                    <  8:85068/146‑1 (MQ=25)
                                                                cTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACATTATATTGAATTAGCACGATTTCTTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTa                                                                                    <  1:263269/149‑1 (MQ=255)
                                                                    ctgctgCGCTACGCTTCTCCGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGTTTTGGTAGTACGGATAATGCGTTCAAGCCGCATCCGGCGTGTACAACTATCACTTGTTCAAGAATCTGAGGGCTGTCAGTAATGTTACTGc                                                                                <  4:133088/149‑1 (MQ=21)
                                                                     tGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGTCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTgac                                                                               >  6:140018/1‑147 (MQ=255)
                                                                       aTGCGCTACTTTTATAATGGCTACGTGGAGCCCACATTATATTGAATTAGCACTATTTCGTAGGCTTGATAAGGCTTTCACGATGAATCCGGTGTTTTCTTCGCTAAATTTGTATAAATTATCAGGCCAGTCAGTAATGTTACTGCCTg                                                                             <  6:270064/149‑1 (MQ=255)
                                                                          cgcTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTTGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGc                                                                           <  1:201901/148‑1 (MQ=255)
                                                                           gcTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                     >  6:243713/1‑89 (MQ=18)
                                                                           gcTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                     <  5:243713/89‑1 (MQ=18)
                                                                             tACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCtt                                                                        >  2:351257/1‑148 (MQ=255)
                                                                               cGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCACTCAGTAATGTTACTGCCTGGCCTTTaa                                                                     <  2:136194/149‑1 (MQ=255)
                                                                                 cTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTaa                                                                     <  2:48282/147‑1 (MQ=255)
                                                                                 cTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTaa                                                                     >  1:48282/1‑147 (MQ=255)
                                                                                    aTCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACTATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGGGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCaa                                                                 <  5:117230/148‑1 (MQ=255)
                                                                                    aTCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCaa                                                                 <  1:321382/148‑1 (MQ=255)
                                                                                     tCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAAt                                                                >  8:157062/1‑148 (MQ=255)
                                                                                     tCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATg                                                               <  7:157062/149‑1 (MQ=255)
                                                                                       aGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGaa                                                             >  3:239845/1‑149 (MQ=255)
                                                                                       aGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGaa                                                             <  1:21282/149‑1 (MQ=255)
                                                                                       aGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGaa                                                             <  1:14878/149‑1 (MQ=255)
                                                                                        ggCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAAt                                                            >  1:192016/1‑149 (MQ=255)
                                                                                        ggCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAAt                                                            >  8:122275/1‑149 (MQ=255)
                                                                                          ccTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCa                                                          <  7:253057/149‑1 (MQ=255)
                                                                                              cGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGaa                                                       <  6:401849/148‑1 (MQ=255)
                                                                                               gTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAAt                                                                <  3:341748/138‑1 (MQ=255)
                                                                                               gTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAAt                                                                >  4:341748/1‑138 (MQ=255)
                                                                                                 ggAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGc                                                   >  5:398302/1‑149 (MQ=255)
                                                                                                       cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATg                                                                                       <  1:13605/107‑1 (MQ=255)
                                                                                                       cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATg                                                                                       >  2:13605/1‑107 (MQ=255)
                                                                                                       cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTa                                                >  3:305924/1‑146 (MQ=255)
                                                                                                       cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACa                                              <  4:305924/148‑1 (MQ=255)
                                                                                                       cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAg                                             >  6:223436/1‑149 (MQ=255)
                                                                                                           atatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATTAATCAGAAACGCTTACAGCGcc                                         <  4:239845/149‑1 (MQ=255)
                                                                                                           atatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGcc                                         <  5:223436/149‑1 (MQ=255)
                                                                                                            tatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCa                                        <  7:47339/149‑1 (MQ=255)
                                                                                                              tatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCAt                                       >  4:125330/1‑148 (MQ=255)
                                                                                                              tatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATc                                      <  2:192016/149‑1 (MQ=255)
                                                                                                               atTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCAGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATc                                      >  8:75398/1‑148 (MQ=255)
                                                                                                                ttGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACTATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATGAATGAATAAGGAACGCTTACCGCGCCCTCaa                                    >  8:385636/1‑149 (MQ=37)
                                                                                                                ttGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCaa                                    >  2:340151/1‑149 (MQ=255)
                                                                                                                   aaTTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgc                                            <  5:84602/138‑1 (MQ=255)
                                                                                                                   aaTTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgc                                            >  6:84602/1‑138 (MQ=255)
                                                                                                                       aGCACGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTc                              >  6:264578/1‑148 (MQ=255)
                                                                                                                        gCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTc                              >  4:331201/1‑147 (MQ=255)
                                                                                                                        gCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTc                              >  5:294648/1‑147 (MQ=255)
                                                                                                                             aTTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTTGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGAc                       <  7:306901/149‑1 (MQ=255)
                                                                                                                                tCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCa                            >  2:107687/1‑141 (MQ=255)
                                                                                                                                tCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCa                            <  1:107687/141‑1 (MQ=255)
                                                                                                                                      gCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGaaa              <  5:282174/149‑1 (MQ=255)
                                                                                                                                        cGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAaat             >  1:44893/1‑148 (MQ=255)
                                                                                                                                            tAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAg        >  8:248827/1‑149 (MQ=255)
                                                                                                                                            tAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAg        <  7:350260/149‑1 (MQ=255)
                                                                                                                                             aaGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAgc       >  5:314350/1‑149 (MQ=255)
                                                                                                                                                gCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAgccgc    >  6:239020/1‑149 (MQ=255)
                                                                                                                                                  gTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGccgccg  >  1:61420/1‑149 (MQ=255)
                                                                                                                                                    tCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGa                                                        <  4:304622/93‑1 (MQ=255)
                                                                                                                                                    tCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGa                                                        >  3:304622/1‑93 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCG  >  NC_000913/2549400‑2549694

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