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

GCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCC  >  NC_000913/2228126‑2228419
                                                                                                                                                         |                                                                                                                                                
gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCt                                                                                                                                                  <  1:156122/150‑1 (MQ=255)
     ttttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCg                                                                                                                                               <  1:736339/148‑1 (MQ=255)
     ttttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTTCCCCGCCCCGGCTTTTAACGGACCCAAGCGGGGAGGAAACTcccccc                                                                                                                                            >  1:1012037/1‑151 (MQ=255)
      tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGTTTGCTCGCGGTTCGCCCCCCGTGCTTTTACCCTACCCCAGCGGGCAGTAAACTCcccc                                                                                                                                             >  2:448819/1‑149 (MQ=255)
      tttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGTTCGCGTTCCGACGCCCGGGCTTTAAACCTAACCCAGGGGGCAGTAAACTCCGccc                                                                                                                                            >  2:359501/1‑150 (MQ=255)
       ttCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCg                                                                                                                                           <  1:1120649/150‑1 (MQ=255)
              tgctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCttt                                                                                                                                    <  1:960596/150‑1 (MQ=255)
                ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTTCCCCCGGTTCCCCGCCCCGGGTTTTAACCTGCCCCGGGGGGCAGTAAACCCCCCCCCAtttttt                                                                                                                                    >  2:1218362/1‑148 (MQ=255)
                ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCCCCCGTGCTTTTACCGTAACCAAGCCGGCAGTAAACCCCGCCCGATTCTTTtcc                                                                                                                                 >  1:528289/1‑151 (MQ=255)
                ctgcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGCGTTTGTCTGCGGTTCCCCCCCCCTGCTTTTACCGCACCCACGCGGGGCGCAAACTCCCCCCCGTTCTTTtcc                                                                                                                                 >  2:602768/1‑151 (MQ=255)
                  gcAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCa                                                                                                                               <  2:461332/151‑1 (MQ=255)
                   cAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGGTCCGCCCCCCGTGCTTTAACCCCACCCCGGCGGGGAGTAACCCCCGCCCCGTTCTTTcccc                                                                                                                                >  1:1264067/1‑149 (MQ=255)
                    aaGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCGCCCCTGCTTTAAACGTACCCCGGCGGGGAGTAAACTCCGCCCGGTTTTTTtccc                                                                                                                                >  1:911125/1‑148 (MQ=255)
                    aaGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCCCGCGTTTCCCCCCCCCTGCTTTAACCGTACCCCGGCGGGCAGTAACCTCCCCCCGATTCTTTtccc                                                                                                                                >  1:848893/1‑148 (MQ=255)
                       aCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTCCGCCCCCCCTGCTTTAAACCGACCCCGGCGGGGCGTAACCTCCGCCCGATTCTTTCCCCCGat                                                                                                                            >  1:416028/1‑149 (MQ=255)
                           tgtTTGCCGCCGCCTGGGTTTGCCACGCCCAGTACAGGCGGCCTCAGCTGCCAGCCATCAATTGGCCGGTTTGCTCCCGTTCCGCCGCCTGTGCTTGAACCGTACCCAGGCGGTCAGTAACCTCCGCCCGATTCTTTCCCCAGATATc                                                                                                                         <  2:971847/148‑1 (MQ=255)
                              ttAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGAt                                                                                                                    <  2:1018704/150‑1 (MQ=255)
                              ttAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCGCCGCCCCGGCTTTAAACGTACCCAGGCGGGGAGGAAACCCCCCCCCGGTCTTTCCCCAGa                                                                                                                             >  2:246648/1‑141 (MQ=255)
                                aaCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGGTCCCCCCCCCCTGCCTTTAACCGACCCCAGCGGGCCGGAAACCCCCCCCCCTTCTTTTCCCCAAAATCGAGAAGcc                                                                                                                 >  1:1052309/1‑151 (MQ=255)
                                   gccgccTGGGTTTGCCACTTCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGAAGTAACCTCCGCCCGATTCTTGCCCCAGATATCGAGATGCCTGc                                                                                                              <  2:312935/151‑1 (MQ=255)
                                   gccgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGc                                                                                                              <  2:114578/151‑1 (MQ=255)
                                    ccgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGcc                                                                                                             <  1:1301760/151‑1 (MQ=255)
                                     cgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                            <  1:245585/151‑1 (MQ=255)
                                     cgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGccc                                                                                                            <  2:493048/151‑1 (MQ=255)
                                      gccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccg                                                                                                           <  1:187696/151‑1 (MQ=255)
                                          gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccgccg                                                                                                        <  1:246705/150‑1 (MQ=255)
                                          gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccgcc                                                                                                         <  1:720260/149‑1 (MQ=255)
                                          gggTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCccgcc                                                                                                         <  2:225974/149‑1 (MQ=255)
                                               tGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGTTTCCCCCCCCCGTCTTTTACCGTACCCCGGCGGGCAGGAAACCCCCCCCCGTTCTTTTCCCCGATATCCCAGTGCCAGCCCGCCGCTaaac                                                                                                  >  2:952281/1‑149 (MQ=255)
                                                 cccttCCCAGTACGGGGGGGCTCCGCGCCCGGCCACCTAGGGGGGGGTTGGCCCGCGTCCCGCCGCCCGTGCTTTACCCGTCCCCGGGCGGGCAGTACCCCCCCCCCGATTCTTCCCCCAGATATGGGGATGCCAGCCCGCCGTTAAg                                                                                                   <  2:893423/144‑1 (MQ=255)
                                                 ccACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCTTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGAAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGc                                                                                                  <  1:776956/149‑1 (MQ=255)
                                                     tCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGGTCCGCCGCCCGTGCTTTAACCCTACCCAGGCGGGCAGTAACCTCCGCCCCATTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCaaaa                                                                                             >  1:414431/1‑150 (MQ=255)
                                                      cccAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCCTTCCCCCGCCCGTGCCTTTACCCTAACCCAGCGGGCCGGAACCCCCACCCGGTTTTTTTCCCCCATAACCAGATTCCCGACCCCCGTTAAGCCcaaa                                                                                             >  2:453063/1‑149 (MQ=255)
                                                      cccAGTACAGGCGGGCTACGCTGCCAGCCAGAAATTGGCGGGTTTGCTCGCGTCCCGCCGCCCGTGTTTTACCCTTCCCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAg                                                                                            <  1:1355195/150‑1 (MQ=255)
                                                         aGTACAGTCGGGCTCCTGCGCTGGCAAGAAATTGGCGGTTTGGCGCGCTTCCCCCCCCCTGTATTTCAGCCGCCCCCAGGGGGGAAGAACCCCCCCCCATATTCTTCCCCCAGATATGAAGATTCCGGCCCGCCGTTAAGCCAAAAGTAcc                                                                                        <  2:530685/151‑1 (MQ=255)
                                                           tACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCGCGTTCCGCCCCCCCGGCCTTTACCCGACCCCGGGGGGGAGTAAACCCCGCCCCATTTTTTCCCCCGCAATCGCGATGGCCCCCCCCCCGTAAGCCcaaaagaac                                                                                        >  2:628169/1‑144 (MQ=255)
                                                              aGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTACGCCGCCCGTGCTTTAACCGTACCCCGGCGGGCAGGAACCCCCGCCCGCTTTCTTCCCCAGATACCGAGATGCCAACCCGCCGCTAAGCCacaa                                                                                             >  2:1074414/1‑141 (MQ=255)
                                                                  ggcctccGCTGCCACCCAGCAATTGGCGGTTTTGCTCGCGTCCCCCCCCCCTTGCTTTAACCGTCCCCAGGCGGGAAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTg                                                                                    <  2:1219172/140‑1 (MQ=255)
                                                                        cGCTGCCAGCCAGCAATTGGCGGTTTTGCTCGCGTTCCGCCGCCCGTGTTTTACCCGTACCCAGGCGGGCAGTACCCCCCGCCCGATCCTTCCCCCAGATATCGAGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACg                                                                          <  1:260715/150‑1 (MQ=255)
                                                                           tGCCAGCCGGCATTGGGGGGGTTTGCTCGCGTGCCGCCGCCCGTGCTTTACCCGTCCCCGGGCGGGAGGTACCCCCCGCCCGATTCTTTCCCCAGATGTCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCAc                                                                           <  2:528115/146‑1 (MQ=255)
                                                                             ccagccagCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTACCCCCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGcc                                                                      <  1:1010925/149‑1 (MQ=255)
                                                                                  cagcAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTc                                                                  >  1:277285/1‑148 (MQ=255)
                                                                                     cAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCCGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAATACCGTTGGTGTACCACGGGCCGGGCGTAcc                                                             >  2:1098734/1‑150 (MQ=255)
                                                                                     acaTTGGCGGGTTTGTTCGCGTGGCGCCACACGTTCTTTAACCGTACCGGTGCGGCGAATTTTCTCCGCCCGATTCTGTCCCCTTGTCTCGGTTTCGCCGTCCGCCGTTAACCCAAAAGTACCGTTGGTGTACCGCGCGCCGGTCGTa                                                               <  2:501255/146‑1 (MQ=255)
                                                                                          ggCGGTTTTGCGCGCGTCCCCCCGCCGGTTCTTTACCCGTCCCCGGGGGGGCGGTACCCCCCGCCCGTTTCTTCCCCCAGATGTCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTAcc                                                             <  2:1286528/145‑1 (MQ=255)
                                                                                            gggtgTTTGCCCGCGTCCCGCCGCCCGTGCTTTACCCGCCCCGAGGCGGAAATTACCCCCCCCCCGATTTTTCCCCCAGTGCTCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGcc                                                      <  2:635796/147‑1 (MQ=255)
                                                                                                 ttcgatCGCGTCCCCCCGCCTGTGCTTAAACCGCCCCCAGGCGGAGATTACCCCCCCCCCGATTCTTCCCCAAGATACCGAGATCCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGAt                                                  <  1:834855/145‑1 (MQ=255)
                                                                                                   tccccgcgTTCCGCCGCCCGTGTTTTAACCGTACCCAGGCGGGCAGTACCCCCCCCCCGATTCTTTCCCCAGATATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCg                                                <  2:941200/146‑1 (MQ=255)
                                                                                                     cccgcgTCCCCCCGCCCGTCTTTTACCGGTCCCCGGGGGGGAAGAACCCCCCGCCGGTTTTTTTCCCCAGAAACCGAGAGGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCg                                                <  2:849953/146‑1 (MQ=255)
                                                                                                      cgcccttcccgcccCCCTGGTTTTACCCGTCCCCGGGGGGGCATAACCCCCCCCCCGTTTCTTCCCCCAGATTTCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTc                                             <  2:303800/143‑1 (MQ=255)
                                                                                                       cgcgTTCCCCCCCCCGTGTTTTAACCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAAATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTc                                             <  1:1379116/149‑1 (MQ=255)
                                                                                                        gcgTTCCGCCGCCCGTGTTTTACCCGTCCCCAGGCGGGAAGTACCCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCag                                           <  1:1020602/150‑1 (MQ=255)
                                                                                                        cgcttccccccccccGTGTTTTAACCCAACCCGGGCGGGCGGTCCCCTCCGCACGTTTTTTCCCCAAGAAATCGAGATGCCACCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTc                                             <  2:1050219/137‑1 (MQ=255)
                                                                                                         cgTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCag                                           >  2:611405/1‑149 (MQ=255)
                                                                                                                  cgccCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGc                                 >  1:590007/1‑151 (MQ=255)
                                                                                                                   gcacGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcc                                >  2:157714/4‑151 (MQ=255)
                                                                                                                    cccGTGCTTTACCCTTACCCAGGCGGACAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGAGGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGccc                               <  1:1074729/151‑1 (MQ=255)
                                                                                                                    cccGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGCTCAGAGCAAACGGccc                               >  1:1000803/1‑151 (MQ=255)
                                                                                                                     ccGTGTTTTACCCGTCCCCAGGCGGGCAGTACCCCCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcccc                              <  1:8395/151‑1 (MQ=255)
                                                                                                                     ccGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcccc                              >  1:94052/1‑151 (MQ=255)
                                                                                                                      cGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGcccc                              >  1:972345/1‑150 (MQ=255)
                                                                                                                        tGCTTTAACCGTCCCCGGGCGGGCAGTACCCCCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCa                             <  1:1159449/149‑1 (MQ=255)
                                                                                                                            tttCACCGTCCCCAGGCGGGCAGTACCCCCCGCCCGTTTCTTTCCCCAGATATCGAGAGGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTaa                         <  2:129249/150‑1 (MQ=255)
                                                                                                                             ttACCCGTCCCCAGGCGGGCAGTACCCTCCGCCCGATTCTTCCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTa                          <  1:459641/148‑1 (MQ=255)
                                                                                                                             ttACCCGTCCCCAGCCGGGCGGCACCCCCCGCCGGTTTCTTCCCCCAGATATCGATATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCa                             <  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)
                                                                                                                                                         |                                                                                                                                                
GCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCC  >  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: