Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,295,3200CT28.2% 133.5 / 59.4 103intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base T (11/18);  total (49/54)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 103 total aligned reads displayed.

TAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913_3_pae_pgi/4295178‑4295466
                                                                                                                                              |                                                                                                                                                  
tAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                              <  1:377976/149‑1 (MQ=32)
tAAGGCGTGTATTCCGCACACGAGATGAACGTCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                              <  2:252377/149‑1 (MQ=11)
 aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTGCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGCTCGGATAAGGCGTTTACGCTGCATCCGACATCAATGCCTGac                                                                                                                                             >  1:104127/1‑148 (MQ=255)
      gTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTGAGCGTGATTTTGTACGTCGGATAAGGCGTTTACGCCGCATCCGACTTCAATCCCTGATGcacc                                                                                                                                        >  1:293103/1‑146 (MQ=11)
      gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACTGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCggc                                                                                                                                        >  2:516078/1‑147 (MQ=17)
         tATTCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCt                                                                                                                                     <  2:1011121/149‑1 (MQ=17)
         tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGGGTGATTTCGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGAAGAGACGCt                                                                                                                                     >  1:958993/1‑149 (MQ=255)
          aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                                                                                                    >  1:632846/1‑149 (MQ=17)
            gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                                  <  1:1044952/149‑1 (MQ=34)
              cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATCCGACCCTTGccc                                                                                                                                >  2:1023685/1‑148 (MQ=255)
               gaaTCCGACATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                               <  1:608477/147‑1 (MQ=37)
               cgaTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                               <  1:160713/147‑1 (MQ=25)
                 aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTTCGGCAGACAGCGCAATAGCGTGATTTAGCGTGAGTTTGTGGGTCGGATAAGGCGTTTAGGACGCTTCCAACATCAATGCCTGATGCGACGCTTGCGGCGt                                                                                                                             >  2:287654/1‑149 (MQ=11)
                 aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                             <  1:1025866/149‑1 (MQ=25)
                 aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                             <  1:881127/149‑1 (MQ=21)
                  tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                            >  1:80074/1‑149 (MQ=38)
                    cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGGCTTAGCGTCAATTTCTAGGCCGGACAAGGCGATTACACCACAATCGACATCAACGCCTCATGCGCCGCCTGCCGCGTcct                                                                                                                          >  2:728723/1‑147 (MQ=11)
                    cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                          >  1:663448/1‑149 (MQ=25)
                    cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGACCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTTTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCTACGCTTGCCGCGGCtt                                                                                                                          >  1:640836/1‑149 (MQ=11)
                     gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                         >  2:41428/1‑149 (MQ=37)
                       cATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                       >  1:763880/1‑149 (MQ=25)
                         tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                     <  2:80074/149‑1 (MQ=39)
                         gcAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                     <  1:340762/148‑1 (MQ=25)
                         ccAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACCCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                     <  1:287654/148‑1 (MQ=255)
                          caccGCCTGATGCGACGCTTAACCCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                    <  1:395197/146‑1 (MQ=255)
                           aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                   >  1:558767/1‑149 (MQ=21)
                           aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                   >  1:365991/1‑149 (MQ=37)
                            aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGa                                                                                                                                         >  1:52847/1‑126 (MQ=17)
                            aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGa                                                                                                                                         <  2:52847/126‑1 (MQ=18)
                              gCCTGATGCGACGCTGAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                <  2:360477/149‑1 (MQ=255)
                               ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                               >  1:796002/1‑148 (MQ=18)
                                tctaTGCAACGCTTAACGCGTCTTATCAGGCATAAGCCAGACAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTCATGCGACGCTTGCCGCGTCTTATCAGGCCTACg                                                                                                              <  1:71804/146‑1 (MQ=11)
                                cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                              <  1:620063/149‑3 (MQ=17)
                                cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                              <  2:365991/149‑3 (MQ=18)
                                cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                              <  1:759404/149‑3 (MQ=18)
                                 tGATGGGACGCGTAACGCGTCTTATCAGGCCGACGCCAGGCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                             <  2:23607/149‑1 (MQ=11)
                                 ggATGGGACGCCTAACGCGTTTTATCAGGCTTACTACTGACAGGGCAATCGCCTAAATTTGCGTGATTTTGTAGGTCGTTTACGTCCTTTACGCCGTATCGGTCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                             <  2:1038668/148‑4 (MQ=2)
                                     gCGACGCTTGACGCGGGTTACAAGGCCTCCGCCAGAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatcttaac                                                                                                         <  1:225478/149‑8 (MQ=11)
                                          gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCCATACCCTGATTTAGCGTAATTATGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCACGCCTACGCCAGACAgcg                                                                                                    >  2:606352/1‑149 (MQ=11)
                                          gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                                  >  2:80039/1‑119 (MQ=35)
                                          gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                                  <  1:80039/119‑1 (MQ=35)
                                                cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                  <  1:311423/97‑1 (MQ=21)
                                                cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                                                                                                  >  2:311423/1‑97 (MQ=33)
                                                 gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                             >  2:813243/1‑149 (MQ=37)
                                                        aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGCCGGGTAAGGGGTTTACGCCGCATCCGGCATGAACGCCTGATGCGACGCTGGCCGCGTATTATCAGGCCTACGCCAGACAGCGAAACAGCCTGAttt                                                                                      >  2:87946/1‑149 (MQ=11)
                                                           aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAGAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                   <  1:87946/149‑1 (MQ=38)
                                                            ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                  <  1:41428/149‑1 (MQ=38)
                                                                tACGCCAGACAGCGCAATAGCCTGAATTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCATGCCTACCCCAACCAGCGCAGTAGCCTGATTTATCgagaa                                                                              >  2:225478/1‑148 (MQ=11)
                                                                  cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                             <  2:723173/100‑1 (MQ=255)
                                                                  cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                             >  1:723173/1‑100 (MQ=255)
                                                                   gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                    <  2:890345/108‑1 (MQ=12)
                                                                     cAGACAGTGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAcc                                                                                                                                 >  2:71802/1‑93 (MQ=12)
                                                                          aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                      >  1:86185/1‑99 (MQ=255)
                                                                          aGCGCAATAGCCTGATTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                      <  2:86185/99‑1 (MQ=255)
                                                                               aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTACGCCTGATGTCGGAGGTCGGAtc                                                               >  2:558227/1‑148 (MQ=11)
                                                                                 tggtcTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAg                                                             <  2:705725/145‑1 (MQ=37)
                                                                                     cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                            <  1:1052314/82‑1 (MQ=255)
                                                                                     cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                            >  2:1052314/1‑82 (MQ=255)
                                                                                       gATTTAGCGTGTTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                           <  1:957110/97‑1 (MQ=255)
                                                                                       gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                           >  2:957110/1‑97 (MQ=255)
                                                                                        aTTGAGCGTAATTTTGTAGTCGGGAGAGGGCGTCTACGCCGCATCCGGGATCAACGCCTGATGCGACGCTTGCCTCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTa                                                      <  2:37562/149‑1 (MQ=11)
                                                                                         tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTc                                                                    >  2:804956/1‑134 (MQ=39)
                                                                                           tAGCGTGATTTTGTGGGGCTGTTACGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGc                                                   <  2:958993/149‑1 (MQ=11)
                                                                                                tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATCAGGCGTTTATGCCGCAt                                              >  1:972917/1‑149 (MQ=255)
                                                                                                   ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                     <  1:462539/75‑1 (MQ=255)
                                                                                                   ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                            <  2:905299/68‑1 (MQ=255)
                                                                                                   ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                            >  1:905299/1‑68 (MQ=255)
                                                                                                   ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                     >  2:462539/1‑75 (MQ=255)
                                                                                                          ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCaa                                    >  1:980056/1‑149 (MQ=39)
                                                                                                           gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                         <  1:775691/63‑1 (MQ=255)
                                                                                                           gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                         >  2:775691/1‑63 (MQ=255)
                                                                                                           gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAttt                                                                            <  1:862598/108‑1 (MQ=255)
                                                                                                              ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                             <  1:888639/72‑1 (MQ=255)
                                                                                                              ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                             >  2:888639/1‑72 (MQ=255)
                                                                                                               gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCt                               >  2:407258/1‑149 (MQ=38)
                                                                                                                 cagggCGTTTCCGCCGCACCCGGTATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGa                             <  2:692243/146‑1 (MQ=11)
                                                                                                                  aaGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                                                         <  1:744753/40‑1 (MQ=40)
                                                                                                                   aGGCGTTTCCGCCGCACCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATg                           <  1:285478/149‑1 (MQ=255)
                                                                                                                   aGGCGTTTACGCCCCCTCCGCCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATg                           <  1:287208/149‑1 (MQ=255)
                                                                                                                    ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGAGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGAATAAGGCGTTGATGCCCCATCCGACATCAACGCCTGATGc                          >  2:550584/1‑149 (MQ=255)
                                                                                                                     gCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                               <  2:729915/63‑2 (MQ=255)
                                                                                                                           aCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTg                                                                                >  2:1030874/1‑88 (MQ=255)
                                                                                                                             gccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATg                                                    >  1:317110/1‑114 (MQ=38)
                                                                                                                              ccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACa                                        <  1:25176/125‑1 (MQ=37)
                                                                                                                              ccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACa                                        >  2:25176/1‑125 (MQ=37)
                                                                                                                               cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                             >  1:772623/1‑55 (MQ=255)
                                                                                                                               cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                             <  2:772623/55‑1 (MQ=255)
                                                                                                                               cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGa                                                                               >  2:396468/1‑85 (MQ=40)
                                                                                                                               cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcg               >  1:922023/1‑149 (MQ=35)
                                                                                                                                gcACCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc              <  1:990738/149‑1 (MQ=21)
                                                                                                                                 cacccGACAGCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgcg             <  1:229242/146‑1 (MQ=21)
                                                                                                                                  tggagACATCGGCGCCTGATGGGACGCTTGGCGTGTATTATCAGGGAGACACAAGACAGAGCAATAGGCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGGGTTTATGCCGCATCCGATATCAACGCCTGATGCGACGCTTAACGCGt            <  2:252503/145‑1 (MQ=11)
                                                                                                                                  aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGt            <  1:813243/149‑1 (MQ=34)
                                                                                                                                    ccGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCt          >  1:139546/1‑149 (MQ=35)
                                                                                                                                     cGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                 <  1:686856/61‑1 (MQ=255)
                                                                                                                                     cGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                 >  2:686856/1‑61 (MQ=255)
                                                                                                                                         ctCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGAGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCa     >  2:820036/2‑149 (MQ=18)
                                                                                                                                            aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGc                                                           <  1:1054061/92‑1 (MQ=33)
                                                                                                                                            aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAAGCCTGATGCGACGCTTACCGCGTCTTATCAGGc  >  1:1033710/1‑149 (MQ=17)
                                                                                                                                            aaCGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc  <  1:1047324/149‑1 (MQ=17)
                                                                                                                                              |                                                                                                                                                  
TAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913_3_pae_pgi/4295178‑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: