Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2900AC43.3% 76.1 / 23.3 60intergenic (+96/+133)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base A (6/28);  new base C (26/0);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.45e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCT  >  NC_000913/3423142‑3423439
                                                                                                                                                    |                                                                                                                                                       
cgcATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCcccc                                                                                                                                                       >  2:191352/1‑151 (MQ=255)
  cATCCTCAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCt                                                                                                                                                     <  2:100712/151‑1 (MQ=255)
       tAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCTccccg                                                                                                                                                >  1:254478/1‑148 (MQ=255)
        aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTGCcccccccccccc                                                                                                                                                   >  2:244907/1‑147 (MQ=255)
        aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCcccccccccccc                                                                                                                                                   >  1:724831/1‑147 (MQ=255)
             agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTg                                                                                                                                          <  2:61830/151‑1 (MQ=255)
                   tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGTCTGCCCCCCCCCCCCCCGGCCGGTG‑Ct                                                                                                                                       >  2:349972/1‑147 (MQ=255)
                        gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGt                                                                                                                              <  2:834496/151‑1 (MQ=255)
                        gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCTCCCGGCCGCTG‑CTTCGCAACGt                                                                                                                              >  2:769973/1‑151 (MQ=255)
                              ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTG‑CTTCGCaaactttaaa                                                                                                                         >  2:261562/1‑142 (MQ=255)
                                  aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGc                                                                                                                    <  2:377805/151‑1 (MQ=255)
                                        ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTGGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCCCCcg                                                                                                               <  2:739103/150‑1 (MQ=255)
                                          tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCGGCCGTTG‑CTTCGCAAAATTCAAAAcccccccc                                                                                                                >  2:476923/1‑147 (MQ=255)
                                           aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                            <  2:254449/150‑1 (MQ=255)
                                             tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCTCCCGGCCGTTG‑CTTCGCAACATTCAAAACCGCTCCccgcgg                                                                                                           >  1:118189/1‑149 (MQ=21)
                                                aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                       <  2:34346/150‑1 (MQ=255)
                                                aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTG‑CTTCCCCACGATCAAATTCGCCCCCCGCGGGttt                                                                                                       >  1:515170/1‑150 (MQ=255)
                                                   aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCCGGCCCTTG‑CTTTGCAAAGTTCAAAAACGCCCCccgcggggtt                                                                                                       >  2:54042/1‑143 (MQ=18)
                                                    aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGACCCGTCCTCCCCCCTCCCGGCCGCTG‑TTTTTCAACGATAAAAACCCCTCCGCGCGGGttt                                                                                                       >  2:886541/1‑146 (MQ=255)
                                                    aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCGGGCCGGTG‑CTTCGCAAAGTTCAAAACCGCCCCCCGCGGGttt                                                                                                       >  1:690532/1‑146 (MQ=255)
                                                      tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGGCTTCCCCCCCCCCTCCCGGCCGGTT‑CTTCTCTACGTTTCAAACCCCCCCCCGCGGGGTTGTCCttc                                                                                                >  2:648490/1‑149 (MQ=255)
                                                      tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCCTTG‑CTTCGCAACATTCAAAACCCCCCCCCGCGGGTTTTTTCTAc                                                                                                >  1:394699/1‑151 (MQ=18)
                                                       gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                <  2:108776/150‑1 (MQ=255)
                                                        aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                              <  1:760000/151‑1 (MQ=255)
                                                        aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCCGTG‑CTTTGCAACGTTTCAATCCGCCCCCCGCGGGCTTGTCCTACTc                                                                                              >  2:432763/1‑151 (MQ=255)
                                                          tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCTCCGGGGCGGTT‑CTTTCCAACATTCAAATCCCCTCCCGGCGGGTTTTTCCttc                                                                                                >  1:412681/1‑145 (MQ=255)
                                                                 tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCAGCCCTCCCCCCCCCCCCCCGGCCGGTT‑CTTCGCCACACTCTAAATCGCCCCCCGCGGAGTTGTTCTACTCCGGcga‑ggg                                                                                    >  2:822593/1‑149 (MQ=255)
                                                                  cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGGCCCCCCCCCCCCCGGCCGCTG‑CTTTGCAACCATCAAATACCCCCCCCGCGGGTTTGTTCTACTCCGGGGA‑gtgg                                                                                   >  1:321549/1‑150 (MQ=255)
                                                                  cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCGCCCCCCACCCCCCCGGGCGGTG‑CTTCGTAACGGTCGAATCCCCCCCCCGCGGATTTGTCCTACTCggggg                                                                                         >  1:852688/1‑144 (MQ=255)
                                                                  cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGGCGTTG‑CTTCGCAAAGATTAAATTCGCCCCCCGCGGGTTTGTTCCACTCggg                                                                                           >  1:424813/1‑144 (MQ=17)
                                                                    aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTc                                                                                 <  1:157722/151‑1 (MQ=255)
                                                                    aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGGTG‑CTTTGCCACCCTCCAATCCCCCCCCCGCGGGGTTTTCCTCCTCcggggggggtt                                                                                   >  2:771806/1‑142 (MQ=12)
                                                                               ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGTTG‑CTTTTCAACGTTCAAATCCGCCCCCCGCGGGGTTTTCCTACTCGGGGGA‑GGGGTCACCGacaaa                                                                        >  1:871492/1‑149 (MQ=255)
                                                                                   aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGTTG‑CTTTTCAAACTTCAAATCCGCCCCCGGCGGGTTTGTCCTACTCCGGGGA‑GTGTTCACCCACAAacaaca                                                                   >  1:54210/1‑150 (MQ=17)
                                                                                      cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGTTT‑TCTCCCCACCTTTCAATCCGCTCCCGGCGGGTTTTGCCTTCTCCCGCGA‑GTCTTCAccccca                                                                          >  2:483703/1‑136 (MQ=9)
                                                                                          cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCCGGGCGGTG‑CTTCGCAAAGATCAAATCCCCTCCCGGCGGAGTTGTTCTTCTCCGGGGA‑GGGTTTCCCCACAAACAACCGATaaaa                                                            >  1:462834/1‑150 (MQ=14)
                                                                                                gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCGGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAgg                                                     <  1:886622/151‑1 (MQ=255)
                                                                                                 ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTCCCCCCACCCCCCCGGCCCGTG‑CTTCGCAACGTTCAAATCCGCCCCCCGCGGGTTTTTTCTACTCCGGAGA‑GGGTTTACCGACAAAAAAAAGATAaaaaaaaa                                                       >  1:447981/1‑148 (MQ=14)
                                                                                                   tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGGTT‑CTTCGCAACCGTCAAAACCCCCCCCCGCGGGGTTGTCCTACTCGGGGGA‑GTGGTCACCCACCAAAAACAAATAaaaaaaaaagacc                                                  >  1:91497/1‑146 (MQ=9)
                                                                                                    ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAGGGcc                                                   <  2:167035/149‑1 (MQ=255)
                                                                                                      cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                <  1:318029/150‑1 (MQ=255)
                                                                                                       gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                                              <  1:404126/151‑1 (MQ=40)
                                                                                                       gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCGGGCCGTTGCCTTTGAAAGTTTCAAA‑CCGCTCCCGGGGGGTTTTTCCTACCCTGGGGG‑GGGGTCTCCTACCAACAAAAAATAaaaaaaaaaagcccggc                                              >  2:105935/1‑142 (MQ=2)
                                                                                                        cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCCCCCCCCTCCGGGGCGGTG‑CTTTGCAACCTTCAAATTCGCCCCCCGCGGGTTTGTTCTACTCCCGAGA‑GTGCTCCCCCACAAAAAACAGAAAAAACAAAAAGccc                                                  >  1:104760/1‑146 (MQ=9)
                                                                                                          tgatttgatGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCGGGCGTTG‑CTTTGCCAAGTTCCAATCCGCCCCCCGCCGGTTTTTTCTCCCCGGGAGG‑GGGTTTTCCCCCAAAAAACCGAAAAAAAAAAAAGcccccgctt                                            >  2:240972/1‑144 (MQ=255)
                                                                                                            atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                         <  2:52911/151‑1 (MQ=38)
                                                                                                             tttgatGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCACCCTTCCGGCCGGTG‑CTTCGCAAAATTCTAAACCCCCCCCCGCGGGGTTTTTCTACTCGGGGGAGGGGTTGACC‑ACAAAAAACAaaaaaaaaaaaaaacaccagc                                              >  2:835463/1‑136 (MQ=2)
                                                                                                                       ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                 <  1:403449/148‑1 (MQ=35)
                                                                                                                       ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                <  2:515115/149‑1 (MQ=35)
                                                                                                                              tATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                        <  1:660368/150‑1 (MQ=38)
                                                                                                                              tATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                        <  1:19024/150‑1 (MQ=38)
                                                                                                                               aTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                       <  2:48662/150‑1 (MQ=38)
                                                                                                                                 ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTTCTCGGGAGA‑GTGTTCACCGGCAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                    <  1:270736/151‑1 (MQ=255)
                                                                                                                                 ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                    <  2:690456/151‑1 (MQ=38)
                                                                                                                                    gggcgTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                 <  1:206561/151‑1 (MQ=38)
                                                                                                                                       cgTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGt               <  2:462784/150‑1 (MQ=39)
                                                                                                                                        gTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTc              >  1:648543/1‑150 (MQ=40)
                                                                                                                                          ccTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg            <  1:476974/150‑1 (MQ=40)
                                                                                                                                             gCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc          >  1:334155/1‑149 (MQ=40)
                                                                                                                                             gCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc          >  2:421651/1‑149 (MQ=40)
                                                                                                                                             gCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa         <  2:58458/150‑1 (MQ=40)
                                                                                                                                                  ccACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc   <  1:2275/151‑1 (MQ=40)
                                                                                                                                                   cACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt  <  1:350015/151‑1 (MQ=40)
                                                                                                                                                    |                                                                                                                                                       
CGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGA‑GTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCT  >  NC_000913/3423142‑3423439

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