Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi526,8000AG24.6% 111.2 / 30.8 71E1180E (GAA→GAG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (24/28);  new base G (5/12);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCG  >  NC_000913_3_pae_pgi/526656‑526940
                                                                                                                                                |                                                                                                                                            
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                                                                                                                                              >  2:207478/1‑125 (MQ=255)
                    aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGaa                                                                                                                                              <  1:207478/125‑1 (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)
                                                                                                                                                |                                                                                                                                            
GTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCG  >  NC_000913_3_pae_pgi/526656‑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: ‑