Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,245,2470GT40.3% 95.8 / 30.7 68E30* (GAA→TAA) ycgYuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (10/30);  new base T (26/1);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 2.06e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTC  >  NC_000913/1245102‑1245397
                                                                                                                                                  |                                                                                                                                                      
tCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTgtaa                                                                                                                                                     >  1:483046/1‑146 (MQ=255)
tCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGTACa                                                                                                                                                    >  2:937607/1‑150 (MQ=255)
     tCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAg                                                                                                                                               <  2:657243/150‑1 (MQ=255)
       ccGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAAt                                                                                                                                            <  2:1076253/151‑1 (MQ=255)
        cGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAAAAttt                                                                                                                                                 >  2:949690/1‑145 (MQ=255)
            aacaacaacAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTAGCTGCGCTTATTTAACATTTAGAAtt                                                                                                                                           >  2:683190/1‑147 (MQ=255)
             acaacaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCa                                                                                                                                      <  2:181370/151‑1 (MQ=255)
              caacaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAAAATTTAGAATTTATCa                                                                                                                                      >  2:719174/1‑150 (MQ=255)
                acaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTTAAaatttataatttattca                                                                                                                                     >  1:902072/1‑146 (MQ=255)
                 caacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATcc                                                                                                                                  <  2:279357/151‑1 (MQ=255)
                  aacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTTAAATTTTGAAATTATTAAtac                                                                                                                                  >  2:394007/1‑148 (MQ=255)
                   acaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCCCTTTTTTAACAATTTTAATTTATTAATCCtt                                                                                                                                >  2:418721/1‑151 (MQ=255)
                   acaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCtt                                                                                                                                <  2:30578/151‑1 (MQ=255)
                    caGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAACATTTATAATTTATTAAACCtt                                                                                                                                >  2:773532/1‑150 (MQ=255)
                    caGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTTAACATTTAGAATTGATCAATCCTTc                                                                                                                               <  1:898611/151‑1 (MQ=255)
                      gAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAAAATTTAGAATTTATTAATCCTTCt                                                                                                                              >  1:387765/1‑150 (MQ=255)
                      gAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTg                                                                                                                             <  2:38288/151‑1 (MQ=255)
                        aTCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTGACATTTAGAATTTGTCAATCCTTCttg                                                                                                                            >  2:452582/1‑148 (MQ=255)
                        aTCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGc                                                                                                                           <  1:1215578/151‑1 (MQ=255)
                          ccTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCa                                                                                                                          <  1:905407/150‑1 (MQ=255)
                               tttAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGGTTTTGGCTGGCTTTTTTAAAATTTTGAATTTATCAATCCTTCTGGCAGAtt                                                                                                                      >  2:298251/1‑149 (MQ=255)
                                  aaTTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGCGCTTATTTAACATTTAGAAATGATCAATCCTTCTTGCAGATTAAGcc                                                                                                                 >  2:650718/1‑151 (MQ=255)
                                   attattGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCa                                                                                                                <  1:830347/151‑1 (MQ=255)
                                       ttGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTTAAAATTTTAAATTGATCAATACTTTTGGGAGATTTAACCaaaa                                                                                                             >  1:1150598/1‑150 (MQ=255)
                                         gTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATAc                                                                                                          <  2:1091225/151‑1 (MQ=255)
                                            tCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTTAAAAATTAGAATTTATTAAAACTTCTGGCAGATTAAGCCaaaaaaacc                                                                                                        >  1:683095/1‑145 (MQ=255)
                                              gttgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGtt                                                                                                     <  1:1236096/151‑1 (MQ=255)
                                                tgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGATTTTGACGAGCTTTTTATAAAATTAAAAATGATCAATTCTTATTGCGAAATTAGCTAAAATaatcca                                                                                                      >  1:178773/1‑143 (MQ=255)
                                                   ttttGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGGGGTTATTTAAAAATTTAAATTTATAAAATCTTCTGGCAGATTTAAACAAAAAAcccccttgg                                                                                                  >  2:982839/1‑143 (MQ=255)
                                                   ttttGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGc                                                                                                 <  2:325537/150‑1 (MQ=255)
                                                    tttGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTAACATTTTGAATTTATCAATCCTTTTTGCAGATTAAGCCAAAAAAACCCTTTGGCa                                                                                                >  1:431017/1‑150 (MQ=255)
                                                       ggTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGCGGTTTTTGAAAAATTAGAATTGATAAATCCCTTTGGCAGATTAAGCCAAAAAAACCCTTTGGCaa                                                                                               >  2:975938/1‑148 (MQ=255)
                                                        gtgtGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAAAATTTAGAAATTATCAATCCTTCTTGCAGAATAAGCCAAAATACCCCTTCGGCAATggt                                                                                           >  1:896474/1‑149 (MQ=255)
                                                           tgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGGTTTGGCTGCGCTTTTTTAAAATTTAGAATTGATCAATCCTTTTTGCAGATTAAGCCAAAAAACCCGGTTGGCaaaggtt                                                                                          >  1:228988/1‑144 (MQ=255)
                                                               tCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCTGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGtt                                                                                     <  1:937810/150‑1 (MQ=255)
                                                                 aGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGATTTTGGTGCGCCTTTTTGACAATTTGAATTTGTTAAACCTTCTTGCAGATTTAGACAAAATACCCGCTCGGGAATTCttttgctac                                                                                   >  2:906430/1‑143 (MQ=255)
                                                                  gAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTAcgc                                                                                 <  2:762279/151‑1 (MQ=255)
                                                                   aaaaTAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGGGCTTAATGTACAATTAGAAATGTATAATCCTTCTTGCAGATTAAGCCAACAAAACCCTTTGGCAAAACTTCACTTAagc                                                                                 >  2:715317/1‑150 (MQ=255)
                                                                     aataaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGaa                                                                              <  2:915842/151‑1 (MQ=255)
                                                                      ataaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAag                                                                             <  2:796067/151‑1 (MQ=255)
                                                                      ataaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAag                                                                             <  1:967697/151‑1 (MQ=255)
                                                                       taaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGCGCGCTTATTTAAAATTTAGAAATGATTAAACCTTCTGGGGGATTTAGCCAAAAAAACCGTTCGGCCATGCTTCCGTTAagcg                                                                                >  1:471153/1‑147 (MQ=255)
                                                                         aataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGCGCTTATTTAACATTTAGAATTTATCAATCCTTTTTGCAGAATAAGACAAAAAAACCGTTCGGCAATGCTTCCGTTAAGCGGag                                                                             >  2:518236/1‑148 (MQ=255)
                                                                          atatTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGCGGCTATTTAACAATTTGAATTTATCAATACTTCTTGCAGATTTAGCCAAAATACCCCGTCGGGAAAGCTTCCGTTACGGGGAGAgg                                                                          >  1:216901/1‑150 (MQ=255)
                                                                              tttAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGTTGCGCTTTTTGAACATTTTGAATTTATAAAACCTTCTGGCAGATTAAACCAAAAAAACCGGTTGGCAAAGCTTTCGTTaaacggaaagggaagg                                                                     >  2:805692/1‑136 (MQ=255)
                                                                               ttAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAAAAATTAGAATTGATCAAACCTTCTGGCAGATTAAGCCAAAAAAACCGGTCGGCAATGCTTTAGTTAAGCGAAGAGGAATgg                                                                     >  1:150034/1‑150 (MQ=255)
                                                                               ttAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGATTTTGCTGGGCCTATTGAACATTTAAAAATTATCAAAACTTTTTGCAGATTTAGCCAAAATACCCCCTCGGCAATGCTTCAGTTACGCGAagag                                                                           >  1:1112787/1‑144 (MQ=255)
                                                                                tAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTGAACATTTAGAATTTATCAAACCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATgg                                                                     >  2:216030/1‑149 (MQ=255)
                                                                                tAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCCGGGCTTCTTTAAAAATTAGAATTGATCAATCCTTCTGGCGGATTAAGCCAAAAAAACCGTTTGGCAATGCTTTAGTTACGCGAAGAGGAATgg                                                                     >  2:828853/1‑149 (MQ=255)
                                                                                 aGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGCGCTTATTTAACATTTAGAATTGATTAATCCTTCTTGCAGATTAAAACAAAATAACCGTTCGGCAAATCTTTAGTTACGCGaaaggggaagg                                                                     >  2:816386/1‑139 (MQ=255)
                                                                                 aGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCt                                                                  <  1:325617/151‑1 (MQ=255)
                                                                                  gACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTg                                                                 <  2:824688/151‑1 (MQ=255)
                                                                                  gACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTg                                                                 <  1:1012112/151‑1 (MQ=255)
                                                                                    cAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGGGCGCTTATTTAACATTTAGAATTGATCAATCCTTCTGGGAGATTAAGGCAAAAAAACCCTTCGGGAATGCTTCAGTTTCGCGGAGAGGAATGGTCt                                                                  >  1:424520/1‑148 (MQ=255)
                                                                                    cAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGaa                                                               <  2:845266/151‑1 (MQ=255)
                                                                                               cAGGCTTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCttttttt                                                     <  2:1154115/150‑1 (MQ=255)
                                                                                                      tAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTTTTTAACATTTAGAATTGATCAAACCTTCTGGCAGATTAAGCCAAAAAAACCGTTCGGCAATGCTTCAGTTACGCGAAGAGGGATGGTCTGAATTCTTTTTTttgttggt                                              >  2:974530/1‑148 (MQ=255)
                                                                                                             ggggCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTTAAAAATTTGAAATGATTAAACCTTTCGGCAGATTTAGCCCAAAAACCCCGTCGGCCATGCTTCCGGTTCGCCGAGAGGGAAGGGTTTAATTCTTTTTTTGGTTGGTaaaattt                                       >  2:233488/1‑146 (MQ=255)
                                                                                                             ggggCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATtct                                       <  2:51410/150‑1 (MQ=255)
                                                                                                                   gATAGTGAGTCGGTATTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAgt                                <  2:544599/151‑1 (MQ=255)
                                                                                                                         gAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTgg                          <  1:719131/151‑1 (MQ=255)
                                                                                                                         gAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTgg                          <  2:362896/151‑1 (MQ=255)
                                                                                                                            gCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGAtt                       <  1:623453/151‑1 (MQ=255)
                                                                                                                             cGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTa                      <  2:277989/151‑1 (MQ=255)
                                                                                                                                gTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTAtg                    <  1:566752/151‑1 (MQ=255)
                                                                                                                                   ttGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTAtgtg                  <  2:840530/150‑1 (MQ=255)
                                                                                                                                    tGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGat                >  1:1091076/1‑151 (MQ=255)
                                                                                                                                       cTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGattat             >  1:210701/1‑151 (MQ=255)
                                                                                                                                         gcgcTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGattatt            >  2:829211/1‑150 (MQ=255)
                                                                                                                                                  gAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTc  <  2:239202/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
TCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGATTATTAATTTAACTC  >  NC_000913/1245102‑1245397

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