Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,781,1520AG25.4% 111.8 / 13.3 60intergenic (+135/‑63)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (35/9);  new base G (0/15);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 3.28e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCT  >  NC_000913_3_pae_pgi/3781005‑3781300
                                                                                                                                                     |                                                                                                                                                     
gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAgg                                                                                                                                                     <  1:529480/150‑1 (MQ=255)
gAATGCGGCATAGTCGGTTGCGCCCCCGCACTCTATCCGCCTTCCTGTGGGAGGGGTGAGCGGTGGGCGGTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAgg                                                                                                                                                     <  3:648595/150‑1 (MQ=255)
 aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAggg                                                                                                                                                    <  3:455992/150‑1 (MQ=255)
  aTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAgggg                                                                                                                                                   >  2:918853/1‑150 (MQ=255)
     cGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTAAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑aaa                                                                                                                                               >  2:825650/1‑150 (MQ=255)
      ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑aaaa                                                                                                                                              >  2:682994/1‑150 (MQ=255)
      ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑aaaa                                                                                                                                              >  1:967180/1‑150 (MQ=255)
        cATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAgg                                                                                                                                            <  4:184772/150‑1 (MQ=255)
                tttGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTg                                                                                                                                    >  2:714277/1‑150 (MQ=255)
                     gcCCCCTGACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCAc                                                                                                                               >  2:1155489/1‑150 (MQ=255)
                          cTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAAtt                                                                                                                          >  2:998185/1‑150 (MQ=255)
                                            cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCGTCCCTCTGAAACGACACCGCGCTCTTTTTTTCTCCCTCGCCCCTCC‑GGGCAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                                                                                        >  4:158061/1‑150 (MQ=255)
                                                       gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                                                                             >  4:253937/1‑150 (MQ=255)
                                                            cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                                                                                        >  3:588831/1‑150 (MQ=255)
                                                              ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AACAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACCACGGGGGCAGTAt                                                                                      >  2:1063342/1‑150 (MQ=255)
                                                                  gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTTTTTTTTTTTTCTCC‑GGAGAAAGGGCA‑AAGATGAGGGA‑AAAGGGCCGCAATGCTCACCAATTTTCTTTTCTCCTGCCCCCGCACTAAGGGGGCATTCTGCTa                                                                                  >  4:250398/1‑150 (MQ=255)
                                                                  gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                                                                                  >  3:423368/1‑150 (MQ=255)
                                                                     cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                                                                               >  1:763771/1‑150 (MQ=255)
                                                                             aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTTTTTTTTTCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt                                                                       >  2:721985/1‑150 (MQ=255)
                                                                             aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt                                                                       >  1:33994/1‑150 (MQ=255)
                                                                               tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt                                                                     >  1:581921/1‑150 (MQ=255)
                                                                               tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt                                                                     >  1:145724/1‑150 (MQ=255)
                                                                                 cctccctcGGAACCGCCCCCGCCCTTTTTTTTCCCCCCCCCCCCCCCC‑GGGGGGGGGGCG‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc                                                                   <  4:82479/150‑1 (MQ=255)
                                                                                 ccccACCCTGAACCGACCCCGCCCTCTTTTTTTCCCCCCGGCCCCTCG‑GGGGGGGGGGCG‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc                                                                   <  4:83087/147‑1 (MQ=255)
                                                                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGTCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc                                                                   >  2:524649/1‑150 (MQ=255)
                                                                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc                                                                   >  4:132053/1‑150 (MQ=255)
                                                                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GCGGAGAGGGCC‑GGGGGGAGGTG‑ACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCCCTAAGCGGGCAGTAGGCTAAACATCGTACTTGAc                                                                   >  4:524412/1‑150 (MQ=255)
                                                                                   tCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa                                                                 >  3:462934/1‑150 (MQ=255)
                                                                                   tCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa                                                                 >  3:693076/1‑150 (MQ=255)
                                                                                   cccctctGAACCGCCCCCCCCCTTTTTTTTCCCCCCCCCCCCCCCG‑GGGGGGGGGGCG‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa                                                                 <  1:1002393/147‑1 (MQ=255)
                                                                                     ccccggAAACGACCCCGCACTCTTTTTTTCCCCCCCCCCCGCCG‑GGGGGGGGGGCG‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAc                                                               <  4:154266/145‑1 (MQ=255)
                                                                                         tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGAGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTAGCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGa                                                           >  1:315515/1‑150 (MQ=255)
                                                                                          gAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCACTTTTCTGCTATCATGCCCCCGCCCTACGGGGGCAGCATGATAACCCTCGTACTTTACGAACCAGaa                                                          >  1:928411/1‑150 (MQ=255)
                                                                                               gACCCCGCTCTTTTTTTTCCCCCCTCCCCCCCCG‑GGGGAGGGGGCG‑GGGGGGGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTcc                                                     <  1:825650/150‑1 (MQ=255)
                                                                                               gACCCCGCACTCTTTTTTTCCCCCTCGCCCCCCG‑GGGGGGAGGGCG‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTcc                                                     <  4:125242/150‑1 (MQ=255)
                                                                                               gACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTcc                                                     >  2:389707/1‑150 (MQ=255)
                                                                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCa                                                 >  3:370915/1‑150 (MQ=255)
                                                                                                    ccccctttttttttCCCCCCCCCCCCCCCGGGGGAGAGGGCG‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCa                                                 <  3:587531/143‑1 (MQ=255)
                                                                                                     gccctctTTTTTGCTCCCCCGCCCCCCG‑GGGGGGGGGGCG‑GGGGTGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCaaa                                               <  4:693076/147‑1 (MQ=255)
                                                                                                          tcttttttCCCCCCCCCCCCCCC‑GGGGGGGGGGCCGGGGGTGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACt                                           <  1:1155489/148‑1 (MQ=255)
                                                                                                           tttttttCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGTCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTgg                                         <  1:524649/150‑1 (MQ=255)
                                                                                                            ttttttCCCCCTCCCCCCCCCGGGGGAGGGGCCC‑GGGGGGGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTgg                                         <  1:918853/150‑1 (MQ=255)
                                                                                                                 tctcCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAAcat                                   >  3:517209/1‑150 (MQ=255)
                                                                                                                      cgcGCCCCTGC‑GGGGGGAGGGCC‑GGGGGGGGGGT‑AAAGGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATcc                              <  4:525452/148‑1 (MQ=255)
                                                                                                                       tCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCg                             >  3:607957/1‑150 (MQ=255)
                                                                                                                           cccTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCt                                                                                                                                  <  1:453621/45‑1 (MQ=255)
                                                                                                                           cccTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCt                                                                                                                                  >  2:453621/1‑45 (MQ=255)
                                                                                                                             cagc‑GGGGAGAGGGGC‑GGGGGGAGGGT‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACTAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTg                       <  3:250398/147‑1 (MQ=255)
                                                                                                                               gcgggggAGGGGGCC‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTgc                      <  4:462934/147‑1 (MQ=255)
                                                                                                                                ggggggAGGGGGCG‑GGGGGGGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTgcg                     <  3:132601/148‑1 (MQ=255)
                                                                                                                                c‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTgcgc                    >  3:453763/1‑150 (MQ=255)
                                                                                                                                 gggggAGGGGGCG‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTgcgc                    <  4:404939/149‑1 (MQ=255)
                                                                                                                                   gggAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTa                  >  2:69738/1‑150 (MQ=255)
                                                                                                                                    ggggggggCG‑GGGGGGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTaa                 <  4:370915/145‑1 (MQ=255)
                                                                                                                                        agGGCT‑GGGGTGAGGGG‑AAGAGGCCGCCCTGCCCACCAATTTTCTGCTATCCTGCGCCCGCACTAACGTCGCGGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTCGCAACATCATCCGTCTTTGCGCTAATAcc             >  4:645130/1‑150 (MQ=255)
                                                                                                                                          ggCC‑GGGGTTCGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCg            <  3:712403/150‑1 (MQ=255)
                                                                                                                                          ggCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCgg           >  1:573677/1‑150 (MQ=255)
                                                                                                                                                gggTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTc      <  2:763771/150‑1 (MQ=255)
                                                                                                                                                   tGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTc   >  1:342042/1‑150 (MQ=255)
                                                                                                                                                    gAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCt  <  4:588831/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCT  >  NC_000913_3_pae_pgi/3781005‑3781300

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