Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi526,8000AG26.8% 154.7 / 43.8 82E1180E (GAA→GAGrhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (25/35);  new base G (12/10);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATAT  >  NC_000913_3_pae_pgi/526655‑526946
                                                                                                                                                 |                                                                                                                                                  
aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAgtat                                                                                                                                                 >  2:212623/1‑149 (MQ=255)
  tGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatga                                                                                                                                               <  1:970353/149‑1 (MQ=255)
  tGCGGGCTGACGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga                                                                                                                                               <  1:666370/149‑1 (MQ=255)
   gCGGGCTGACGCTGGACAAACTAACCAGACGACTGAGGACGGAACACACACCGGCGAGAAAAGCCCATCTTTATCACCACGACCACCGGGCACTGCCACTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatgat                                                                                                                                              <  2:256805/149‑1 (MQ=2)
      ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTACCAGCGAAGACGACACTACGGCATGGAACGCGCAAGatgatgga                                                                                                                                           >  2:192861/1‑147 (MQ=14)
         tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                                                        <  2:425349/149‑1 (MQ=255)
           gcGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgggg                                                                                                                                      <  1:764175/148‑1 (MQ=255)
           aCGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgggg                                                                                                                                      <  1:764243/149‑1 (MQ=255)
              gTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCaa                                                                                                                                   >  1:859756/1‑149 (MQ=255)
               tGGGGCAACTGGCCACACAGGGGGAGCGTGAATACACACCGTCGCGAAAATCTCATCTTTATCACTGCGACCACCGGGGACTTCCGCTTGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                                  <  2:592015/149‑1 (MQ=17)
                 gAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGGGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCa                                                                                                                                <  2:41030/149‑1 (MQ=255)
                  aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAg                                                                                                                               <  1:856017/149‑1 (MQ=255)
                     aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCtt                                                                                                                            <  2:979799/149‑1 (MQ=255)
                         ggCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATg                                                                                                                        >  2:825741/1‑149 (MQ=255)
                          gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATga                                                                                                                       <  2:663109/149‑1 (MQ=255)
                           ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                      <  1:212623/149‑1 (MQ=255)
                           ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                                      <  1:783398/149‑1 (MQ=255)
                            cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGCATGATGAATGGGGCCACCAGCTTACTgagg                                                                                                                     >  2:544662/1‑149 (MQ=11)
                             aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                    <  1:121445/149‑1 (MQ=255)
                             aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGGATGATGAATGGGGCGACCAGCTTAATgagga                                                                                                                    >  1:779585/1‑149 (MQ=11)
                             aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                                                    <  1:844904/149‑1 (MQ=255)
                              gACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGACTGGGGCAACCAGCTTAaggagggg                                                                                                                   >  2:963286/1‑147 (MQ=11)
                              gACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgaggag                                                                                                                   >  1:658949/1‑149 (MQ=255)
                                 aGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATTGAGAATGGTGCAACCAGCATAATGAGGAGAAc                                                                                                                >  2:922965/1‑149 (MQ=11)
                                   gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGGATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                              >  2:57730/1‑149 (MQ=255)
                                   gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                                              >  1:250037/1‑149 (MQ=255)
                                   gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCAGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACAAGCTTAATGAGGAGAAcac                                                                                                              >  1:938114/1‑147 (MQ=11)
                                     ggAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATAACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGc                                                                                                            <  2:795394/149‑1 (MQ=255)
                                      gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                           <  1:739077/149‑1 (MQ=255)
                                          cGGAATACACACCGGCGCGAAAAGCGCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcac                                                                                                       >  1:977653/1‑148 (MQ=11)
                                           ggAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacc                                                                                                      >  2:868625/1‑147 (MQ=255)
                                            gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacct                                                                                                     <  1:584356/149‑4 (MQ=255)
                                            gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtgt                                                                                                     <  1:944148/149‑1 (MQ=255)
                                                acacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctgcac                                                                                                 <  2:1009806/149‑8 (MQ=11)
                                                acacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTAt                                                                                                 <  1:1041901/149‑1 (MQ=255)
                                                    cccGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGc                                                                                             <  1:556178/148‑1 (MQ=255)
                                                       ggCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                                      <  2:587948/105‑1 (MQ=14)
                                                       ggCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgggg                                                                                                                                      >  1:587948/1‑105 (MQ=14)
                                                          gcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTAcc                                                                                       <  1:595776/149‑1 (MQ=255)
                                                                  cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGc                                                                                                            <  2:968642/120‑1 (MQ=11)
                                                                  cTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGc                                                                                                            >  1:968642/1‑120 (MQ=11)
                                                                   tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGTGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                                                              >  2:738481/1‑149 (MQ=255)
                                                                   tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTGCCGTCTGCCAg                                                                              <  1:664570/149‑1 (MQ=255)
                                                                    cATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAgg                                                                             >  1:208309/1‑149 (MQ=255)
                                                                    cATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAgg                                                                             >  1:332759/1‑149 (MQ=255)
                                                                     aTCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                                     <  2:984285/140‑1 (MQ=255)
                                                                     aTCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                                     >  1:984285/1‑140 (MQ=255)
                                                                     aTCGTTATCAGTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAggg                                                                            <  2:1011112/149‑1 (MQ=255)
                                                                       cTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgca                                                                          <  2:114699/149‑1 (MQ=255)
                                                                        tttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCTGAAGATGATGAATCAGGCAACCAGCCTAAAGAGGAGAACCCGCATCAAGGGACACAGCCCTACCGTCTGCCATGgcac                                                                         >  2:412639/1‑148 (MQ=14)
                                                                          tATCACCGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagca                                                                       <  1:1009566/149‑1 (MQ=255)
                                                                                     ccaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATgtgt                                                                                                   <  1:539458/110‑1 (MQ=255)
                                                                                     ccaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATgtgt                                                                                                   >  2:539458/1‑110 (MQ=255)
                                                                                     ccaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCATGGCGCACATAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAgg                                                            >  1:580720/1‑149 (MQ=17)
                                                                                      caccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                                  <  1:516602/78‑1 (MQ=255)
                                                                                      caccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                                                  >  2:516602/1‑78 (MQ=255)
                                                                                                  ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAc                                                                                                                >  2:249490/1‑84 (MQ=255)
                                                                                                  ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAc                                                                                                                <  1:249490/84‑1 (MQ=255)
                                                                                                  ccGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                               >  1:606396/1‑149 (MQ=255)
                                                                                                   cGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTAc                                              >  1:827374/1‑149 (MQ=255)
                                                                                                        gcgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATaa                                         >  2:116164/1‑149 (MQ=255)
                                                                                                            ttATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAat                                                                                                                                                   <  1:61401/39‑1 (MQ=255)
                                                                                                            ttATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAat                                                                                                                                                   >  2:61401/1‑39 (MQ=255)
                                                                                                            ttATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTTCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAGCCGt                                     <  1:854746/149‑1 (MQ=255)
                                                                                                               tCAGCGAAGACGGCAATACGGGGTGGAGCGCGGAATATTATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCGTCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCAc                                  <  2:585518/149‑1 (MQ=255)
                                                                                                               tCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCAc                                  >  1:231681/1‑149 (MQ=255)
                                                                                                                cAGCGAAGACGGCAATACGGTGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCAcc                                 <  1:738481/149‑1 (MQ=255)
                                                                                                                 aGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATCAAGGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCg                                <  2:888062/149‑1 (MQ=255)
                                                                                                                  gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                       <  1:653962/141‑1 (MQ=255)
                                                                                                                  gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAcc                                       >  2:653962/1‑141 (MQ=255)
                                                                                                                    gAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCCTCTGCCAGGGCAGCAGAATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGtc                             >  2:393001/1‑148 (MQ=255)
                                                                                                                     aaGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGc                                                                                                            <  2:469689/69‑1 (MQ=255)
                                                                                                                     aaGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGc                                                                                                            >  1:469689/1‑69 (MQ=255)
                                                                                                                      aGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGTCTGTACTATAACCGTCCACGGtact                           >  2:706616/1‑149 (MQ=255)
                                                                                                                       gACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCCGCCGTACCGTCTGCCAGGGCCGCAGCATGATGAGGAATGCGGGCTGTCCTATATCTGTCACCGGtacca                          >  2:36800/1‑147 (MQ=17)
                                                                                                                        aCGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                           <  2:123125/67‑1 (MQ=255)
                                                                                                                        aCGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGca                                                                                                           >  1:123125/1‑67 (MQ=255)
                                                                                                                          ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAcc                                                                                                               >  1:420371/1‑61 (MQ=255)
                                                                                                                          ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAAcc                                                                                                               <  2:420371/61‑1 (MQ=255)
                                                                                                                          ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGa                       >  2:426721/1‑149 (MQ=255)
                                                                                                                           gCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGAt                      >  1:71692/1‑149 (MQ=255)
                                                                                                                           gCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGAt                      <  2:208309/149‑1 (MQ=255)
                                                                                                                                             cgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGtat    <  2:606396/149‑1 (MQ=255)
                                                                                                                                             cgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGtat    <  1:561830/149‑1 (MQ=255)
                                                                                                                                               gAATATGATGAATGGGGCAACCAGCTTAATGAGTAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAAGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTatat  <  1:766071/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATAT  >  NC_000913_3_pae_pgi/526655‑526946

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: