Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,727,7680GT36.4% 130.2 / 120.3 108intergenic (+12/‑69)nemA/gloAchromate reductase, quinone reductase, FMN‑linked; N‑Ethylmaleimide reductase; old yellow enzyme/glyoxalase I, Ni‑dependent
Reads supporting (aligned to +/- strand):  ref base G (38/30);  new base T (23/16);  total (62/46)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 109 total aligned reads displayed.

GTGGCATTTGGTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATG  >  NC_000913_3_pae_pgi/1727619‑1727914
                                                                                                                                                     |                                                                                                                                                  
gTGGCATTTGGTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCg                                                                                                                                                    >  4:552882/1‑150 (MQ=255)
  ggCATTTGGTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAg                                                                                                                                                  >  2:929829/1‑150 (MQ=255)
         ggTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTa                                                                                                                                           >  4:211351/1‑150 (MQ=255)
          gTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTaa                                                                                                                                          >  2:207468/1‑150 (MQ=255)
            cGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAg                                                                                                                                        >  4:512573/1‑150 (MQ=255)
                  tGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGccgccg                                                                                                                                  >  3:50545/1‑150 (MQ=255)
                   ggATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCcgccgc                                                                                                                                 <  2:1069422/150‑1 (MQ=255)
                    gATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCt                                                                                                                                <  3:306499/150‑1 (MQ=255)
                       tGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCtata                                                                                                                             <  3:211351/150‑1 (MQ=255)
                       tGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCtata                                                                                                                             <  1:327686/150‑1 (MQ=255)
                        gCGAGCCCGGATCTGGTCGCCCGCTTGAAGCGGAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATAc                                                                                                                            <  2:1125986/150‑1 (MQ=255)
                          gAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGATCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAgcggcggcttt                                                                                                                                        <  2:735612/136‑6 (MQ=255)
                          gAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAgcggcggcttt                                                                                                                                        >  1:735612/1‑131 (MQ=255)
                          gAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTa                                                                                                                          >  1:929463/1‑150 (MQ=255)
                          gAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTa                                                                                                                          >  1:929279/1‑150 (MQ=255)
                           aaCCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTaa                                                                                                                         <  4:140442/150‑1 (MQ=255)
                             cccGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTaaaa                                                                                                                       <  1:100485/150‑1 (MQ=255)
                             cccGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTaaaa                                                                                                                       <  1:103273/150‑1 (MQ=255)
                             cccGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTaaaa                                                                                                                       <  3:147267/150‑1 (MQ=255)
                              ccGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAaac                                                                                                                      >  3:88892/1‑150 (MQ=255)
                              ccGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTACTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAaac                                                                                                                      >  1:270412/1‑150 (MQ=255)
                               cGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAaaca                                                                                                                     >  2:724465/1‑150 (MQ=255)
                                 gATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAaacaac                                                                                                                   <  3:704571/150‑1 (MQ=255)
                                      ggTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACAtttt                                                                                                              <  1:832985/150‑1 (MQ=255)
                                      ggTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACAtttt                                                                                                              <  2:783676/150‑1 (MQ=255)
                                            ccGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCt                                                                                                        <  2:869159/150‑1 (MQ=255)
                                              gCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGt                                                                                                      >  1:942625/1‑150 (MQ=255)
                                              gCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGt                                                                                                      >  1:982539/1‑150 (MQ=255)
                                                ttGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTa                                                                                                    <  1:603435/150‑1 (MQ=255)
                                                ttGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTa                                                                                                    <  1:75155/150‑1 (MQ=255)
                                                   cAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGcc                                                                                                 >  1:208548/1‑150 (MQ=255)
                                                    aGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTg                                                                                                             <  1:652047/137‑1 (MQ=255)
                                                    aGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTg                                                                                                             >  2:652047/1‑137 (MQ=255)
                                                    aGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCa                                                                                                >  1:1043526/1‑150 (MQ=255)
                                                       gcAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCAttt                                                                                             >  3:376979/1‑150 (MQ=255)
                                                        cAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCAtttt                                                                                            >  2:712644/1‑150 (MQ=255)
                                                          aaGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGa                                                                                          <  3:735336/150‑1 (MQ=255)
                                                           aGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAg                                                                                         <  1:747/150‑1 (MQ=255)
                                                           aGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAg                                                                                         <  1:97490/150‑1 (MQ=255)
                                                             cTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGa                                                                                       >  4:553427/1‑150 (MQ=255)
                                                             cTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGa                                                                                       >  1:12992/1‑150 (MQ=255)
                                                              tGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGAt                                                                                      >  1:643860/1‑150 (MQ=255)
                                                                  cTTAACCCACAGCGTGCCGAAAGTTTCTACGGAGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATaaaa                                                                                  >  4:302739/1‑150 (MQ=255)
                                                                     aaCCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGa                                                                               <  2:211527/150‑1 (MQ=255)
                                                                       cccACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATg                                                                             >  1:380902/1‑150 (MQ=255)
                                                                         cacaGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCg                                                                           >  2:572193/1‑150 (MQ=255)
                                                                         cacaGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCg                                                                           >  1:333578/1‑150 (MQ=255)
                                                                         cacaGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCg                                                                           >  1:1097123/1‑150 (MQ=255)
                                                                              cGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGtcttc                                                                      >  2:242231/1‑150 (MQ=255)
                                                                              cGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGtcttc                                                                      <  2:352461/150‑1 (MQ=255)
                                                                              cGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGtcttc                                                                      >  1:472468/1‑150 (MQ=255)
                                                                                  ccGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCa                                                                  >  2:124432/1‑150 (MQ=255)
                                                                                  ccGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCa                                                                  <  1:1090075/150‑1 (MQ=255)
                                                                                     aaaGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATAc                                                               >  2:929406/1‑150 (MQ=255)
                                                                                         tttCTACGGTGGCGGCGCGGGAGGCTATACCGAGTACCCGCGGGTGTAATCCAAGATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGACGCTAAAAAGATGCCTCTCCTCCATACCAtg                                                           >  2:112216/1‑150 (MQ=255)
                                                                                           tCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCAtgct                                                         >  1:1116044/1‑150 (MQ=255)
                                                                                           tCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCAtgct                                                         >  3:595413/1‑150 (MQ=255)
                                                                                             tACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCAtgctgc                                                       <  1:381957/150‑1 (MQ=255)
                                                                                                ggTGGCGGCGCTAACAGTTATACAGAATACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGcgcg                                                    >  4:426325/1‑150 (MQ=255)
                                                                                                 gTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGt                                                   >  1:966162/1‑150 (MQ=255)
                                                                                                   ggcggcGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTg                                                 <  2:7990/150‑1 (MQ=255)
                                                                                                    gcggcgCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCCGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTgg                                                <  3:209401/150‑1 (MQ=255)
                                                                                                      ggcgCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCg                                              >  2:1019064/1‑150 (MQ=255)
                                                                                                       gcgcGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGa                                             <  1:724465/150‑1 (MQ=255)
                                                                                                        ggcgGAAGGCAATGCCGATTACCCGACGTGGTAATCCAACATTGCGAGCGGCGTGATGCCGCCGCTATACTAAAACAACATGTTGAATGTTTTAGCCATTTTGAGGATAAACAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGAt                                            <  2:489921/149‑1 (MQ=255)
                                                                                                          cgGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGAttt                                          <  2:12992/150‑1 (MQ=255)
                                                                                                              aGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCaa                                      <  1:712644/150‑1 (MQ=255)
                                                                                                                  tataCCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCt                                  >  4:347147/1‑150 (MQ=255)
                                                                                                                     aCCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCa                               >  1:645695/1‑150 (MQ=255)
                                                                                                                      ccGATTACCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCAt                              <  1:1019064/150‑1 (MQ=255)
                                                                                                                        gATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCg                            >  2:659252/1‑150 (MQ=255)
                                                                                                                            aCCCGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGAttt                        <  1:242231/150‑1 (MQ=255)
                                                                                                                            aCCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGAttt                        <  1:442434/150‑1 (MQ=255)
                                                                                                                             cccGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAGGTATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGAtttt                       >  3:389717/1‑150 (MQ=255)
                                                                                                                             cccGACGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGAtttt                       >  1:231604/1‑150 (MQ=255)
                                                                                                                             cccGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGCTCAAACGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGAtttt                       >  1:1181084/1‑150 (MQ=255)
                                                                                                                             cccGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGa                                                                               <  1:463855/94‑1 (MQ=255)
                                                                                                                             cccGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGa                                                                               >  2:463855/1‑94 (MQ=255)
                                                                                                                             cccGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGAtttt                       <  2:333578/150‑1 (MQ=255)
                                                                                                                             cccGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGAtttt                       >  4:399784/1‑150 (MQ=255)
                                                                                                                              ccGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGAttttt                      >  1:687620/1‑150 (MQ=255)
                                                                                                                               cGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGc                                                                                                  <  4:431096/73‑1 (MQ=255)
                                                                                                                               cGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAg                                                                                <  3:302175/91‑1 (MQ=255)
                                                                                                                               cGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAg                                                                                >  4:302175/1‑91 (MQ=255)
                                                                                                                                  cGTTGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATAc                  <  3:553427/150‑1 (MQ=255)
                                                                                                                                     tGTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCaa               >  3:505914/1‑150 (MQ=255)
                                                                                                                                     tGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCaa               >  2:1061682/1‑150 (MQ=255)
                                                                                                                                      gTAATCCAACATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCaaa              <  3:512573/150‑1 (MQ=255)
                                                                                                                                         aTCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTCGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTg           >  3:638128/1‑150 (MQ=255)
                                                                                                                                         aTCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTCGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTg           >  3:638165/1‑150 (MQ=255)
                                                                                                                                           ccAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGa                                             <  1:337017/114‑1 (MQ=255)
                                                                                                                                           ccAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGa                                             >  2:337017/1‑114 (MQ=255)
                                                                                                                                             gaCATTGCAATGGGCGTAAAGCCGCCGCTATACTACCACAACATTTAGAATCTGTTAGCCCTTTTGAGGCTAAACAGATGCTACTTCTTCATACCATGCTGCGCGTTGGCGATTAGCAACGCTCCATCGATTTTTATACCAAAGTGCTgg       >  4:40428/2‑150 (MQ=255)
                                                                                                                                             aaCATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTgcgc                                                     <  2:141671/104‑1 (MQ=255)
                                                                                                                                             aaCATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTgcgc                                                     >  1:141671/1‑104 (MQ=255)
                                                                                                                                             aaCATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTgg       >  1:450810/1‑150 (MQ=255)
                                                                                                                                             aaCATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTgg       <  2:204856/150‑1 (MQ=255)
                                                                                                                                             aaCATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTgg       >  2:827777/1‑150 (MQ=255)
                                                                                                                                              aCATTGCTAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTggg      <  2:232915/150‑1 (MQ=255)
                                                                                                                                                  tGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATg  >  3:233944/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTGGCATTTGGTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATG  >  NC_000913_3_pae_pgi/1727619‑1727914

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