Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi526,7920AC24.4% 147.3 / 45.4 82S1178R (AGC→CGC) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (25/37);  new base C (8/12);  total (33/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGC  >  NC_000913_3_pae_pgi/526647‑526939
                                                                                                                                                  |                                                                                                                                                   
gCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcgg                                                                                                                                                 >  1:927515/1‑149 (MQ=255)
gCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgcgc                                                                                                                                                 >  1:59882/1‑150 (MQ=255)
gCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGCCCAGACAGGTGGAGCCGGAATAAACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgcgc                                                                                                                                                 >  2:779011/1‑150 (MQ=255)
 cTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTTCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                                <  2:502817/150‑3 (MQ=255)
 cTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                                <  1:658400/150‑3 (MQ=255)
 cTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                                <  2:927515/150‑3 (MQ=255)
   tGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGCACATACGGCGTGGCGCGGGGa                                                                                                                                              >  1:617241/1‑150 (MQ=11)
    gcgcAGGGCGGGCTGACGGTGGAGCAACTGGCCAGGCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGaa                                                                                                                                             <  2:960314/150‑1 (MQ=255)
      gcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACCCCACACGTGTGATATCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGTAATACGGCGTGGCGCGGGGAgta                                                                                                                                           <  1:478479/150‑1 (MQ=11)
        agggcgggcggACGGTGGAGCAACGGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatg                                                                                                                                         <  2:648075/147‑1 (MQ=255)
         gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatga                                                                                                                                        <  1:809602/150‑1 (MQ=255)
         gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatga                                                                                                                                        >  2:278948/1‑150 (MQ=255)
           gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgatg                                                                                                                                      <  2:804686/150‑1 (MQ=255)
            cGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGGCCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAtgatga                                                                                                                                     <  1:241770/150‑1 (MQ=11)
              ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                                                                   <  2:142158/150‑1 (MQ=255)
                cTGACGGGCGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                 <  2:643015/150‑1 (MQ=255)
                 tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                                <  2:930831/150‑1 (MQ=255)
                 tGACGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg                                                                                                                                <  1:662839/150‑1 (MQ=255)
                 tGACGGGGGAGCAACGGGCCAGGCAGGTGGAGCCGGAATGCACACCGGCGCGAAAAGCTCATTTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg                                                                                                                                <  2:915036/150‑1 (MQ=14)
                 ggACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg                                                                                                                                <  1:439845/149‑1 (MQ=255)
                 ggACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg                                                                                                                                <  1:317926/149‑1 (MQ=255)
                  gACGGTGGAGCAACGGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                               <  2:540705/150‑1 (MQ=11)
                  gACGGTGCAGCAACTGGCCAGCCAGGTGGAGCCGGAACACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                               <  1:207090/150‑1 (MQ=11)
                     ggTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCATTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                 >  2:762985/1‑145 (MQ=255)
                     ggTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCATTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                 <  1:762985/145‑1 (MQ=255)
                     ggTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCaa                                                                                                                            >  1:555380/1‑150 (MQ=255)
                     ggTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGACGTGGCGCGGGGAGTATGATGAATGGGGCaa                                                                                                                            >  1:359528/1‑150 (MQ=255)
                      gTGGAGCAACGGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc                                                                                                                           <  1:844095/150‑1 (MQ=11)
                          ggCAACTGGCCAGGCGGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGc                                                                                                                       <  2:116957/149‑1 (MQ=255)
                           gCAAGTGGCCAGAGAGGTGTAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCt                                                                                                                      <  1:506195/150‑1 (MQ=11)
                           gCAACTGGCGAGGCAGGTGGAGCCGGAATACAAACCTGCGCGAAAAGCTCATCTTGATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACTAGCt                                                                                                                      <  2:751054/150‑1 (MQ=17)
                                tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATg                                                                                                                 >  1:385316/1‑150 (MQ=255)
                                 ggCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                <  1:4417/150‑1 (MQ=255)
                                 ggCCAGACAGGTGGAGACGGATTCGATAGCGTCGCGGACAGATCTTTGAAAACACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                <  1:959694/150‑1 (MQ=14)
                                  gacAGACAGGTGGAGCCGGAAAACACACCGGCGCGAAAAGTTCATTTTTATCGCTGCGAGCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                               <  2:551920/148‑1 (MQ=14)
                                  gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgag                                                                                                               <  1:175277/150‑1 (MQ=255)
                                   ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                              <  1:760997/150‑1 (MQ=255)
                                    cAGACAGGTGGAGCCGGAATGCACGCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGGTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGAAACCGGCTTAATgagga                                                                                                             >  2:506453/1‑150 (MQ=17)
                                    cAGACAGGTGGAGCCGGAATACACACCGGCGCCAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                             <  1:129683/150‑1 (MQ=255)
                                       aCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                                <  2:470900/112‑1 (MQ=255)
                                       aCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                                >  1:470900/1‑112 (MQ=255)
                                          ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAACTTAATGAGGAGAAccc                                                                                                       >  1:13371/1‑150 (MQ=255)
                                           gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                      >  1:840164/1‑150 (MQ=255)
                                             ggAGCGGGAACACACGCCGGCGTGAAAGGGTTATCTTTATGACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATTAGGAGAAACCGca                                                                                                    <  2:127244/150‑1 (MQ=14)
                                              gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcat                                                                                                   <  1:945553/150‑1 (MQ=255)
                                                 ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatga                                                                                                                                        <  1:552425/110‑1 (MQ=11)
                                                 ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatga                                                                                                                                        >  2:552425/1‑110 (MQ=11)
                                                    gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTCTCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctg                                                                                             >  1:136998/1‑146 (MQ=11)
                                                          cacacCGGCGCGAAAAGCTCATCTTTATGACTGCGACGACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                        <  2:281365/150‑1 (MQ=255)
                                                          cacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                                                        >  1:942627/1‑150 (MQ=255)
                                                            cacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGc                                                                                      <  2:445194/150‑1 (MQ=255)
                                                                ggCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTa                                                                                  >  1:859650/1‑150 (MQ=255)
                                                                 gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATAAGCCGTAc                                                                                 >  1:427319/1‑150 (MQ=255)
                                                                     gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCCTAATGAGGAGAACCCGCATCATGTGTATAAGCCGTACCGTc                                                                             >  1:814822/1‑150 (MQ=255)
                                                                     gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGAGGCAACCAGGTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGGACCGTc                                                                             >  1:3072/1‑150 (MQ=255)
                                                                      aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAACACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCt                                                                            <  1:384907/150‑1 (MQ=255)
                                                                       aaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                                <  2:677126/81‑1 (MQ=255)
                                                                       aaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                                >  1:677126/1‑81 (MQ=255)
                                                                          gCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCa                                                                        <  2:859650/150‑1 (MQ=255)
                                                                           cttaTCTTTAACACTGCGACCACCAGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGACTATGATGAATGGGGCAACCAGCTTAATGGGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                       <  1:687963/147‑1 (MQ=17)
                                                                            tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATCCGGCGTGGAGCGCGGAATATGATGAATTGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAgg                                                                      >  1:770660/1‑150 (MQ=255)
                                                                            tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAgg                                                                      >  2:574055/1‑150 (MQ=255)
                                                                              aTCTTTTTCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                                <  1:408696/74‑3 (MQ=11)
                                                                              aTCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCCGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                                                >  2:408696/1‑72 (MQ=14)
                                                                                  ttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCACCATGAATATCAGCCGTACCGTCTGCCAGGgcagca                                                                <  1:733093/150‑1 (MQ=255)
                                                                                    aTCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                >  2:167472/1‑100 (MQ=255)
                                                                                    aTCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                                <  1:167472/100‑1 (MQ=255)
                                                                                            gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                      <  2:614432/150‑1 (MQ=255)
                                                                                               caccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTGCCGTCTGCCAGGGCAGCAGCATAATGAGGaa                                                   >  2:138666/1‑150 (MQ=255)
                                                                                                  cGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTa                                                                                  <  2:601583/116‑1 (MQ=255)
                                                                                                  cGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTa                                                                                  >  1:601583/1‑116 (MQ=255)
                                                                                                     ggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                             <  1:883094/150‑1 (MQ=255)
                                                                                                       aCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCt                                           <  2:59882/150‑1 (MQ=255)
                                                                                                          gCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                        <  2:556334/150‑1 (MQ=255)
                                                                                                           ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgc                                                                                                                                                   >  1:892162/1‑42 (MQ=255)
                                                                                                           ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgc                                                                                                                                                   >  1:409772/1‑42 (MQ=255)
                                                                                                           ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgc                                                                                                                                                   <  2:892162/42‑1 (MQ=255)
                                                                                                           ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgc                                                                                                                                                   <  2:409772/42‑1 (MQ=255)
                                                                                                                  ggcTTATCAGCGAAGGCGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                <  2:513246/149‑1 (MQ=255)
                                                                                                                   gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCg                               >  1:260145/1‑150 (MQ=255)
                                                                                                                          aGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCgg                        >  2:164712/1‑150 (MQ=255)
                                                                                                                           gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGt                       >  2:26262/1‑150 (MQ=255)
                                                                                                                                gACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAGTGAGGAGAACCCGCATCATGTGTATCAGCCGTAc                                                                                 >  1:370555/1‑87 (MQ=255)
                                                                                                                                gACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAGTGAGGAGAACCCGCATCATGTGTATCAGCCGTAc                                                                                 <  2:370555/87‑1 (MQ=255)
                                                                                                                                   ggCAGTACGGCGTGGAGCGCGGAATATGATGGATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGAt               >  2:828936/1‑150 (MQ=255)
                                                                                                                                   ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGAt               >  2:549263/1‑150 (MQ=255)
                                                                                                                                               gggAGCGCGGAATATGATGAATGGCGCACTGTGTTTAATGTCGACAATAAAAAACATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAgggg   <  1:960205/149‑1 (MQ=17)
                                                                                                                                                ggAGCGGGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGc  <  2:672605/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGC  >  NC_000913_3_pae_pgi/526647‑526939

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