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

CGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCT  >  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)
                   tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGTCTGCCCCCCCCCCCCCCGGCCGGTGCt                                                                                                                                       >  2:349972/1‑147 (MQ=255)
                        gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                                                              <  2:834496/151‑1 (MQ=255)
                        gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCTCCCGGCCGCTGCTTCGCAACGt                                                                                                                              >  2:769973/1‑151 (MQ=255)
                              ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGCaaactttaaa                                                                                                                         >  2:261562/1‑142 (MQ=255)
                                  aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                    <  2:377805/151‑1 (MQ=255)
                                        ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTGGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCCCCcg                                                                                                               <  2:739103/150‑1 (MQ=255)
                                          tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCGGCCGTTGCTTCGCAAAATTCAAAAcccccccc                                                                                                                >  2:476923/1‑147 (MQ=255)
                                           aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                            <  2:254449/150‑1 (MQ=255)
                                             tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCTCCCGGCCGTTGCTTCGCAACATTCAAAACCGCTCCccgcgg                                                                                                           >  1:118189/1‑149 (MQ=21)
                                                aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                       <  2:34346/150‑1 (MQ=255)
                                                aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCCCCACGATCAAATTCGCCCCCCGCGGGttt                                                                                                       >  1:515170/1‑150 (MQ=255)
                                                   aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCCGGCCCTTGCTTTGCAAAGTTCAAAAACGCCCCccgcggggtt                                                                                                       >  2:54042/1‑143 (MQ=18)
                                                    aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGACCCGTCCTCCCCCCTCCCGGCCGCTGTTTTTCAACGATAAAAACCCCTCCGCGCGGGttt                                                                                                       >  2:886541/1‑146 (MQ=255)
                                                    aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCGGGCCGGTGCTTCGCAAAGTTCAAAACCGCCCCCCGCGGGttt                                                                                                       >  1:690532/1‑146 (MQ=255)
                                                      tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGGCTTCCCCCCCCCCTCCCGGCCGGTTCTTCTCTACGTTTCAAACCCCCCCCCGCGGGGTTGTCCttc                                                                                                >  2:648490/1‑149 (MQ=255)
                                                      tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCCTTGCTTCGCAACATTCAAAACCCCCCCCCGCGGGTTTTTTCTAc                                                                                                >  1:394699/1‑151 (MQ=18)
                                                       gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                <  2:108776/150‑1 (MQ=255)
                                                        aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                              <  1:760000/151‑1 (MQ=255)
                                                        aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCCGTGCTTTGCAACGTTTCAATCCGCCCCCCGCGGGCTTGTCCTACTc                                                                                              >  2:432763/1‑151 (MQ=255)
                                                          tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCTCCGGGGCGGTTCTTTCCAACATTCAAATCCCCTCCCGGCGGGTTTTTCCttc                                                                                                >  1:412681/1‑145 (MQ=255)
                                                                 tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCAGCCCTCCCCCCCCCCCCCCGGCCGGTTCTTCGCCACACTCTAAATCGCCCCCCGCGGAGTTGTTCTACTCCGGcgaggg                                                                                    >  2:822593/1‑149 (MQ=255)
                                                                  cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGGCCCCCCCCCCCCCGGCCGCTGCTTTGCAACCATCAAATACCCCCCCCGCGGGTTTGTTCTACTCCGGGGAgtgg                                                                                   >  1:321549/1‑150 (MQ=255)
                                                                  cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCGCCCCCCACCCCCCCGGGCGGTGCTTCGTAACGGTCGAATCCCCCCCCCGCGGATTTGTCCTACTCggggg                                                                                         >  1:852688/1‑144 (MQ=255)
                                                                  cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGGCGTTGCTTCGCAAAGATTAAATTCGCCCCCCGCGGGTTTGTTCCACTCggg                                                                                           >  1:424813/1‑144 (MQ=17)
                                                                    aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTc                                                                                 <  1:157722/151‑1 (MQ=255)
                                                                    aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGGTGCTTTGCCACCCTCCAATCCCCCCCCCGCGGGGTTTTCCTCCTCcggggggggtt                                                                                   >  2:771806/1‑142 (MQ=12)
                                                                               ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGTTGCTTTTCAACGTTCAAATCCGCCCCCCGCGGGGTTTTCCTACTCGGGGGAGGGGTCACCGacaaa                                                                        >  1:871492/1‑149 (MQ=255)
                                                                                   aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGTTGCTTTTCAAACTTCAAATCCGCCCCCGGCGGGTTTGTCCTACTCCGGGGAGTGTTCACCCACAAacaaca                                                                   >  1:54210/1‑150 (MQ=17)
                                                                                      cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGTTTTCTCCCCACCTTTCAATCCGCTCCCGGCGGGTTTTGCCTTCTCCCGCGAGTCTTCAccccca                                                                          >  2:483703/1‑136 (MQ=9)
                                                                                          cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCCGGGCGGTGCTTCGCAAAGATCAAATCCCCTCCCGGCGGAGTTGTTCTTCTCCGGGGAGGGTTTCCCCACAAACAACCGATaaaa                                                            >  1:462834/1‑150 (MQ=14)
                                                                                                gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCGGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAgg                                                     <  1:886622/151‑1 (MQ=255)
                                                                                                 ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCTCCCCCCACCCCCCCGGCCCGTGCTTCGCAACGTTCAAATCCGCCCCCCGCGGGTTTTTTCTACTCCGGAGAGGGTTTACCGACAAAAAAAAGATAaaaaaaaa                                                       >  1:447981/1‑148 (MQ=14)
                                                                                                   tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGGTTCTTCGCAACCGTCAAAACCCCCCCCCGCGGGGTTGTCCTACTCGGGGGAGTGGTCACCCACCAAAAACAAATAaaaaaaaaagacc                                                  >  1:91497/1‑146 (MQ=9)
                                                                                                    ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAGGGcc                                                   <  2:167035/149‑1 (MQ=255)
                                                                                                      cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                <  1:318029/150‑1 (MQ=255)
                                                                                                       gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                                              <  1:404126/151‑1 (MQ=40)
                                                                                                       gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCGGGCCGTTGCCTTTGAAAGTTTCAAACCGCTCCCGGGGGGTTTTTCCTACCCTGGGGGGGGGTCTCCTACCAACAAAAAATAaaaaaaaaaagcccggc                                              >  2:105935/1‑142 (MQ=2)
                                                                                                        cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTCCCCCCCCCCCTCCGGGGCGGTGCTTTGCAACCTTCAAATTCGCCCCCCGCGGGTTTGTTCTACTCCCGAGAGTGCTCCCCCACAAAAAACAGAAAAAACAAAAAGccc                                                  >  1:104760/1‑146 (MQ=9)
                                                                                                          tgatttgatGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCGGGCGTTGCTTTGCCAAGTTCCAATCCGCCCCCCGCCGGTTTTTTCTCCCCGGGAGGGGGTTTTCCCCCAAAAAACCGAAAAAAAAAAAAGcccccgctt                                            >  2:240972/1‑144 (MQ=255)
                                                                                                            atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                         <  2:52911/151‑1 (MQ=38)
                                                                                                             tttgatGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCACCCTTCCGGCCGGTGCTTCGCAAAATTCTAAACCCCCCCCCGCGGGGTTTTTCTACTCGGGGGAGGGGTTGACCACAAAAAACAaaaaaaaaaaaaaacaccagc                                              >  2:835463/1‑136 (MQ=2)
                                                                                                                       ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                 <  1:403449/148‑1 (MQ=35)
                                                                                                                       ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                <  2:515115/149‑1 (MQ=35)
                                                                                                                              tATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                        <  1:660368/150‑1 (MQ=38)
                                                                                                                              tATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                        <  1:19024/150‑1 (MQ=38)
                                                                                                                               aTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                       <  2:48662/150‑1 (MQ=38)
                                                                                                                                 ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTTCTCGGGAGAGTGTTCACCGGCAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                    <  1:270736/151‑1 (MQ=255)
                                                                                                                                 ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                    <  2:690456/151‑1 (MQ=38)
                                                                                                                                    gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                 <  1:206561/151‑1 (MQ=38)
                                                                                                                                       cgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGt               <  2:462784/150‑1 (MQ=39)
                                                                                                                                        gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTc              >  1:648543/1‑150 (MQ=40)
                                                                                                                                          ccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg            <  1:476974/150‑1 (MQ=40)
                                                                                                                                             gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc          >  1:334155/1‑149 (MQ=40)
                                                                                                                                             gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc          >  2:421651/1‑149 (MQ=40)
                                                                                                                                             gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa         <  2:58458/150‑1 (MQ=40)
                                                                                                                                                  ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc   <  1:2275/151‑1 (MQ=40)
                                                                                                                                                   cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt  <  1:350015/151‑1 (MQ=40)
                                                                                                                                                    |                                                                                                                                                       
CGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCT  >  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: