Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8160GA43.3% 110.3 / 132.3 134intergenic (+75/+144)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (73/3);  new base A (0/58);  total (73/61)
Fisher's exact test for biased strand distribution p-value = 8.00e-35
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

TTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622666‑1622966
                                                                                                                                                           |                                                                                                                                                        
ttttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCg                                                                                                                                                          >  1:359324/1‑151 (MQ=255)
  ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGa                                                                                                                                                         >  2:774049/1‑150 (MQ=255)
  ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGa                                                                                                                                                         <  2:1232249/150‑1 (MQ=255)
    aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaa                                                                                                                                                       <  2:730224/150‑1 (MQ=255)
    aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaa                                                                                                                                                       >  2:353963/1‑150 (MQ=255)
      ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaa                                                                                                                                                     >  4:6839/1‑150 (MQ=255)
      ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaa                                                                                                                                                     >  2:68771/1‑150 (MQ=255)
       gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaaaa                                                                                                                                                   >  1:343847/1‑151 (MQ=255)
          gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAcc                                                                                                                                              >  4:697405/1‑151 (MQ=255)
          gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAcc                                                                                                                                              >  4:573311/1‑151 (MQ=255)
          gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAcc                                                                                                                                              >  2:981903/1‑151 (MQ=255)
               gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAg                                                                                                                                         >  1:402083/1‑151 (MQ=255)
                   ataatTTTATTTTCTGGTTTTCGGGCCGGGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGACCCTTGTTACAGTTGTGTACAAGGGGAATGATTTTTATGACGGCGAAAAAAAACCGCCAGTaa                                                                                                                                      >  4:337079/1‑150 (MQ=255)
                           attttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGt                                                                                                                             >  2:997579/1‑151 (MQ=255)
                           attttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGt                                                                                                                             >  2:170838/1‑151 (MQ=255)
                             tttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                            >  4:1148652/1‑150 (MQ=255)
                              ttctgttttttgggCCGTTTATTTCAGCCAAAAGGGGGTTGACAAATGAAGGGAGATAACCCCCCCCCAACCCGCCCCTTTTTTACGGTTGTGTAAAGGGGGCTTGTTTTTTAGGCCGGCAAAAAAAACCCGCCAGTAAACCGGCGGTGa                                                                                                                           <  2:842255/150‑1 (MQ=255)
                               tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAAt                                                                                                                         >  2:964486/1‑151 (MQ=255)
                               tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGGAAt                                                                                                                         >  2:910887/1‑151 (MQ=255)
                                ttgtttttcgggCGGCTTTTTTCAGCCAAAAAGGGGATGACTAATGAGGGGAGATAACCCCCCCCCTACCCGCCCCTTTGTTACGGTGGTGTACAGGGGGCCTTTTTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAAt                                                                                                                         <  2:98463/146‑1 (MQ=255)
                                      ttctggCCGCTTTTTCCAGCCAAAAAGGGAATGACAAAGGAAGGGAAAAAACCCCTCCCCTAACCGCCCCTTTGTTACGGTGGTGAAAAGGGGGCTTGTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                                   <  4:550056/150‑1 (MQ=255)
                                      ttcgggCCGCTTTTTCCACCCAAAAGGGGGATGACTAAGGAGGGGAAATAACCCCCCCCCAACCCGCCCCTTTGTAACGGTTGTGTACAGGGGGCTTATTTTTTAGGAGGGGAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                                   <  4:120253/150‑1 (MQ=255)
                                      ttcggcccGCTTATTCCACCCAAAAAGGGGATGATTAATGAGCGGAGATAACCCCCCCCCTACCCGCCCCTTTGTTCCGGTTGTGAAAAGGGGGCTTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                                   <  4:1426373/144‑1 (MQ=255)
                                          gcccGCGTATTCCAGCCACAAGTGGGAGGAAAAAGGACGGAAAACAGCCCCTCTCCAATCCCCCCTCTTGTCACTGTTGTTAAAAAGGGTCCTGTTTTTGATGGCGGAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                               <  4:216147/148‑1 (MQ=255)
                                           ggcgcgTATTCCAGCAAGAAGGGGGTTAACTAAAGAAGGAGAATAAGCCCTCCCCTACCCGCCCCCTGGTTACTTTTGTTTAAAGGGGCCCTAATTTTGAGGCGGAAAAAAAAATGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                             <  2:1437036/149‑1 (MQ=255)
                                           cccGCGTTTTACACCCAAAAAGGGGATAACTAAGAAACGGAGAAACCCCCTCACCTACCCGGCCCTTTGTTACAGTTGTGAAAAGGGGGCTTGATTTTTGGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                             <  3:64256/150‑1 (MQ=255)
                                            ccGCTTATTTAGCCAAAAAGTGGGTTAATAAAAGAGCGGAGAAACCCCCTCCCCTACCCGCCCCTTTTTTAGATTGGTAAACAGGGGTCGTTTTTTTTTTCACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATa                                                                                                            <  4:854493/151‑1 (MQ=255)
                                                ttATTCCAGCCAAAAATGGGATGACAAATGAAGGGAGAAAACCCCTCCCCTACCCGCCCCCTTTTTACAGTTGTGTACAGGGGGCTTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                                        <  1:210178/150‑1 (MQ=255)
                                                gTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                                        >  2:1255087/1‑151 (MQ=255)
                                                  tttccAGCCAAAAAGGGGATGAATAATGAAGGGAAATAACCCCCCCCCTACCCGGCCCTTTGTTACGTTTGTGTACAAGGGGCCTGTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                                       <  2:1226022/146‑1 (MQ=255)
                                                   tttCACCCAAAAAGGGGATGACTAATGAACGGAAATAACCCCCCCCCAACCCGCCCCCTTGTTACATTTGTGTAAAGGGGGCTTTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtg                                                                                                      <  4:1400193/150‑1 (MQ=255)
                                                    ttcagatccaAATGGGATGACTAATGAACGGAGATAACCCCTCACCTAACCGGCCGCTTGTTATCGTTGTGTACAAGGGGCCTGATTTTTATGACGCCGAAAAAAGACCGCCAGTAAACCGGCGGAGAATGCTTGCATGGATAGATTggtg                                                                                                    >  4:177752/9‑151 (MQ=255)
                                                     tcaacaacaAATGGGATGACTAATGAGGGGAAATATCCCCCCCCCTACCCGGCCCTTTGTTACGGTTGTTTACAGGGGCCTTTTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTgtgt                                                                                                   <  3:86704/145‑1 (MQ=255)
                                                            caAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                            >  2:942717/1‑151 (MQ=255)
                                                            caAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                            >  1:316396/1‑151 (MQ=255)
                                                              aaagggggttaatAAAGAACGGAAAAAACCCCCCCCCCAACCCGCCCCTTTTTTACTTTTTTGTAAAGGGGGCTTTTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                            <  4:1115962/139‑1 (MQ=255)
                                                              aaTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGTAAGATTTGTGTTTTGCtttt                                                                                          >  4:269821/1‑151 (MQ=255)
                                                               aTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGGTTTGCtttt                                                                                          >  4:216585/1‑150 (MQ=255)
                                                                    gATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                                     >  2:474828/1‑151 (MQ=255)
                                                                         aaaaTGAAGGGAAAAAACCCCCCCCCTACCCGCCCCTTTTTTACGTTTGTGTAAAGGGGGCTTTTTTTTTATGAGGGGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAg                                                                                <  4:1179985/149‑1 (MQ=255)
                                                                            aTGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGc                                                                              >  2:714780/1‑150 (MQ=255)
                                                                              gAACGGAAATACCCCCTCCCCAACCCGCCCCCTTGTTACGTTTGTGTACAAGGGGCCTTTTTTTTAGGGCGGAAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtt                                                                           <  4:766095/151‑1 (MQ=255)
                                                                               aaCGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAttt                                                                          >  1:257964/1‑151 (MQ=255)
                                                                                cgggAGATAACCCCTCCCCTACCCGGCCCTTTTTTACAGTTGTGTAAAAGGGGCTTGTTTTTTATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt                                                                         <  2:269256/149‑1 (MQ=255)
                                                                                 cGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGGGTTTTGCTTTTACGCTAACAGGCATTTTc                                                                        >  2:1032955/1‑151 (MQ=255)
                                                                                   gagaTAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTc                                                                        >  3:450009/1‑149 (MQ=255)
                                                                                    gaaaTAACCCCCCCCCAACCCGCCCCTTTTTTACGGTTGTGTAAAGGGGGCTTGTTTTTATGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTcc                                                                       <  3:536328/147‑1 (MQ=255)
                                                                                        aactCCCCCCCCTACCCGCCCCTTTTTTACAGTTTTGAAAAAGGGGCTTGTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                                  <  3:327623/148‑1 (MQ=255)
                                                                                         aaTCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                                 >  4:1124353/1‑151 (MQ=255)
                                                                                           tCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGCATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGAATTTTCCTGCACTg                                                                >  4:1346397/1‑150 (MQ=255)
                                                                                           tCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGa                                                               >  3:178821/1‑151 (MQ=255)
                                                                                            ccccccccaacccGCCCCTTTTTTACATTTGTTAACAGGGGGCTTTTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                              <  4:980710/140‑1 (MQ=255)
                                                                                            cccTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                              >  2:1234594/1‑151 (MQ=255)
                                                                                            cccTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                              >  2:7973/1‑151 (MQ=255)
                                                                                             ccacACCAACCCGGCCTGTGGTTAAATTTGTTTACAGGGGCCTTTTTTTTTTTGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa                                                             <  4:581070/148‑1 (MQ=255)
                                                                                             ccTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                              >  4:1436817/1‑150 (MQ=255)
                                                                                               cccccTACCCGCCCCTTTTTTAATTTTGTGTAAAAGGGCCTTGTTTTTAATGGCGCAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                           <  2:426347/148‑1 (MQ=255)
                                                                                                ccccTACCCGCCCCTTTTTTACATTTGTGTACAGGGGGCTTGTTTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACg                                                          <  2:487063/149‑1 (MQ=255)
                                                                                                ccccTACCCGCCCCTTTGTTCCATTTTTGAAAAAGGGCCCTGTTTTGACGGGGGGAAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                           <  2:204665/148‑1 (MQ=255)
                                                                                                 aCCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACg                                                          >  4:1364682/1‑150 (MQ=255)
                                                                                                  ccTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                                         >  4:717273/1‑150 (MQ=255)
                                                                                                     acccGCCCCCTTTTTACATTTGTGAAAAGGGGCCTTGTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCg                                                     <  2:284655/149‑1 (MQ=255)
                                                                                                        cgcccctTTTTTACGGTTGTGTAAAGGGGGCTTGTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                  <  2:697408/148‑1 (MQ=255)
                                                                                                        cGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                  >  2:1205148/1‑151 (MQ=255)
                                                                                                          cccctTTTTTACAGTTGTGTACAGGGGGCTTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGac                                                <  2:274692/150‑1 (MQ=255)
                                                                                                          cccctTTGTTACAGTTGTGAAAAGGGGGCTTTTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGac                                                <  2:406074/150‑1 (MQ=255)
                                                                                                           ccccTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGac                                                >  4:678067/1‑150 (MQ=255)
                                                                                                             cctttGTTAGGGTGGTGTACAGGGGGCTTTTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                              <  4:769344/148‑1 (MQ=255)
                                                                                                                ttGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt                                           >  3:414531/1‑151 (MQ=255)
                                                                                                                 tGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                                          >  3:286866/1‑151 (MQ=255)
                                                                                                                   ttACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc                                        >  4:360050/1‑151 (MQ=255)
                                                                                                                     acatttgttaaaaggggcCCTGTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCa                                       <  2:1191047/136‑1 (MQ=255)
                                                                                                                     aCAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                                      >  2:547/1‑151 (MQ=255)
                                                                                                                      gaGTTGTATAAAGGGGGCCTGATTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                                     <  2:1032329/150‑1 (MQ=255)
                                                                                                                        gTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAg                                   >  3:540267/1‑151 (MQ=255)
                                                                                                                         ttGTGTAAAAGGGCCCGTTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                                      <  3:377062/147‑1 (MQ=255)
                                                                                                                           gtgtACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt                                >  3:208472/1‑151 (MQ=255)
                                                                                                                           gtgtAAAGGGGGCTTTTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt                                <  2:359620/151‑1 (MQ=255)
                                                                                                                           gtgtAAAAGGGGCTTTTTTTTTAGGAGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtt                                 <  4:528552/150‑1 (MQ=255)
                                                                                                                            tgtACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGAGAATGCTTGCATGGATAGGTTTGTGGTTTGCTTTTACGCTAACAGGCATTTTCCCGCACTGATAAAGAATCGTTGTCACAGTAGCATCAGtttt                               >  2:88920/1‑151 (MQ=255)
                                                                                                                               aCAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                            >  2:631093/1‑151 (MQ=255)
                                                                                                                                  gggcgCCTTATTTTTATGCGGACAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaat                         <  2:1149326/147‑1 (MQ=255)
                                                                                                                                   ggggCCTTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatg                        <  2:406651/151‑1 (MQ=255)
                                                                                                                                    gggCTTTATTTTTATGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTNNCaatga                       <  4:1259002/151‑1 (MQ=255)
                                                                                                                                    gggCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga                       >  1:299376/1‑151 (MQ=255)
                                                                                                                                     gggctttttttttATGCGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga                       <  4:484534/143‑1 (MQ=255)
                                                                                                                                        cTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGt                   >  2:715086/1‑151 (MQ=255)
                                                                                                                                        cTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAAACGGCGGTGAATGCTTTCATGGAAAGAATTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTTCACTGATAAAGAATCGTTGACACAGTAGCATCAATTTTTTCAATGAATGt                   >  4:1043680/1‑151 (MQ=255)
                                                                                                                                           gATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt                  >  2:1068172/1‑150 (MQ=255)
                                                                                                                                            ttttttATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTa                 <  2:6839/149‑1 (MQ=255)
                                                                                                                                            ttttttATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTa                 <  4:596864/149‑1 (MQ=255)
                                                                                                                                             tttttATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaa                <  2:573313/150‑1 (MQ=255)
                                                                                                                                              ttttATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAc              >  4:428955/1‑151 (MQ=255)
                                                                                                                                              ttttATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAc              >  2:727480/1‑151 (MQ=255)
                                                                                                                                              ttttATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaaa               <  4:1349718/150‑1 (MQ=255)
                                                                                                                                               tttATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACg             >  4:479765/1‑151 (MQ=255)
                                                                                                                                                 tATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg            >  2:340076/1‑150 (MQ=255)
                                                                                                                                                 aaTGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg            <  2:1367548/149‑1 (MQ=255)
                                                                                                                                                       gggCAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTa      <  4:1155498/150‑1 (MQ=255)
                                                                                                                                                           gAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACt  >  2:1175794/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                        
TTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAACT  >  NC_000913/1622666‑1622966

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