Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,295,3200CT22.4% 169.2 / 51.3 98intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (44/32);  new base T (10/12);  total (54/44)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913_3_pae_pgi/4295174‑4295466
                                                                                                                                                  |                                                                                                                                                   
cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                   <  2:442576/149‑1 (MQ=17)
cGGATAAGGCCTTAATGCCGCATCCTACATCATCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGTTTTCGCCGCATCCGACATCAACGc                                                                                                                                                   <  2:853677/149‑1 (MQ=11)
  tgtAAGGCGTTTATGCCCCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                                 <  2:76019/147‑1 (MQ=12)
   aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTg                                                                                                                                                >  1:442576/1‑149 (MQ=32)
       ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                            >  2:373575/1‑149 (MQ=17)
       ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                                                                                                            >  1:577111/1‑149 (MQ=32)
       ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACAGCAACGCCTGATGc                                                                                                                                            >  1:708307/1‑149 (MQ=17)
                gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTTCGTGGTGGAACCCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGc                                                                                                                                  <  1:472130/149‑1 (MQ=11)
                gCCCCATCCCACATCAACGCCTGATGCGACGCTTAACGCTTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGc                                                                                                                                  <  1:828308/149‑1 (MQ=255)
                 ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGcc                                                                                                                                 <  1:474267/149‑1 (MQ=21)
                      tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTc                                                                                                                            >  1:604566/1‑149 (MQ=25)
                        cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCtt                                                                                                                          >  1:650846/1‑149 (MQ=38)
                         gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTa                                                                                                                         >  2:668314/1‑149 (MQ=21)
                         gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑GGCCGCGTCTTa                                                                                                                         >  1:143132/1‑149 (MQ=37)
                          aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTAt                                                                                                                        <  1:236803/149‑1 (MQ=21)
                           cATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATc                                                                                                                       <  2:929196/149‑1 (MQ=38)
                            aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTATcc                                                                                                                      >  1:516873/1‑148 (MQ=25)
                             tCAATGACTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTATCAg                                                                                                                     <  1:48951/149‑1 (MQ=255)
                              cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCTATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTATCAgg                                                                                                                    >  1:642548/1‑149 (MQ=21)
                              cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTATCAgg                                                                                                                    <  1:324005/149‑1 (MQ=21)
                              cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTATCAgg                                                                                                                    <  1:111950/149‑1 (MQ=21)
                               aaggtcTGATTCGACGCTTGAGGTGTCTTATCAGGCCTAGGCCCGTCAGCGCAATAGCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTGCGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGc                                                                                                                   <  2:869368/144‑1 (MQ=11)
                               aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGc                                                                                                                   >  1:311479/1‑149 (MQ=21)
                                aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGcc                                                                                                                  >  1:927265/1‑149 (MQ=21)
                                aCGCCTGATGAGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGcg                                                                                                                  >  1:944896/1‑148 (MQ=37)
                                 cGCCTGATTTGTCGCTAAAGGTGTCTTATCAGGCCTACGCCAGACAGCGCAAAAGCCTCATTTCGTGTAATTTTGTAGGTCGACTCAGGCGTTTACGCCGCATCCGCCACCACTGCCTGCTGCGATGCT‑TGCCGCGTTTTATCAGTCCt                                                                                                                 <  2:577111/149‑1 (MQ=2)
                                  gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTa                                                                                                                >  1:671389/1‑149 (MQ=25)
                                  gCCAGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCTCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTa                                                                                                                >  2:608359/1‑149 (MQ=255)
                                     tGATGCGGCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTatct                                                                                                             >  2:281109/1‑146 (MQ=12)
                                           gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACa                                                                                                       >  1:707298/1‑149 (MQ=39)
                                              gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCCGCCAgcg                                                                                                    >  1:811177/1‑149 (MQ=255)
                                                ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCtt                                                                                                                          <  2:216445/125‑1 (MQ=40)
                                                ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCa                                                                                                  <  2:707298/149‑1 (MQ=39)
                                                ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCAGCGTCtt                                                                                                                          >  1:216445/1‑125 (MQ=38)
                                                 tttcgcgTCTTTTCGGGTCTAAGCCTGACAGCGCAATTGCCTGATTCAGCGTGATTTTGTAGGGTGGGTAAGGCGTTTACGCCGCGTCCGACATCAGCGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                 <  2:664563/146‑1 (MQ=11)
                                                  aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCTAGACAGCGCAAt                                                                                                >  2:900149/1‑149 (MQ=37)
                                                    cgcgTCTTATCAGGCCTACCCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATACGACGCT‑TGCCGCGTCTTCTCAGGCCTACCCCAACCAGCGCAATac                                                                                              >  2:473148/1‑148 (MQ=11)
                                                     gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATg                                                                                                                                             >  2:110507/1‑102 (MQ=32)
                                                      cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCATACAGCGCAATATcc                                                                                            >  1:77688/1‑149 (MQ=255)
                                                      cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                            >  1:724140/1‑149 (MQ=38)
                                                             tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                                                          <  2:321799/97‑1 (MQ=33)
                                                               aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTCTGCC‑CGTCCTATGAAGCCTACGCCAGGCCGCGCCGGAGGCTGAtttcac                                                                                   >  1:689587/1‑146 (MQ=1)
                                                                ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGt                                                                                                                             >  2:197684/1‑106 (MQ=255)
                                                                ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGt                                                                                                                             <  1:197684/106‑1 (MQ=255)
                                                                 gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCTGCGTCTTATc                                                                                                                       <  1:944159/111‑1 (MQ=255)
                                                                 gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATc                                                                                                                       >  2:944159/1‑111 (MQ=40)
                                                                    tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAt                                                                              <  1:385536/149‑1 (MQ=38)
                                                                    tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TACCGCGTCTTATCAGGCCTACGCTAGACAGTGCAATAGCCTGGTGTAGCGTgct                                                                              >  2:687488/1‑147 (MQ=11)
                                                                       gtcAGAGAGCGCAATAGCCTGATTGAGCGTGATTTTGTAGGTGGGATAGGGCGTTTCCGCCGGGTCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAtttt                                                                           <  2:782709/147‑1 (MQ=11)
                                                                        ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGcc                                                                                                                  >  1:566433/1‑109 (MQ=255)
                                                                        ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                     <  2:36669/138‑1 (MQ=38)
                                                                        ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                     >  1:36669/1‑138 (MQ=37)
                                                                         cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCa                                                                                                           >  1:696232/1‑115 (MQ=40)
                                                                          aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCTTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                        <  1:407578/149‑1 (MQ=37)
                                                                               gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATGAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCg                                                                   <  1:190308/149‑1 (MQ=38)
                                                                                        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATGGCCTGATTTAGCGTGAACTGGTAGGACGGCTAAGGCg                                                          >  2:307907/1‑149 (MQ=11)
                                                                                        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                    >  2:685082/1‑139 (MQ=39)
                                                                                        ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                    <  1:685082/139‑1 (MQ=40)
                                                                                        ccTGATTTAGCGTGATTTTGGAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCg                                                          <  1:686757/149‑1 (MQ=255)
                                                                                            aTTTAGCGTGATTTTGTAGGTCGGATAAGGCTTTCGTTCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCGTTAGCGTCATTTTGTAGGTCGCATAAGCCGTTTa                                                      >  2:790561/1‑149 (MQ=11)
                                                                                            aTTGAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTa                                                      <  1:823682/149‑1 (MQ=38)
                                                                                                  cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAtt                                                                                       <  1:671324/110‑1 (MQ=255)
                                                                                                  cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAtt                                                                                       >  2:671324/1‑110 (MQ=255)
                                                                                                   gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGc                                                                                                                                  <  2:845024/66‑1 (MQ=40)
                                                                                                   gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGc                                                                                                                                  >  1:845024/1‑66 (MQ=40)
                                                                                                    tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTat                                                                                                               <  2:504505/84‑2 (MQ=255)
                                                                                                    tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTAc                                                                                                               >  1:504505/1‑84 (MQ=255)
                                                                                                    tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCg                                                          >  1:492991/1‑137 (MQ=40)
                                                                                                    tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCg                                                          <  2:492991/137‑1 (MQ=39)
                                                                                                      aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCa                                                                                                           >  1:835021/1‑86 (MQ=255)
                                                                                                         ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCATACACCGCAATAGCCTGCTTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCgcc                                         <  1:672486/149‑3 (MQ=255)
                                                                                                               ttCGGATAAGGCGTTTACCCCGCATCCGACATGAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAc                                   <  1:44221/148‑1 (MQ=255)
                                                                                                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCT‑TGcc                                                                                                                                 <  1:613647/55‑1 (MQ=255)
                                                                                                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGcc                                                                                                                                 >  2:613647/1‑55 (MQ=255)
                                                                                                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATg                                                    <  1:908449/132‑1 (MQ=40)
                                                                                                               gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATg                                                    >  2:908449/1‑132 (MQ=40)
                                                                                                                tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCTGCCGCATCCGACCTCaagg                                  >  1:856051/1‑147 (MQ=255)
                                                                                                                tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGcc                                                  >  2:900744/1‑133 (MQ=37)
                                                                                                                tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGcc                                                  <  1:900744/133‑1 (MQ=38)
                                                                                                                tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACg                                  >  2:886482/1‑149 (MQ=38)
                                                                                                                  ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGcc                                >  1:728565/1‑149 (MQ=39)
                                                                                                                     tAAGGCGTGTACGCCGCATAAGGTATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGa                             <  2:711414/149‑1 (MQ=11)
                                                                                                                           gTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAc                       >  1:582428/1‑149 (MQ=38)
                                                                                                                            tttACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGCCATCAACCCCTGATGCgccc                      >  2:207544/1‑146 (MQ=255)
                                                                                                                            gttctgccgcACTCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCGGGCCTACGCCAGACAGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACg                      <  2:143132/144‑1 (MQ=11)
                                                                                                                             ttACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGc                     <  2:650846/149‑1 (MQ=37)
                                                                                                                              tACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCt                    >  1:582810/1‑149 (MQ=38)
                                                                                                                                 tcagcAGCCGATATCTACGCCTGATGCGACGCT‑TGCGGCGTCTTATCAGGCCTACGCCAGCCAGCGCTACAGCCCGATCTTGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAGCCGACATCAACGCCTGATGCGACGCTTaa                 <  2:469554/146‑1 (MQ=11)
                                                                                                                                  ccccacccGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAc                <  1:214188/143‑1 (MQ=21)
                                                                                                                                   cgcATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCGTTTATGCCGCATCCGACATAAACGCCTGATGCGCCGCTTATcg               >  2:622148/1‑149 (MQ=255)
                                                                                                                                   cgcATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATCGGACGCTTCAcg               >  2:317298/1‑149 (MQ=255)
                                                                                                                                         cGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCtt         >  2:472130/1‑149 (MQ=35)
                                                                                                                                             aTCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTCACGCGTCTTATCa     >  2:324005/1‑149 (MQ=21)
                                                                                                                                             aTCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCa     <  1:307907/149‑1 (MQ=34)
                                                                                                                                              tCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAg    <  1:751450/149‑1 (MQ=33)
                                                                                                                                                agcGCCTGATGCGATGCT‑TGCCGTGTCTTATCAGGCCTAGGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCCCCCCCGAACTCAAACCCAGACGCAACGCTTAACGCGTCTTATCAGGc  <  2:480144/147‑1 (MQ=11)
                                                                                                                                                aaCGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGCCACTGCAATAGCGTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGc  >  1:189854/1‑149 (MQ=11)
                                                                                                                                                 aCGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATcc                                            <  1:617613/106‑1 (MQ=31)
                                                                                                                                                 aCGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATcc                                            >  2:617613/1‑106 (MQ=31)
                                                                                                                                                  |                                                                                                                                                   
CGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCT‑TGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913_3_pae_pgi/4295174‑4295466

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