Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,995,4570TG41.9% 169.7 / 19.5 93G232G (GGT→GGGdapFdiaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base T (9/45);  new base G (39/0);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 6.86e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTG  >  NC_000913/3995311‑3995606
                                                                                                                                                  |                                                                                                                                                     
cGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGATCAGCAAGGGTTg                                                                                                                                                    >  4:64124/1‑150 (MQ=255)
 gTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCt                                                                                                                                                  <  4:780851/151‑1 (MQ=255)
 gTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGCGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGGAGGGAATCAGCAAAGttttc                                                                                                                                                   >  4:879600/1‑148 (MQ=255)
   tttCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGGAGGGGTTCCGCCAGGGTTTCTgg                                                                                                                                                >  4:1047777/1‑151 (MQ=255)
   tttCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTg                                                                                                                                                 <  1:1041077/150‑1 (MQ=255)
   tttCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTg                                                                                                                                                 <  4:836663/150‑1 (MQ=255)
     tCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTGCTgggc                                                                                                                                              >  1:1125487/1‑149 (MQ=255)
      ccGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCg                                                                                                                                             <  1:1086421/151‑1 (MQ=255)
           gcgcgcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGTTTCAGCAAGGTTTGCTGGCCgaagaa                                                                                                                                        <  1:1451369/151‑1 (MQ=255)
               gcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTAc                                                                                                                                    <  2:108561/151‑1 (MQ=255)
               gcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGAAAGGGTTTCTGGCCGgagaagaa                                                                                                                                     >  4:224995/1‑148 (MQ=255)
                 cAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGaagaagga                                                                                                                                     >  1:86762/1‑146 (MQ=255)
                 cAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTAcgc                                                                                                                                  <  1:102216/151‑1 (MQ=255)
                 cAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTAcgc                                                                                                                                  <  1:1182469/151‑1 (MQ=255)
                 cAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGATTAGCAAGGGTTGCTGGCCGGAGAAGGAcg                                                                                                                                   >  4:106742/1‑150 (MQ=255)
                 cAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGGTTGAGGAGGGGGTCAGCAAAGGTTGCttggcggaaaa                                                                                                                                        >  4:346450/1‑135 (MQ=255)
                      tCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGTAGGGGTTTAGGAAAGGTTGGTGGGCggagaaagaagcgggg                                                                                                                              >  1:790726/1‑140 (MQ=255)
                      tCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAAGGTTTCTGGCCGAAGAAATAAGCgggg                                                                                                                              >  1:912838/1‑150 (MQ=255)
                       cGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTTCTGGCCGGAGAAGGACGCgggg                                                                                                                              >  2:2347/1‑149 (MQ=255)
                       cGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGaa                                                                                                                            <  4:772451/151‑1 (MQ=255)
                       cGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTGCTGGCCGAAGAAGTACGCGTgg                                                                                                                              >  1:438846/1‑149 (MQ=255)
                          ttttATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGAGTTAGGAAAGGTTGCTGGCGGGAAAAGGAAGCGTTGAAATc                                                                                                                         >  1:298402/1‑151 (MQ=255)
                           tttATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGGAGGGGTTCAGCAAGGGTTGCTGGCCGAAGAAGAAAGCGGGGAACTcc                                                                                                                        >  4:364381/1‑151 (MQ=255)
                            ttATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTTCCGGAAGGAATCAGCAAAGGTTGCCGGGCGAAGAAAGACGCGGGGGAATccc                                                                                                                       >  4:349341/1‑151 (MQ=255)
                              aTGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTCAGCAAGGGTTGCTGGCCGGAGAAGTACGCGTGGGACTcccc                                                                                                                      >  1:1103749/1‑150 (MQ=255)
                              aTGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGATGCAGGAGGGAGTTCGCAAAGTGTGGCGGCCGAAAAAGAACGCGGGGGACTcccc                                                                                                                      >  4:123814/1‑150 (MQ=255)
                               tGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCcgg                                                                                                                    >  1:873678/1‑151 (MQ=255)
                                gCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCcggc                                                                                                                   <  4:1217874/151‑1 (MQ=255)
                                 cAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCGGTAGGGATTCAGCAAGGGTTTTTGGGCGGAGAAAAACCCGTgggaggcccccg                                                                                                                    >  1:122067/1‑139 (MQ=255)
                                  aaGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGTAGGGAGTCAGGAAGGGGTGATGGCCGGAGAAGTACGCGTGGGAGTCCCccgcgg                                                                                                                 >  1:111749/1‑151 (MQ=255)
                                  aaGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCCAGGGTTTCTGGCCGAAGAAGGACGCGTGGGAATCCCcggcgg                                                                                                                 >  1:49660/1‑151 (MQ=255)
                                    gTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGTAGGGATTCAGCAAAGGTTGCTGGCCGAAAAAAGACGCGGGGGACCCCCccgcgg                                                                                                                 >  4:800802/1‑149 (MQ=255)
                                     tGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGATGCAGGGGGGGGTCATCAAAGTTTTGTTGCCCCAGAAGGATGCGCGGGACCCCCCggcgggcg                                                                                                              >  4:1472717/1‑151 (MQ=255)
                                      ggTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGGAGGGATTCAGCCAGGGTTGCTGGCCGAAGAAGGACGCGTGGGACTCCCcggccg                                                                                                                 >  4:251700/1‑145 (MQ=255)
                                      ggTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTTCTGGCCGAAGAAGAACGCGGGGGAATCCCCCGCGggcg                                                                                                              >  4:1329334/1‑150 (MQ=255)
                                        ttAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTCAGCAAGGGTTGCTTGCCGCAGAAGAACGCGGGGGACTCCCCCGCGggcg                                                                                                              >  4:141921/1‑148 (MQ=255)
                                         tAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGGAGGGGGTCTGCAAGGGTTGGCTGCCGGAGAAGAAAGCGCGGGACTCCCccgcgggggcctt                                                                                                          >  4:202670/1‑144 (MQ=255)
                                          aaGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCtt                                                                                                          <  1:1233943/150‑1 (MQ=255)
                                          aaGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCtt                                                                                                          <  1:1407896/150‑1 (MQ=255)
                                          aaGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCtt                                                                                                          <  4:1261326/150‑1 (MQ=255)
                                            gcgcGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGTAGGGGTTTAGCCAGGGTTGCTGGCCGAAGAAGGACGCGGGGGACTCCCCggcgggcggttt                                                                                                          >  1:1085591/1‑144 (MQ=255)
                                             cgcgAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGata                                                                                                      <  1:906572/151‑1 (MQ=255)
                                             cgcgAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGata                                                                                                      <  1:422100/151‑1 (MQ=255)
                                              gcgAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGatat                                                                                                     <  1:1459185/151‑1 (MQ=255)
                                                  gCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGcc                                                                                                 <  4:1434943/151‑1 (MQ=255)
                                                    atatTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTCAGCAAGGGTTGCTGGGCGAAGAAGGAAGCGGGGAACTCCCCGGCGGGCGTCTTGATATCGcc                                                                                                 >  1:632845/1‑149 (MQ=255)
                                                     tatTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGAAGAAGTACGCGTGGGACTCCCCGGCGGGCGGCTTTATATCGCCTgg                                                                                              >  1:123291/1‑151 (MQ=255)
                                                      atTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTgg                                                                                              <  4:605015/150‑1 (MQ=255)
                                                         cGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAg                                                                                          <  4:121123/151‑1 (MQ=255)
                                                         cGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAg                                                                                          <  1:1195762/151‑1 (MQ=255)
                                                          gTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAAGGTTGCTGGCCGGAGAAGGACGCGGGGGAATCCCCGGCGGGCGTCTTGATATCGCCCGGgaa                                                                                           >  1:835212/1‑149 (MQ=255)
                                                            ttACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGTAGGGGATTAGCAAGGTTTTCTGGCCCGAGAAGGACGCGGGGGAATCCCCGGGGGGCGGCTTGTTATCGCCCGGGAAgg                                                                                         >  4:548569/1‑149 (MQ=255)
                                                            ttACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGt                                                                                        <  4:538238/150‑1 (MQ=255)
                                                            ttACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTCAGCAAGGGTTTCTGGCCGAAGAAGTACGCGGGGGACTCCCCCGCGGGCGTCTTGGTATCGCCCGGgaa                                                                                           >  1:281838/1‑147 (MQ=255)
                                                              aCGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGAAGAAAGACGCGTGGAACTCCCCGGCGGGCGGCTTGGTATCGCCTGGaaaag                                                                                         >  1:1309714/1‑145 (MQ=255)
                                                                  gTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTc                                                                                 <  4:1125397/151‑1 (MQ=255)
                                                                  gTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGGTGGAGGAGGGGATCAGCAAGGGTTGCTGGCCGGAGAAGGACGCGTGGGACTCCCCGGCGGGCGGCTTGATATCGCCCGGGAAAGGCCggg                                                                                   >  4:242385/1‑149 (MQ=255)
                                                                   tttATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGAATCAGGAAGGTTTGCTTGGCGAAGAAATACGCGGGGGACTCCCCCGCGGGCGGCTTTGTATCGCCCGGGAAGGGCCcggg                                                                                  >  1:362714/1‑148 (MQ=255)
                                                                    ttATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGTAGGAGGGGTTCAGCAAGGGTTTCTGGCCCGAGAAGGACGCGTGGGACTCCCCGGCGGGCGTCCTTTTATTGCCCGGGAAAGGCCggggga                                                                                >  4:938493/1‑147 (MQ=255)
                                                                     tATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCAc                                                                               <  1:1251103/150‑1 (MQ=255)
                                                                          gCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTta                                                                         <  4:209415/151‑1 (MQ=255)
                                                                            gTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGGTGCCGGAGGGGATCAGCCAGGGTTTCTGGCCGAAGAAGGACGCGGGGGGCCCCCCGGCGGTCGGCTTTTTATTTCCTTGAAAGGGCCGGGGCACCCggt                                                                          >  4:396941/1‑146 (MQ=255)
                                                                             tGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCCGCGGGCGGCTTTGTATCGCCTGGAAAAGGCCGGGGCACCCGTtat                                                                        >  1:1150762/1‑149 (MQ=255)
                                                                             tGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGGTGCAGGAGGGATTTAGCAAGGGTTGCTGGCCGAAGAAAGAAGCGGGGGACTCCCCCGCGGGCGGCTTTGTATTCCCTGGGAAAGGCCGGGGCACCCCTtat                                                                        >  4:952823/1‑149 (MQ=255)
                                                                              ggggCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTtatata                                                                     <  4:1390437/151‑1 (MQ=255)
                                                                                  cAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGAc                                                                 <  4:244345/151‑1 (MQ=255)
                                                                                      agaAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGAATTAGCAAGGGTTTCTGGCCGGAGAAGGACGCGGGGGACTCCCCCGCGGGCGGCTTTGTTTCCCCCGGGAAAGGCCGGGGCACCCCGTTTATATGAATgg                                                              >  4:142804/1‑150 (MQ=255)
                                                                                        aaaCCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGccc                                                           <  1:133717/151‑1 (MQ=255)
                                                                                        aaaCCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGccc                                                           <  3:15906/151‑1 (MQ=255)
                                                                                           cccAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCcggc                                                        <  1:813249/151‑1 (MQ=255)
                                                                                                  cTGCGGCAGCGGCGCGTGTGCGGCGGTTTCCGGAGGGATTCAGCCAGGGTTGCTGGGCGGAAAAAGAAGCGGGGGACCCCCCGGCGGGCGTTTTTTTATTGCCTGGGAAAGAGCGGGGGACCCCTTATATATGTATGGCCCCGGGGGACa                                                  >  4:777602/1‑150 (MQ=255)
                                                                                                   tGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGAGTCAGCAAGGGTTGCTGGCCGAAGAAGGACGCGGGGGACTCCCCCGCGGGCGGCTTGATATCGCCTGGGAAGGTCCGGGGCACCCCCTATATATGACTGGCCccgcggga                                                    >  4:387162/1‑145 (MQ=255)
                                                                                                    gCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATg                                                <  2:94396/150‑1 (MQ=255)
                                                                                                    gCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGt                                               <  4:557858/151‑1 (MQ=255)
                                                                                                        cAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTac                                           <  1:1204823/151‑1 (MQ=255)
                                                                                                         aGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTacg                                          <  1:1331946/151‑1 (MQ=255)
                                                                                                         aGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTCAGCAAGGGTTGCTGGCCGAAAAAGTAAGCGGGGGAATCCCCCGCGGGCGGCTTGATATCGCCTGGAAAAGGCCCGGGCACCCGGTATATATGAATGGGCCGGCGGGACa                                                  >  1:111955/1‑143 (MQ=255)
                                                                                                          gcggcgCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTacga                                         <  4:487084/151‑1 (MQ=255)
                                                                                                           cggcgCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTacgac                                        <  4:153068/151‑1 (MQ=255)
                                                                                                           cggcgCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTGCTGGCCGAAGAAGTACGCGGGGGACTCCCCCGCGGGCGGCTTGATATCGCCCGGGAAAGGCCGGGGGACCCCTTATATATGACCGGCCcggggg                                                      >  4:511115/1‑137 (MQ=255)
                                                                                                            ggcgCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTacgac                                        <  3:193191/150‑1 (MQ=255)
                                                                                                              cgcgTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGa                                     <  2:28542/151‑1 (MQ=255)
                                                                                                                   gtgCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTAt                                <  1:342995/151‑1 (MQ=255)
                                                                                                                        gcggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATc                           <  1:975537/151‑1 (MQ=255)
                                                                                                                         cggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTcct                            >  4:783530/1‑147 (MQ=255)
                                                                                                                         cggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCt                          <  4:801716/151‑1 (MQ=255)
                                                                                                                          ggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCt                          >  4:925068/1‑150 (MQ=255)
                                                                                                                           gTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTAt                        <  4:364651/151‑1 (MQ=255)
                                                                                                                              gCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGaa                     <  4:431587/151‑1 (MQ=255)
                                                                                                                              gCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGaa                     <  4:691754/151‑1 (MQ=255)
                                                                                                                                 gTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCa                  <  1:724223/151‑1 (MQ=255)
                                                                                                                                            gCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGga         >  2:177325/1‑149 (MQ=255)
                                                                                                                                              aaGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGgaa        >  4:257614/1‑148 (MQ=255)
                                                                                                                                                 gTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTg  <  2:85339/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTG  >  NC_000913/3995311‑3995606

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