Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2750TC42.5% 94.8 / 14.3 73pseudogene (335/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base C (9/22);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 9.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCC  >  NC_000913/2228126‑2228419
                                                                                                                                                         |                                                                                                                                                
gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCt                                                                                                                                                  <  1:156122/150‑1 (MQ=255)
     ttttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCg                                                                                                                                               <  1:736339/148‑1 (MQ=255)
     ttttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGCG‑TTTCCCCGCCCCGG‑CTTT‑TAACGGACCCAAGCGGGGAGGAAACTcccccc                                                                                                                                            >  1:1012037/1‑151 (MQ=255)
      tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGC‑GGGTTTGCTCGCG‑GTTCGCCCCCCGTG‑CTTT‑TACCCTACCCCAGCGGGCAGTAAACTCcccc                                                                                                                                             >  2:448819/1‑149 (MQ=255)
      tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGTTCGCG‑TTCCGACGCCCGGG‑CTTTAAACCTAACCCAGG‑GGGCAGTAAACTCCGccc                                                                                                                                            >  2:359501/1‑150 (MQ=255)
       ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCg                                                                                                                                           <  1:1120649/150‑1 (MQ=255)
              tgctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCttt                                                                                                                                    <  1:960596/150‑1 (MQ=255)
                ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTTCCCCCGG‑TTCCCCGCCCCGGG‑TTTT‑AACC‑TGCCCCGGGGGGCAGTAAACCCCCCCCCAtttttt                                                                                                                                    >  2:1218362/1‑148 (MQ=255)
                ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCGCG‑TTCCGCCCCCCGTG‑CTTT‑TACCGTAACCAAGCCGGCAGTAAACCCCGCCCGATTCTTTtcc                                                                                                                                 >  1:528289/1‑151 (MQ=255)
                ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GCGTTTGTCTGCG‑GTTCCCCCCCCCTG‑CTTT‑TACCGCACCCACGCGGGGCGCAAACTCCCCCCCGTTCTTTtcc                                                                                                                                 >  2:602768/1‑151 (MQ=255)
                  gcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCa                                                                                                                               <  2:461332/151‑1 (MQ=255)
                   cAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCGCG‑GTCCGCCCCCCGTG‑CTTT‑AACCCCACCCCGGCGGGGAGTAACCCCCGCCCCGTTCTTTcccc                                                                                                                                >  1:1264067/1‑149 (MQ=255)
                    aaGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGCG‑TTCCGCCGCCCCTG‑CTTT‑AAACGTACCCCGGCGGGGAGTAAACTCCGCCCGGTTTTTTtccc                                                                                                                                >  1:911125/1‑148 (MQ=255)
                    aaGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCCCGCG‑TTTCCCCCCCCCTG‑CTTT‑AACCGTACCCCGGCGGGCAGTAACCTCCCCCCGATTCTTTtccc                                                                                                                                >  1:848893/1‑148 (MQ=255)
                       aCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGCG‑TTCCGCCCCCCCTG‑CTTT‑AAACCGACCCCGGCGGGGCGTAACCTCCGCCCGATTCTTTCCCCCGat                                                                                                                            >  1:416028/1‑149 (MQ=255)
                           tgtTTGCCGCCGCCTGGGTTTGCCACGCCCAGTACAGGCGGCCTCAGCTGCCAGCCATCAATTGGC‑CGGTTTGCTCCCG‑TTCCGCCGCCTGTG‑CTTG‑AACCGTACCCAGGCGGTCAGTAACCTCCGCCCGATTCTTTCCCCAGATATc                                                                                                                         <  2:971847/148‑1 (MQ=255)
                              ttAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGAt                                                                                                                    <  2:1018704/150‑1 (MQ=255)
                              ttAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑GTCCGCCGCCCCGG‑CTTT‑AAACGTACCCAGGCGGGGAGGAAACCCCCCCCCGGTCTTTCCCCAGa                                                                                                                             >  2:246648/1‑141 (MQ=255)
                                aaCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCCCG‑GTCCCCCCCCCCTGCCTTT‑AACCG‑ACCCCAGCGGGCCGGAAACCCCCCCCCCTTCTTTTCCCCAAAATCGAGAAGcc                                                                                                                 >  1:1052309/1‑151 (MQ=255)
                                   gccgccTGGGTTTGCCACTTCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGAAGTAACCTCCGCCCGATTCTTGCCCCAGATATCGAGATGCCTGc                                                                                                              <  2:312935/151‑1 (MQ=255)
                                   gccgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGc                                                                                                              <  2:114578/151‑1 (MQ=255)
                                    ccgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGcc                                                                                                             <  1:1301760/151‑1 (MQ=255)
                                     cgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                            <  1:245585/151‑1 (MQ=255)
                                     cgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                            <  2:493048/151‑1 (MQ=255)
                                      gccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccg                                                                                                           <  1:187696/151‑1 (MQ=255)
                                          gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccgccg                                                                                                        <  1:246705/150‑1 (MQ=255)
                                          gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccgcc                                                                                                         <  1:720260/149‑1 (MQ=255)
                                          gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccgcc                                                                                                         <  2:225974/149‑1 (MQ=255)
                                               tGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTTCTCGCGTTTCCCCCCCCCGT‑‑CTTT‑TACCGTACCCCGGCGGGCAGGAAACCCCCCCCCGTTCTTTTCCCCGATATCCCAGTGCCAGCCCGCCGCTaaac                                                                                                  >  2:952281/1‑149 (MQ=255)
                                                 cccttCCCAGTACGGGGGGGCTCCGCGCCCGGCCACCTAGGGGG‑GGGTTGGCCCGCG‑TCCCGCCGCCCGTG‑CTTT‑ACCCGTCCCCGGGCGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGATATGGGGATGCCAGCCCGCCGTTAAg                                                                                                   <  2:893423/144‑1 (MQ=255)
                                                 ccACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCT‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGAAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGc                                                                                                  <  1:776956/149‑1 (MQ=255)
                                                     tCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑GTCCGCCGCCCGTG‑CTTT‑AACCCTACCCAGGCGGGCAGTAACCTCCGCCCCATTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCaaaa                                                                                             >  1:414431/1‑150 (MQ=255)
                                                      cccAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTTCTCGCC‑TTCCCCCGCCCGTG‑CCTT‑TACCCTAACCCAGCGGGCCGGAACCCCCACCCGGTTTTTTTCCCCCATAACCAGATTCCCGACCCCCGTTAAGCCcaaa                                                                                             >  2:453063/1‑149 (MQ=255)
                                                      cccAGTACAGGCGGGCTACGCTGCCAGCCAGAAATTGGC‑GGGTTTGCTCGCG‑TCCCGCCGCCCGTG‑TTTT‑ACCCTTCCCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAg                                                                                            <  1:1355195/150‑1 (MQ=255)
                                                         aGTACAGTCGGGCTCCTGCGCTGGCAAGAAATTGGC‑GGTTTGGCGCGCT‑TCCCCCCCCCTGTA‑TTTC‑AGCCGCCCCCAGGGGGGAAGAACCCCCCCCCATATTCTTCCCCCAGATATGAAGATTCCGGCCCGCCGTTAAGCCAAAAGTAcc                                                                                        <  2:530685/151‑1 (MQ=255)
                                                           tACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGGTTGCTCGCG‑TTCCGCCCCCCCGG‑CCTT‑TACCCGACCCCGGGGGGGAGTAAACCCCGCCCCATTTTTTCCCCCGCAATCGCGATGGCCCCCCCCCCGTAAGCCcaaaagaac                                                                                        >  2:628169/1‑144 (MQ=255)
                                                              aGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTACGCCGCCCGTG‑CTTT‑AACCGTACCCCGGCGGGCAGGAACCCCCGCCCGCTTTCTTCCCCAGATACCGAGATGCCAACCCGCCGCTAAGCCacaa                                                                                             >  2:1074414/1‑141 (MQ=255)
                                                                  ggcctccGCTGCCACCCAGCAATTGGC‑GGTTTTGCTCGCG‑TCCCCCCCCCCTTG‑CTTT‑AACCGTCCCCAGGCGGGAAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTg                                                                                    <  2:1219172/140‑1 (MQ=255)
                                                                        cGCTGCCAGCCAGCAATTGGC‑GGTTTTGCTCGCG‑TTCCGCCGCCCGTG‑TTTT‑ACCCGTACCCAGGCGGGCAGTACCCCCCGCCCGATCCTTCCCCCAGATATCGAGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACg                                                                          <  1:260715/150‑1 (MQ=255)
                                                                           tGCCAGCCGGCATTGGGG‑GGGTTTGCTCGCG‑TGCCGCCGCCCGTG‑CTTT‑ACCCGTCCCCGGGCGGGAGGTACCCCCCGCCCGATTCTTTCCCCAGATGTCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCAc                                                                           <  2:528115/146‑1 (MQ=255)
                                                                             ccagccagCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                      <  1:1010925/149‑1 (MQ=255)
                                                                                  cagcAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTc                                                                  >  1:277285/1‑148 (MQ=255)
                                                                                     cAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAATACCGTTGGTGTACCACGGGCCGGGCGTAcc                                                             >  2:1098734/1‑150 (MQ=255)
                                                                                     acaTTGGC‑GGGTTTGTTCGCG‑TGGCGCCACACGTT‑CTTT‑AACCGTACCGGTGCGGCGAATTTTCTCCGCCCGATTCTGTCCCCTTGTCTCGGTTTCGCCGTCCGCCGTTAACCCAAAAGTACCGTTGGTGTACCGCGCGCCGGTCGTa                                                               <  2:501255/146‑1 (MQ=255)
                                                                                          ggC‑GGTTTTGCGCGCG‑TCCCCCCGCCGGTT‑CTTT‑ACCCGTCCCCGGGGGGGCGGTACCCCCCGCCCGTTTCTTCCCCCAGATGTCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTAcc                                                             <  2:1286528/145‑1 (MQ=255)
                                                                                            gggtgTTTGCCCGCG‑TCCCGCCGCCCGTG‑CTTT‑ACCCGCCCCGAGGCGGAAATTACCCCCCCCCCGATTTTTCCCCCAGTGCTCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGcc                                                      <  2:635796/147‑1 (MQ=255)
                                                                                                 ttcgatCGCG‑TCCCCCCGCCTGTG‑CTTA‑AACCGCCCCCAGGCGGAGATTACCCCCCCCCCGATTCTTCCCCAAGATACCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGAt                                                  <  1:834855/145‑1 (MQ=255)
                                                                                                   tccccgcg‑TTCCGCCGCCCGTG‑TTTT‑AACCGTACCCAGGCGGGCAGTACCCCCCCCCCGATTCTTTCCCCAGATATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCg                                                <  2:941200/146‑1 (MQ=255)
                                                                                                     cccgcg‑TCCCCCCGCCCGTC‑TTTT‑ACCGGTCCCCGGGGGGGAAGAACCCCCCGCCGGTTTTTTTCCCCAGAAACCGAGAGGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCg                                                <  2:849953/146‑1 (MQ=255)
                                                                                                      cgcccttcccgcccCCCTGG‑TTTT‑ACCCGTCCCCGGGGGGGCATAACCCCCCCCCCGTTTCTTCCCCCAGATTTCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTc                                             <  2:303800/143‑1 (MQ=255)
                                                                                                       cgcg‑TTCCCCCCCCCGTG‑TTTT‑AACCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAAATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTc                                             <  1:1379116/149‑1 (MQ=255)
                                                                                                        gcg‑TTCCGCCGCCCGTG‑TTTT‑ACCCGTCCCCAGGCGGGAAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCag                                           <  1:1020602/150‑1 (MQ=255)
                                                                                                        cgcttccccccccccGTG‑TTTT‑AACCCAACCCGGGCGGGCGGTCCCCTCCGCACGTTTTTTCCCCAAGAAATCGAGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTc                                             <  2:1050219/137‑1 (MQ=255)
                                                                                                         cg‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCag                                           >  2:611405/1‑149 (MQ=255)
                                                                                                                  cgccCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGc                                 >  1:590007/1‑151 (MQ=255)
                                                                                                                   gcacGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcc                                >  2:157714/4‑151 (MQ=255)
                                                                                                                    cccGTG‑CTTT‑ACCCTTACCCAGGCGGACAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGccc                               <  1:1074729/151‑1 (MQ=255)
                                                                                                                    cccGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGCTCAGAGCAAACGGccc                               >  1:1000803/1‑151 (MQ=255)
                                                                                                                     ccGTG‑TTTT‑ACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcccc                              <  1:8395/151‑1 (MQ=255)
                                                                                                                     ccGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcccc                              >  1:94052/1‑151 (MQ=255)
                                                                                                                      cGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcccc                              >  1:972345/1‑150 (MQ=255)
                                                                                                                        tG‑CTTT‑AACCGTCCCCGGGCGGGCAGTACCCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCa                             <  1:1159449/149‑1 (MQ=255)
                                                                                                                            ttt‑CACCGTCCCCAGGCGGGCAGTACCCCCCGCCCGTTTCTTTCCCCAGATATCGAGAGGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTaa                         <  2:129249/150‑1 (MQ=255)
                                                                                                                             tt‑ACCCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTa                          <  1:459641/148‑1 (MQ=255)
                                                                                                                             tt‑ACCCGTCCCCAGCCGGGCGGCACCCCCCGCCGGTTTCTTCCCCCAGATATCGATATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCa                             <  2:848638/145‑1 (MQ=255)
                                                                                                                                acccGCACCCAGGGGGGTCCTACCCCCCGCCCGATTCTTCCCCCAGATATGGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGTGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGc                       <  2:363028/147‑1 (MQ=255)
                                                                                                                                    gTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTg                 >  1:681692/1‑151 (MQ=255)
                                                                                                                                        ccAGGCGGGGAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGa             <  2:272450/151‑1 (MQ=255)
                                                                                                                                        ccAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCg              >  1:1379882/1‑150 (MQ=255)
                                                                                                                                          aGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGAc            >  2:326864/1‑150 (MQ=255)
                                                                                                                                            gggggCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAtt         <  1:952564/149‑1 (MQ=255)
                                                                                                                                              gggCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAtttt       >  1:62131/1‑151 (MQ=255)
                                                                                                                                              gggAAGTACCCCCCGCCCGTTTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACAttt        <  1:1009132/150‑1 (MQ=255)
                                                                                                                                                   gTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGcc  >  1:1323377/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                
GCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGC‑GGGTTTGCTCGCG‑TTCCGCCGCCCGTG‑CTTT‑AACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCC  >  NC_000913/2228126‑2228419

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