Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi464,3560AT25.0% 136.6 / 51.8 68intergenic (+48/‑46)ybaV/fadMputative competence‑suppressing periplasmic helix‑hairpin‑helix DNA‑binding protein/long‑chain acyl‑CoA thioesterase III
Reads supporting (aligned to +/- strand):  ref base A (22/29);  new base T (6/11);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGG  >  NC_000913_3_pae_pgi/464207‑464497
                                                                                                                                                     |                                                                                                                                             
cgcgAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCa                                                                                                                                               <  3:95877/150‑1 (MQ=255)
 gcgAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAg                                                                                                                                              <  4:524443/150‑1 (MQ=255)
    aaGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                                                           <  1:520117/150‑1 (MQ=255)
     agagTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGACGAGGTCATACCAGTTAt                                                                                                                                          <  2:917475/150‑1 (MQ=255)
       agTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGa                                                                                                                                        <  1:3536/150‑1 (MQ=255)
          aCGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     <  3:580329/150‑1 (MQ=255)
          aCGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     <  1:983816/150‑1 (MQ=255)
           cGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctc                                                                                                                                    <  2:216602/150‑1 (MQ=255)
            ggTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACAATGTAATCAACTTGCATAGCCGCAATTTTTTGACAGACAGAAGCAGTCATACCAGTTATGACatct                                                                                                                                   >  1:722202/1‑150 (MQ=255)
               ccGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttca                                                                                                                                <  4:120206/150‑4 (MQ=255)
               ccGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACctcttca                                                                                                                                <  3:465216/150‑4 (MQ=255)
                 ggttAAAGCCGTGGAGGATCTAAAGCCGGTGCCGGGGGTGGGCAATTTGCTGGTGGAACGTAATCTGGCGGTATTAACACTCGAATTAATTTGCATAGTGGCAATCTTTTGCCAGACTGCAGAGGTCATCACTGGTATGATctctccagt                                                                                                                              <  2:542585/148‑6 (MQ=255)
                  tttAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACtt                                                                                                                             <  2:554388/150‑1 (MQ=255)
                            tgGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaacaa                                                                                                                   <  2:33568/150‑1 (MQ=255)
                              gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaa                                                                                                                         >  1:179414/1‑142 (MQ=255)
                              gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaa                                                                                                                         <  2:179414/142‑1 (MQ=255)
                              gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACg                                                                                                                 >  1:400717/1‑150 (MQ=255)
                              gaggaTCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACg                                                                                                                 >  1:549866/1‑150 (MQ=255)
                                     tAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTGATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGtt                                                                                                          <  2:885543/150‑1 (MQ=255)
                                      aaaGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTa                                                                                                         >  3:15708/1‑150 (MQ=255)
                                      aaaGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTa                                                                                                         <  4:15708/150‑1 (MQ=255)
                                            ggTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     >  3:42577/1‑116 (MQ=255)
                                            ggTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     >  3:42612/1‑116 (MQ=255)
                                            ggTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     <  4:42577/116‑1 (MQ=255)
                                            ggTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATAACTGGTATGACct                                                                                                                                     <  4:42612/116‑1 (MQ=255)
                                                  ggggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCa                                                                                             <  3:144132/150‑1 (MQ=255)
                                                   gggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaaca                                                                                                                    <  3:434384/126‑1 (MQ=255)
                                                   gggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATaacaaca                                                                                                                    >  4:434384/1‑126 (MQ=255)
                                                   gggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATGACAACAACGTAAGGTTATTGCGCTATGCaa                                                                                            <  1:518588/150‑1 (MQ=255)
                                                    ggATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                           <  2:531923/150‑1 (MQ=255)
                                                          gCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                                                           <  1:704700/96‑1 (MQ=255)
                                                          gCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTa                                                                                                                                           >  2:704700/1‑96 (MQ=255)
                                                              ttCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCaa                                                                                 >  1:21076/1‑150 (MQ=255)
                                                                     gtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACct                                                                                                                                     >  3:146828/1‑91 (MQ=255)
                                                                     gtggAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACct                                                                                                                                     <  4:146828/91‑1 (MQ=255)
                                                                         aaCGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGat                                                                      <  4:514950/150‑1 (MQ=255)
                                                                             tgaTCTGGCGGTCTTAACCCTGTAATTAATTTGCATAGTGGCAATTCTTTGCCAGACTGAAGAGGTCATACTAGTTATGGGCTCGGTAATTATATCAACAACGTAAGGTTATTGCGCTATGCAAACACAAAACACAGTTCGTGGTTatca                                                                  <  2:243916/148‑1 (MQ=255)
                                                                                tCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCAtct                                                               <  2:408660/150‑1 (MQ=255)
                                                                                   ggCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGa                                                            <  2:400717/150‑1 (MQ=255)
                                                                                    gCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGAc                                                           <  3:428019/150‑1 (MQ=255)
                                                                                         aTTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                                                                          >  1:329975/1‑114 (MQ=255)
                                                                                         aTTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAac                                                                                          <  2:329975/114‑1 (MQ=255)
                                                                                           tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTAc                                                                                                                               <  1:374331/75‑1 (MQ=255)
                                                                                           tAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTAc                                                                                                                               >  2:374331/1‑75 (MQ=255)
                                                                                             aCCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAg                                                  >  2:654222/1‑150 (MQ=255)
                                                                                               ccTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCa                                                >  1:849764/1‑150 (MQ=255)
                                                                                                                     gCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCtt                          >  1:1021501/1‑150 (MQ=255)
                                                                                                                       aaTTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGa                        >  3:539428/1‑150 (MQ=255)
                                                                                                                       aaTTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGa                        >  4:304007/1‑150 (MQ=255)
                                                                                                                             ttGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTtct                  <  1:654222/150‑1 (MQ=255)
                                                                                                                                 cAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaa                                          <  3:238319/122‑1 (MQ=255)
                                                                                                                                 cAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaa                                          >  4:238319/1‑122 (MQ=255)
                                                                                                                                  aGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaag             >  3:433819/1‑150 (MQ=255)
                                                                                                                                  aGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaag             <  2:549866/150‑1 (MQ=255)
                                                                                                                                     cTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaa                                          >  1:752958/1‑118 (MQ=255)
                                                                                                                                     cTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTcaa                                          <  2:752958/118‑1 (MQ=255)
                                                                                                                                       acagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                           <  1:827059/65‑1 (MQ=25)
                                                                                                                                       acagagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaaa                                                                                           >  2:827059/3‑67 (MQ=25)
                                                                                                                                         agagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCg      >  4:463709/1‑150 (MQ=255)
                                                                                                                                          gagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaagaa           >  3:104828/1‑144 (MQ=255)
                                                                                                                                          gagGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCgaagaa           <  4:104828/144‑1 (MQ=255)
                                                                                                                                           agGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCt    <  3:438774/150‑1 (MQ=255)
                                                                                                                                           agGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAAcaca                                                                                       >  2:850160/1‑67 (MQ=37)
                                                                                                                                           agGTCATAACTGGTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAAcaca                                                                                       <  1:850160/67‑1 (MQ=37)
                                                                                                                                             gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaa                                                                                            <  2:474472/60‑1 (MQ=255)
                                                                                                                                             gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCaa                                                                                            >  1:474472/1‑60 (MQ=255)
                                                                                                                                             gTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTgg  >  2:925072/1‑150 (MQ=255)
                                                                                                                                                aTACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTc                                                                            <  2:662293/73‑1 (MQ=255)
                                                                                                                                                aTACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTc                                                                            >  1:662293/1‑73 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
CGCGAAGAGTACGGTCCGTTTAAAACTGTGGAGGATCTAAAGCAGGTGCCGGGGATGGGCAATTCGCTGGTGGAACGTAATCTGGCGGTATTAACCCTGTAATTAATTTGCATAGTGGCAATTTTTTGCCAGACTGAAGAGGTCATACCAGTTATGACCTCTGTACTTATAACAACAACGTAAGGTTATTGCGCTATGCAAACACAAATCAAAGTTCGTGGATATCATCTCGACGTTTACCAGCACGTCAACAACGCCCGCTACCTTGAATTTCTCGAAGAAGCCCGCTGG  >  NC_000913_3_pae_pgi/464207‑464497

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