Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8570GA42.7% 53.2 / 85.7 82A278A (GCC→GCTcynRtranscriptional activator of cyn operon, autorepressor
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base A (16/19);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGG  >  NC_000913/357727‑357995
                                                                                                                                  |                                                                                                                                          
tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCggcggc                                                                                                                                       >  2:964765/1‑136 (MQ=255)
tGAACTAGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCggcggc                                                                                                                                       <  1:964765/136‑1 (MQ=255)
 gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTc                                                                                                                                   >  2:218668/1‑139 (MQ=255)
  aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCt                                                                                                                                  <  1:537887/139‑1 (MQ=255)
       gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAg                                                                                                                             >  2:1028672/1‑139 (MQ=255)
          ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTg                                                                                                                          >  1:379041‑M2/1‑138 (MQ=255)
                     cGTAGCGCTTCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgcc                                                                                                               <  1:1110885‑M2/139‑13 (MQ=255)
                     cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgcc                                                                                                               <  2:649485‑M2/139‑13 (MQ=255)
                     cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                               >  1:634011/1‑139 (MQ=255)
                       tAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagt                                                                                                                      <  2:721170‑M2/130‑6 (MQ=255)
                       tAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagt                                                                                                                      >  1:721170‑M2/1‑125 (MQ=255)
                         gcgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggat                                                                                                           >  2:842973‑M2/1‑123 (MQ=255)
                          cgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggat                                                                                                           >  2:237044‑M2/1‑122 (MQ=255)
                          cgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggat                                                                                                           <  1:237044‑M2/138‑17 (MQ=255)
                          cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                >  2:1103806/1‑133 (MQ=255)
                          cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                <  1:1103806/133‑1 (MQ=255)
                           gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                         <  2:634011/139‑1 (MQ=255)
                               cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                         <  2:616846/135‑1 (MQ=255)
                               cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                         >  1:616846/1‑135 (MQ=255)
                               cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                     >  2:85876/1‑139 (MQ=255)
                                     aTTCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGtt                                                                                               <  1:1028672/139‑1 (MQ=255)
                                       tCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtttt                                                                                                                      >  2:219205/1‑114 (MQ=255)
                                       tCCAGCCGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtttt                                                                                                                      <  1:219205/114‑1 (MQ=255)
                                           gCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttat                                                                                         >  1:649485‑M2/1‑105 (MQ=255)
                                            ccGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTc                                                                                                                                   <  1:546205/96‑1 (MQ=255)
                                            ccGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTc                                                                                                                                   >  2:546205/1‑96 (MQ=255)
                                            ccGCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCa                                                                                        >  1:1003374/1‑139 (MQ=255)
                                              gCAGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAGCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                      <  2:914782/139‑1 (MQ=255)
                                                tgatctgtgTCAGCGGTTACCGTGATTGATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccgcatgcggcgtgaacgccttatccggc                                                                                    <  1:974924‑M2/135‑40 (MQ=255)
                                                aGACCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                    >  2:1197296/1‑139 (MQ=255)
                                                       tgtCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                       <  1:171195/113‑1 (MQ=255)
                                                       tgtCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                       >  2:171195/1‑113 (MQ=255)
                                                       tgtCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                            >  1:784065/1‑124 (MQ=255)
                                                       tgtCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                            <  2:784065/124‑1 (MQ=255)
                                                            gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGcc                                                                                                                               <  1:928320/84‑1 (MQ=255)
                                                            gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGcc                                                                                                                               >  2:928320/1‑84 (MQ=255)
                                                            gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggcctacgtagagca                                                                        >  2:586347‑M2/1‑88 (MQ=255)
                                                            gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggcctacgtagagca                                                                        <  2:769941‑M2/139‑52 (MQ=255)
                                                            gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggcctacgtagagca                                                                        <  1:542207‑M2/139‑52 (MQ=255)
                                                             cggttACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                          >  1:62703/5‑120 (MQ=255)
                                                             cGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                          <  2:62703/120‑1 (MQ=255)
                                                             cGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggcctacgtagagcac                                                                       <  1:11479‑M2/139‑53 (MQ=255)
                                                              ggttaccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggcctacgtagagcact                                                                      >  2:489251‑M2/4‑86 (MQ=255)
                                                              ggTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggcctacgtagagcact                                                                      <  1:489251‑M2/139‑54 (MQ=255)
                                                                 tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggcctacgtagagcactgaa                                                                   <  2:210264‑M2/139‑57 (MQ=255)
                                                                 tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtg                                                                                                   >  2:1424209‑M2/1‑83 (MQ=255)
                                                                 tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtg                                                                                                   <  1:1424209‑M2/107‑25 (MQ=255)
                                                                       gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggcctacgtagagcactgaactcgta                                                             >  2:11479‑M2/1‑77 (MQ=255)
                                                                        aTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggc                                                                                    >  2:145179‑M2/1‑76 (MQ=255)
                                                                        aTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggc                                                                                    <  1:145179‑M2/115‑40 (MQ=255)
                                                                          tCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAGGTGCTTATATGCCTTCGCAGCGGCTGTCTACCAGTTgtattttatgccggatgcggcgtgaacgacttatccggcctacgtagagcactgaactcgtaggc                                                          <  2:573045‑M2/139‑66 (MQ=255)
                                                                             tttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                       >  1:1243480/1‑139 (MQ=255)
                                                                              ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGtt                                                                                               >  2:78875/1‑98 (MQ=255)
                                                                              ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGtt                                                                                               <  1:78875/98‑1 (MQ=255)
                                                                                 cGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAACCTTAAGCCCGtt                                                   >  2:602794/1‑138 (MQ=255)
                                                                                   ccAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttat                                                                                         <  1:1051703‑M2/99‑35 (MQ=255)
                                                                                   ccAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttat                                                                                         >  2:1051703‑M2/1‑65 (MQ=255)
                                                                                      acaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAACCCGTCAtgtt                                              >  1:353543/1‑139 (MQ=255)
                                                                                        aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGAAAAAGATTTAAGCCCGTCAtgttgt                                            >  2:181783/1‑139 (MQ=255)
                                                                                          ccGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggcctacgtagagcactgaactcgtaggcctgataagcgtagcgc                                          <  2:921211‑M2/139‑82 (MQ=255)
                                                                                            gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggcctacgtagagcactgaactcgta                                                             >  1:1452854‑M2/1‑56 (MQ=255)
                                                                                            gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCCGCCAGCTgtagtttatgccggatgcggcgtgaacgccttatccggcctacgtagagcactgaactcgta                                                             <  2:1452854‑M2/118‑63 (MQ=255)
                                                                                             cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCgcgg                                                                                                                                           <  2:435191/39‑1 (MQ=255)
                                                                                             cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCgcgg                                                                                                                                           >  1:435191/1‑39 (MQ=255)
                                                                                             cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCt                                                                                                                                      <  2:1145341/44‑1 (MQ=255)
                                                                                             cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCt                                                                                                                                      >  1:1145341/1‑44 (MQ=255)
                                                                                             cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                          >  1:113458/1‑120 (MQ=255)
                                                                                             cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                          <  2:113458/120‑1 (MQ=255)
                                                                                             cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                       <  2:353543/139‑1 (MQ=255)
                                                                                             cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGcc                                                                                                                               <  2:624553/51‑1 (MQ=255)
                                                                                             cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGcc                                                                                                                               >  1:624553/1‑51 (MQ=255)
                                                                                             cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtga                                                                                                  >  2:422125‑M2/1‑55 (MQ=255)
                                                                                             cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTgtagtttatgccggatgcggcgtga                                                                                                  <  1:422125‑M2/80‑26 (MQ=255)
                                                                                                          cGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                          <  2:1243480/139‑1 (MQ=255)
                                                                                                            ccATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                         <  2:82193/58‑1 (MQ=255)
                                                                                                            ccATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                         >  1:82193/1‑58 (MQ=255)
                                                                                                               tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                 <  2:919435/47‑1 (MQ=255)
                                                                                                               tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAATAGCGTTAAGCCCGTTATGTTGTGTGGCAATCGATGCTGGTAaaaa                     >  2:243126/1‑139 (MQ=38)
                                                                                                               tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                    >  2:436650/1‑92 (MQ=255)
                                                                                                               tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                    <  1:436650/92‑1 (MQ=255)
                                                                                                               tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                     >  1:386876/1‑139 (MQ=255)
                                                                                                               tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTCTTccgc                                                                                                                 >  1:919435/1‑47 (MQ=255)
                                                                                                                  gCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                        <  2:610971/101‑1 (MQ=255)
                                                                                                                  gCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                        >  1:610971/1‑101 (MQ=255)
                                                                                                                           ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCag                                                                                       >  1:1321511/1‑61 (MQ=255)
                                                                                                                           ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCag                                                                                       <  2:1321511/61‑1 (MQ=255)
                                                                                                                                  ggcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg  <  2:386876/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                          
TGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGG  >  NC_000913/357727‑357995

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: