Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi526,7960CG28.0% 103.1 / 50.6 76A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (26/28);  new base G (7/14);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCG  >  NC_000913_3_pae_pgi/526649‑526940
                                                                                                                                                   |                                                                                                                                                
ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                                                 <  1:781219/149‑1 (MQ=255)
 cgcgcAGTGCGTGCTGACGGTGAAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCgg                                                                                                                                                <  1:578701/148‑1 (MQ=255)
  gcgcAGTGCGGGCTGACGGCGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGa                                                                                                                                               <  1:481641/149‑1 (MQ=255)
    gcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGagt                                                                                                                                             >  1:44039/1‑147 (MQ=255)
       gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAACCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGGGGAATatg                                                                                                                                          >  1:562493/1‑149 (MQ=255)
        tGCGGGCTGACGGTGGGGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGGTTAGCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                         <  1:185879/149‑1 (MQ=255)
         gCGGGCTGACGGTGGAGCAACTCGCTAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGAGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                        <  1:790058/149‑1 (MQ=255)
          cGGGCTGACGGTGGAGCAACTGGGCAGACAGGCGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACTGTGGTCTGTTGCTGGAGTTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgatg                                                                                                                                       <  1:110724/149‑1 (MQ=14)
                 aCGGTGGAGCACATGGTCAGACAGGTGGAGCAGAAATCGTCACCGGCGCTAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGTGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgggg                                                                                                                                <  1:212087/149‑1 (MQ=17)
                    gTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCaa                                                                                                                             >  2:502994/1‑149 (MQ=255)
                       gAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAAAATGATGAATGGGGACACCa                                                                                                                          >  2:727039/1‑149 (MQ=255)
                       gAGCAACTGGCCAGACAGGTGGAGCCGGAATACAACCCGGCGCGAATAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAGTGGGGCAACCa                                                                                                                          <  2:563999/149‑1 (MQ=255)
                           aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGCATGATGAACGGCGCAACCAGAtt                                                                                                                      >  2:380011/1‑149 (MQ=11)
                           aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGaa                                                                                                                                              <  1:207478/125‑1 (MQ=255)
                           aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGaa                                                                                                                                              >  2:207478/1‑125 (MQ=255)
                            aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCGGCGAAGACGGCAATACAACGTGGGGCGCGGAATATGATGAATGGGGCAACCAGCTTa                                                                                                                     >  1:749393/1‑149 (MQ=255)
                              tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTTGCGCTTATCAGCGAATACGTTAATTTGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAAt                                                                                                                   <  1:200817/149‑1 (MQ=11)
                              tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGCATATGATGGATGGGGCAACCAGCTTatt                                                                                                                   >  2:427576/1‑147 (MQ=255)
                                gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGc                                                                                                                        <  2:120973/142‑1 (MQ=255)
                                gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGc                                                                                                                        >  1:120973/1‑142 (MQ=255)
                                 ccggccAGGTTGACCCGGAATACACAGCGGCGTGAAAAGCTCATCTTTATCACTGCCATCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                <  2:22028/144‑1 (MQ=17)
                                 ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                <  2:689762/149‑1 (MQ=255)
                                  gaGACAGGTGGAGCCGGAAGACACACCGGCGGGAAAAGCTCATCTTTAGCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                               <  2:327549/148‑1 (MQ=11)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGTCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                               <  1:332893/149‑1 (MQ=255)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGGCCGCGGCCTCACAGGAATAGCGCAGCCAGCGTAATgaag                                                                                                               >  2:688785/1‑147 (MQ=14)
                                  cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                               <  1:597276/149‑1 (MQ=255)
                                   aGACAGGTGGAGCCGGAATACCCACCGGCGCGAAAAGCTCATCTTGATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                              <  2:589500/149‑1 (MQ=11)
                                   aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATTCGGTGTGGAGTGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                              <  1:213597/149‑1 (MQ=255)
                                   aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                              <  1:214739/149‑1 (MQ=255)
                                   aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATATTTATCACTGCGACAACCGGGGACTGCCGGTGGCGCGTATCCGCGAAGACGGCCATACGGCGCGGACCGCGGACTATGAAGAATGGGAAACCCACGTTAAggagga                                                                                                              >  2:304949/1‑149 (MQ=14)
                                   aGAAAGGTGGAGCCGGAATACACCCCGGCGCGAAAAGCTCATATTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                              <  2:9302/149‑1 (MQ=255)
                                    gACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGATGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgaggag                                                                                                             <  1:345960/149‑1 (MQ=255)
                                           ggAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGTGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGc                                                                                                      <  1:643913/149‑1 (MQ=255)
                                              gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGCACGCGcatc                                                                                                   >  2:70489/1‑149 (MQ=255)
                                              gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatc                                                                                                   >  1:616962/1‑149 (MQ=255)
                                              gCCGGAATACACACCGGCGCGAAAAGCGCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTCATGAGGAGAACCCGcatc                                                                                                   >  2:520897/1‑149 (MQ=11)
                                                  gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGTCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacct                                                                                               <  1:169803/149‑4 (MQ=11)
                                                  gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGCGGCAACCCGCTTAATGAGGAGAACCCCcatcacct                                                                                               >  2:141824/1‑146 (MQ=11)
                                                      acactcCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctgcac                                                                                           <  2:418001/149‑8 (MQ=11)
                                                            tggCGCGAAAAGATCATCTTGATCTCTGCGACGACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCg                                                                                     <  2:34529/148‑1 (MQ=255)
                                                                gcgAAAAGCTCATCTTTATCACTGCGACCCCCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCAGCATCATGTGTATCAGCCGTaac                                                                                 >  1:304979/1‑147 (MQ=255)
                                                                  gAAAAGCTCATCTTTATCCCTGCGACCACCGGGGACTGCCGCTAAAACTTATCAACGAAGACGGCAATACGGCATGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                               >  1:210693/1‑149 (MQ=18)
                                                                        cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                     <  1:29507/121‑1 (MQ=255)
                                                                        cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                     >  2:29507/1‑121 (MQ=255)
                                                                        cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCa                                                                         <  2:786001/149‑1 (MQ=255)
                                                                         tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAAACGGCAATAAGGCATTGAGCGCGGAATAAGATTAATGGAGCCACCAACGTAATTGCGAGAACCACCATCATGTTTAACAATCGGACCGtattcaca                                                                        >  2:350618/1‑142 (MQ=255)
                                                                         tCATCTTTATCACTACGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCATGAATCGCGGAATATGATGAATGGGGAAACAAGCTTAATGAGGAGAAGCCGCATCAGGTGTCTCAGCCATCCCGTCTGCGAg                                                                        >  1:158647/1‑149 (MQ=17)
                                                                           aTCGTTACCACGGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAggg                                                                      <  2:585013/149‑1 (MQ=255)
                                                                              tttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAATCTGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                                                                   <  1:74665/149‑1 (MQ=255)
                                                                                tATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagca                                                                 >  2:666423/1‑149 (MQ=255)
                                                                                 gtCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAAACGGCAATACGGCGGAGAGCGCAGAATATGATGACTGGGGCAACCAGCTAAACGAGGAGAACCCGCATCATGTGTAACAGCCTTACCGTCTGCCAGGgcagcaa                                                                >  2:603720/2‑148 (MQ=17)
                                                                                    aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGGTGAATGGGGCAACCAGCTTAATGGGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCat                                                             >  2:501456/1‑149 (MQ=255)
                                                                                         gacgagcGGGGACGGGCGCTGGCGCTTATCAGCTAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAGCCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga                                                        <  2:478097/143‑1 (MQ=17)
                                                                                         gACCACCTGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAACCGCATCATGGGTATCAGCCATACCGTCTGCCAGGGCAGCAGCaagatgg                                                        >  2:245309/1‑148 (MQ=17)
                                                                                          accaccGGTGACTGCCGCTGGCGCTTATCAGCGAAAACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGCACCCCCATCATGTTTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAg                                                       <  2:733856/149‑1 (MQ=17)
                                                                                          accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGGCAGGGCAGCAGCATGATGAg                                                       <  2:591568/149‑1 (MQ=255)
                                                                                            taccGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATTAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGa                                                     <  2:8334/148‑1 (MQ=255)
                                                                                            ccccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAATCCGGATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGa                                                     <  1:245309/147‑1 (MQ=255)
                                                                                            caccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAgaa                                                     >  1:368534/1‑147 (MQ=255)
                                                                                            caccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGa                                                     >  1:465642/1‑149 (MQ=255)
                                                                                              ccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGACGTGAAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACAGTCTGGCAGGGCAGCAGCATGATGAGGaaa                                                   >  1:8334/1‑148 (MQ=255)
                                                                                                ggggACTGCCGGTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCa                                                 <  2:116504/149‑1 (MQ=255)
                                                                                                   gACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAggg                                              <  1:679935/149‑1 (MQ=255)
                                                                                                              gcgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATaa                                   <  1:502994/149‑1 (MQ=255)
                                                                                                                gcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                 >  1:321300/1‑149 (MQ=255)
                                                                                                                   tATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGGCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTc                              >  2:565257/1‑149 (MQ=255)
                                                                                                                       agggAAGACGGCAATACGGGGTGGATGGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCTATGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCg                          <  2:158647/146‑1 (MQ=17)
                                                                                                                        gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGATTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCgc                         >  1:177661/1‑148 (MQ=255)
                                                                                                                          gaagaaTGCAATGCGGAGTGGAGCGTGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGACACCGCATCATGTGTATCTGGCGTACCGTCTGCCAGGGCAGTAGCATGATGAGGAATCACGGCTGTACTATAACCGTCACCGGta                       <  1:723997/149‑1 (MQ=17)
                                                                                                                           aaGACGGCAATACGGCCTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtac                      <  1:666423/149‑1 (MQ=255)
                                                                                                                                ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCAc                            >  2:477880/1‑138 (MQ=255)
                                                                                                                                ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCAc                            <  1:477880/138‑1 (MQ=255)
                                                                                                                                ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACgc                 >  2:391841/1‑148 (MQ=255)
                                                                                                                                 tcAATACGCCGGGGAGCGGGTAATATGTGGAATGGGTCAACCAGCTGAATGAGGAGAGCCCACATCATGTGTATCAGCCGCCCCGTCTGCAAGGTGCTCTGGTGGATGAGGAATCAGGGCTGTATTATAACCGTCACCGGTACTACGAt                <  2:323961/148‑1 (MQ=255)
                                                                                                                                     tACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCATCAGCATGATTAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGt            <  1:90139/149‑1 (MQ=255)
                                                                                                                                       cGGCGTGGAGCGCGGAATGTGATGAATGGGGCAACCAGCTAAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTg          >  1:623852/1‑149 (MQ=255)
                                                                                                                                               aGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGATGAGAACCCGCATCACGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTAAAACCGTCACCGGAACTACGATCCGTTGCAGGGGCg  <  1:133766/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCG  >  NC_000913_3_pae_pgi/526649‑526940

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