Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,763,9920GA39.6% 94.5 / 121.3 111A604T (GCC→ACC) rhsARhs protein with putative toxin 55 domain; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base A (27/17);  total (60/51)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 114 total aligned reads displayed.

TTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAC  >  NC_000913_3_pae_pgi/3763927‑3764141
                                                                 |                                                                                                                                                     
ttGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTaa                                                                   <  1:453637/150‑1 (MQ=9)
 tGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAAc                                                                  <  1:2072/150‑1 (MQ=12)
 tGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAAc                                                                  >  2:1501727/1‑150 (MQ=12)
  gCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAAcg                                                                 <  4:47554/150‑1 (MQ=32)
   ccGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAAcgc                                                                <  2:1651797/150‑1 (MQ=255)
   ccGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGTGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAAcgc                                                                <  1:1201871/150‑1 (MQ=17)
   ccGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAAcgc                                                                >  2:169409/1‑150 (MQ=32)
    cGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTCGCACTCAGGGCGGGCTAAcgcg                                                               >  1:14902/1‑150 (MQ=11)
     gTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACgcgc                                                              >  1:1271934/1‑150 (MQ=32)
     gTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCCAACgcgc                                                              >  4:137322/1‑150 (MQ=17)
      tGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCa                                                             >  2:66819/1‑150 (MQ=12)
      tGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCa                                                             <  4:572672/150‑1 (MQ=32)
          aGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTAt                                                         >  4:567409/1‑150 (MQ=12)
          aGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTAt                                                         >  2:1232524/1‑150 (MQ=32)
          aGACACGCAGGGCCATGAAACGAGGGATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTAt                                                         >  1:822507/1‑150 (MQ=11)
           gACACGCAGGGCCATGAAACGCGGTATGAATACAGCATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATg                                                        <  2:939116/150‑1 (MQ=17)
           gACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATg                                                        >  4:51952/1‑150 (MQ=255)
           gACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATg                                                        <  4:338081/150‑1 (MQ=12)
           gACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATg                                                        <  1:475126/150‑1 (MQ=32)
            acacGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATgg                                                       >  1:942302/1‑150 (MQ=32)
            acacGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTACgg                                                       >  4:272674/1‑150 (MQ=17)
             cacGCAGGGCCATGAAACGCGGTATGAATGCAACATCGCCGGTGACCTGACCGCCGTCAGGGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGa                                                      <  4:944675/150‑1 (MQ=12)
             cacGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGa                                                      >  4:666708/1‑150 (MQ=32)
              acGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGaa                                                     <  4:593779/150‑1 (MQ=255)
              acGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGaa                                                     >  1:996567/1‑150 (MQ=255)
                gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATa                                                   >  1:599806/1‑150 (MQ=12)
                gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATa                                                   >  3:4469/1‑150 (MQ=255)
                 cAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTAt                                                         <  2:704944/143‑1 (MQ=32)
                 cAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTAt                                                         >  1:704944/1‑143 (MQ=32)
                  ggggCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACg                                                 <  2:574062/149‑1 (MQ=255)
                    ggCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGAt                                               <  4:534445/150‑1 (MQ=32)
                     gCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGAtg                                              <  1:1182181/150‑1 (MQ=32)
                      ccATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGAtgc                                             >  3:115464/1‑150 (MQ=32)
                         tGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGAtgctgc                                          >  2:1064741/1‑150 (MQ=12)
                         tGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGAtgctgc                                          >  2:1449061/1‑150 (MQ=32)
                         tGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGAtgctgc                                          <  2:800636/150‑1 (MQ=32)
                         tGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGAtgctgc                                          >  1:393349/1‑150 (MQ=32)
                          gAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGcc                                         <  3:666708/150‑1 (MQ=32)
                           aaaCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCg                                        <  2:1295952/150‑1 (MQ=32)
                           aaaCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCg                                        >  3:283026/1‑150 (MQ=32)
                           aaaCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCg                                        >  1:1045723/1‑150 (MQ=32)
                             aCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGa                                      <  3:272674/150‑1 (MQ=32)
                               gcgGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCCGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACg                                    >  4:141394/1‑150 (MQ=255)
                                cgGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACgg                                   <  4:759363/150‑1 (MQ=12)
                                 ggTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACggg                                  <  1:140048/150‑1 (MQ=255)
                                 ggTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACggg                                  <  4:868774/150‑1 (MQ=32)
                                  gTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGt                                 >  2:1859/1‑150 (MQ=255)
                                  gTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGt                                 <  2:583448/150‑1 (MQ=12)
                                  gTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGt                                 <  4:283026/150‑1 (MQ=32)
                                   tATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTc                                >  4:61886/1‑150 (MQ=255)
                                    aTGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCa                               >  2:399812/1‑150 (MQ=255)
                                     tGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCAt                              >  1:1640069/1‑150 (MQ=32)
                                     tGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCAt                              >  3:1027537/1‑150 (MQ=32)
                                     tGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGACGTAAGTACAAGAAAGAGCGGGATAAAGCGTAGTATGAGATAGGATGATGCCGGACGGGTAAt                              >  1:1421144/1‑150 (MQ=11)
                                      gAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCg                                        >  2:337472/1‑139 (MQ=33)
                                      gAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCg                                        <  1:337472/139‑1 (MQ=33)
                                      gAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATc                             >  2:394848/1‑150 (MQ=32)
                                       aaTACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGCACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGACCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATcc                            >  4:612510/1‑150 (MQ=17)
                                         tACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTc                                >  1:1629637/1‑144 (MQ=32)
                                         tACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTc                                <  2:1629637/144‑1 (MQ=32)
                                         tACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGc                          <  1:1299253/150‑1 (MQ=32)
                                          acaacgccGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGcc                         >  1:1183102/8‑150 (MQ=12)
                                          acaacgCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGcc                         <  2:1335196/150‑1 (MQ=12)
                                           caacgccGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCt                        <  2:822507/144‑1 (MQ=12)
                                            aacgccGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCCCGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTg                       >  3:352362/6‑150 (MQ=255)
                                            aacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTg                       >  3:256810/1‑150 (MQ=32)
                                            aacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTg                       >  3:122551/1‑150 (MQ=32)
                                            aacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTg                       >  4:896483/1‑150 (MQ=32)
                                            aacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTg                       >  4:154604/1‑150 (MQ=32)
                                             acgccGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGa                      >  4:879435/5‑150 (MQ=12)
                                              cgccGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAc                     <  2:157713/147‑1 (MQ=12)
                                              caTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAc                     <  2:1271934/150‑1 (MQ=32)
                                              caTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAc                     <  3:896483/150‑1 (MQ=32)
                                               gccGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGAGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAcc                    >  1:385290/3‑150 (MQ=12)
                                                 cGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAg                  >  2:603239/1‑150 (MQ=32)
                                                  gCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGt                 >  2:1597441/1‑150 (MQ=32)
                                                   ccGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTg                <  4:805005/150‑1 (MQ=32)
                                                     ggTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCt                        >  3:511208/1‑140 (MQ=33)
                                                     ggTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCt                        <  4:511208/140‑1 (MQ=33)
                                                     ggTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGaa              <  1:1380749/150‑1 (MQ=32)
                                                       ggacctgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAGCAGTGaaaa            <  1:1146778/149‑1 (MQ=17)
                                                        gacctgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAAc           <  2:516241/150‑1 (MQ=32)
                                                        gacctgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAAc           <  1:710537/150‑1 (MQ=32)
                                                         acctgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACg          >  3:770155/1‑150 (MQ=32)
                                                         acctgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACg          <  2:1443048/150‑1 (MQ=32)
                                                           ctgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGc        >  1:1336462/1‑150 (MQ=32)
                                                           ctgaccACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGc        >  3:366103/1‑150 (MQ=17)
                                                            tgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCa       <  2:437355/150‑1 (MQ=32)
                                                            tgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCa       <  4:770155/150‑1 (MQ=32)
                                                            tgaccACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCa       >  4:640830/1‑150 (MQ=17)
                                                            tgaccACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCa       <  4:29449/150‑1 (MQ=17)
                                                            tgaccACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCa       <  4:654372/150‑1 (MQ=17)
                                                             gaccACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAg      >  2:1321445/1‑150 (MQ=17)
                                                              accaccGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGc     >  2:681037/1‑150 (MQ=17)
                                                              accGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGc     <  3:1034572/150‑1 (MQ=32)
                                                              accGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGc     <  1:316748/150‑1 (MQ=32)
                                                              accGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGc     <  1:1232524/150‑1 (MQ=32)
                                                              accGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGCACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGc     <  3:397272/150‑1 (MQ=17)
                                                                cgccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCca   >  2:684571/1‑150 (MQ=32)
                                                                 gccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCaattaa  <  4:596126/150‑6 (MQ=17)
                                                                 |                                                                                                                                                     
TTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAC  >  NC_000913_3_pae_pgi/3763927‑3764141

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: ‑