Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3500CT26.8% 71.3 / 11.7 41S64F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (2/28);  new base T (11/0);  total (13/28)
Fisher's exact test for biased strand distribution p-value = 2.47e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTAC  >  NC_000913/4604117‑4604597
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
cAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                                                                                                                                                                                          <  1:130874/251‑1 (MQ=255)
cAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                                                                                                                                                                                          <  1:497708/251‑1 (MQ=255)
      aacttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGTAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                                                                                                                                                    <  2:214747/251‑1 (MQ=255)
      aacttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTTTGAAGTGGGCGCagc                                                                                                                                                                                                                                     >  1:226518/1‑249 (MQ=255)
       acttaacttGCCCTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGTATTAAGAAGAGAAATGTCTCACTTATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTGTATTCAGAAAGTGTGCGGGTAACTCTTATAAGCACCATTCCTATACCAGTCAGGAGTCATTTCAGGATTCCATGTCGCGGATATCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag‑a                                                                                                                                                                                                                                  <  2:531126/251‑1 (MQ=255)
         ttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTTCATGAGAAGTGAGTGCAGAG‑Agg                                                                                                                                                                                                                                >  1:269547/1‑251 (MQ=255)
           aacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTGCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGa                                                                                                                                                                                                                              <  1:190673/251‑1 (MQ=255)
                      tGTTTATGTAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATTAGGAGTGAGCGCAGAG‑Aggggttttgttgtcc                                                                                                                                                                                                                   >  1:242383/1‑239 (MQ=255)
                          tatgtagcagtcGCAGAAATGAACACAGCCGAATTAAGAAGAGAAATTTCGCTCTCATAGAAAAAGGCTTCATGAGTAGTAACGTTATTGATAGTGTATTAAGAAAGGTTGCGGGTAACTCTTATCAGCTCCATACCTATATCAGTCAGGAGTCATTTCAGGATTCCATGTCGCGGATTTAGTTTGCTGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAT‑AGGGATTATCTTGCctgac                                                                                                                                                                                                                <  2:136487/239‑1 (MQ=255)
                           aTGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCctgact                                                                                                                                                                                                               <  2:442626/250‑1 (MQ=255)
                                 aTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTTTTGAGAGTTTATTCAGAAAGTTTGCGTGTAACCCTTATAAGCTCCATACCTAGACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATCTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTgg                                                                                                                                                                                                        <  1:261435/251‑1 (MQ=255)
                                             ggAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCAGTTCAGGATGCTATGTCGCGGATCGCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑Gttt                                                                                                                                                                                            <  1:67236/250‑1 (MQ=255)
                                                      ccGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTAccc                                                                                                                                                                                  <  1:126225/251‑1 (MQ=255)
                                                       cGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTGTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCg                                                                                                                                                                                 <  2:226518/251‑1 (MQ=255)
                                                         aaTTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAAGATGCCATGTCGCGGATCTCGTTTGCGGCGGACATTTTTTTTTTTTTTGTCCTGAGAAGTGGGCGCAGAG‑AGGGATTATCTTTGCTTACTTAACTGGCGATTAA‑ATTTTcgcggtccccgc                                                                                                                                                                               >  1:435381/1‑243 (MQ=255)
                                                           ttAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAATTTTTTCTTTTTCTTCCCTGCGAAGGGGGCGCAGGG‑AGGGAATATCTTTGCTGACTGATCTGGCGGATAA‑GGTTTCGCGTAACCggggg                                                                                                                                                                             >  1:300250/1‑248 (MQ=255)
                                                                aagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTAATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGt                                                                                                                                                                        <  1:331485/251‑1 (MQ=255)
                                                                                 cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatga                                                                                                                                                         <  1:46207/249‑1 (MQ=255)
                                                                                  aTAGAAAAATGCGTCATGAGTAGTAGCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCGCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCAGATCTCGTTTGCGGAGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGTTTATCTTGCCTGACTTAACTGGCGATTAA‑TTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat                                                                                                                                                      <  2:528388/251‑1 (MQ=255)
                                                                                     gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGGAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAATGAGGGC‑GAGGAGGGATTTTGTTGCCTGACTTAAACGGGGGTTCA‑GTTTCCGCGTACACGCCGTTTATTTTTTTCGGATGGTGATATTgt                                                                                                                                                   >  1:314845/1‑250 (MQ=255)
                                                                                      aaaaaTGCGTCATTAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAATTTTTTTTTTTTTTCCCTGGGAAGTGAGGGGAGAG‑AGGGGTTTTTTTTCATGCCTGAACTGGCGATTAA‑GTTTCCCGGTCAACGGCGCTTAGTTTTTGCGGGTGATGATATTGa                                                                                                                                                   >  1:359723/1‑250 (MQ=255)
                                                                                          aTGCGTCAGGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTc                                                                                                                                              <  2:163740/251‑1 (MQ=255)
                                                                                          aTGCGTCAAGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGGTTTCGGGTAACCCTTATAAGCTCCATTCCTATACCAGTCAGGAGTTATTTCAGGATGCCATGTCGCGGAGCTCGTTTGCGGCTGTCATTTTTTCTTTTTCTGCCATGATAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTc                                                                                                                                              <  1:524409/251‑1 (MQ=255)
                                                                                                 aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTTCATGGGGAGGGAGCGCCGGG‑GGGGATTATCTTTCCTGACTGAACTGGCGATTAAAGTTTC‑GCGGACCCGCCGGTTTGTTATTGCGGGTGTTGATTTTTGAGGTCGGCTGAt                                                                                                                                       >  1:276755/1‑251 (MQ=255)
                                                                                                 aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAt                                                                                                                                       <  1:180890/251‑1 (MQ=255)
                                                                                                        gtagcGGTATTGAGATTGTAATCATAAAGAGTGCGGGTAACCCTTATAAGCTCCATACATATACCAGTCTGGAGTCATTTCAGGATGCCAGGTTGCGGATATCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGTGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTCCCCGGCGGTTAGTAATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa                                                                                                                                <  2:299606/247‑1 (MQ=255)
                                                                                                              gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGAGGAGTGAGGGCAGAG‑AGGGATTATTTTTCCCGTCTGAACTTGCGATTTA‑GGTTCCGCGTTCCCGGCGTTTAGGTATTGCGGGTGGTGGTATTGAAGCGCgggcggttggtttttttttg                                                                                                                          >  1:380968/1‑233 (MQ=255)
                                                                                                                      agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATTCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc                                                                                                                  <  2:49663/251‑1 (MQ=255)
                                                                                                                      agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc                                                                                                                  <  2:388871/251‑1 (MQ=255)
                                                                                                                        tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCATGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCATGAGAAGTGAGCGCAGAG‑AGGGGTTATCTTTTCTGACTTAACTGGCGATTAA‑GTTTCTGCGGACCCGGCGTTTTGTTTTTTTGGGTGATGATATTTTAGATCGGCTGTGTGGGTCATTTTCGCCAACAcccc                                                                                                                >  1:174792/1‑249 (MQ=255)
                                                                                                                                           ggTAAACCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTACATTGTCGCCATCACCGCTGGACGGTGTATTAAGTaa                                                                                             <  2:106007/251‑1 (MQ=255)
                                                                                                                                                    tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCCTGGGAAGGGAGCGCAGAG‑AGGGGTTTTTTTTCCTTACTGAACTGGCGGTTTA‑GTTTTCGCGGACCCCGCGGTTTTTTTTTTCGGGTTGTGGTATTTAAGCTCCGCTGGTTTGTTTATTTTCGCCATCCCCGCTGGGCGGTGGATTTAGTAAAGCGTccac                                                                                    >  1:64032/1‑248 (MQ=255)
                                                                                                                                                    tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCaa                                                                                     <  1:274699/250‑1 (MQ=255)
                                                                                                                                                     ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAc                                                                                    <  1:166902/250‑1 (MQ=255)
                                                                                                                                                                aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTGTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGAtttt                                                                        <  2:69811/251‑1 (MQ=255)
                                                                                                                                                                            cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAAc                                                             <  1:297935/250‑1 (MQ=255)
                                                                                                                                                                                              gCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCATGAGGAGGGTGCGCCGGG‑AGGGGTTATTTTGCCTGGCTGAACTGGGGGTTAA‑GTTTTCGCGTACCCGGCGGTTTGTTTTTGCGGGTGATGATATTTTAGCCTGGCTGATTGGTTCATTGTTGCCCTCACCGCTGGACGGGGGATTAAGTAAAACGTCAACGCTGGGGATTTTTCATTAAGGACTCTTTTTGTCATTTAAttg                                          >  1:78756/1‑249 (MQ=255)
                                                                                                                                                                                                cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGGGCGCAGGG‑AGGGGTTATCTTGCCTGACTGAAATGGCGGTTAA‑GTTTCCGCGTACCCGGCGGTTTGTTATTGCGGGTGATGATATTTAAGCTCGGCTGGTTGGGTTATTGTTGCCATCACCGCTGGGAGGGGGATTAAGTAAAACGGCAAAGCCGGGGATTTTTTTTTAGGAACACTTTTTTTCATTTAATgggg                                        >  1:110102/1‑249 (MQ=255)
                                                                                                                                                                                                              cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGTCTGAACTGGAGATTAA‑GTTGCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTGTTGTCATTAAATGGTGTACGTCAGGcgacc                          <  1:288784/251‑1 (MQ=255)
                                                                                                                                                                                                                   gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGGGAGCGCAGAG‑AGGGATTATCTTGCCTTACTGAACTTGCTATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACTCTGGAGATTTTGCATCAGGACCTCTTTTGGTCATTAAATGGTTTACGTCAGGccacccaccg                     <  2:154667/251‑1 (MQ=255)
                                                                                                                                                                                                                           cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAAc              <  1:458257/250‑1 (MQ=255)
                                                                                                                                                                                                                                      tttctGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGGGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTAc  >  1:455791/1‑251 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
CAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAG‑AGGGATTATCTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTAC  >  NC_000913/4604117‑4604597

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