Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA46.0% 72.4 / 45.9 74intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (35/5);  new base A (0/34);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 4.36e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  NC_000913/5011‑5306
                                                                                                                                                    |                                                                                                                                                    
tcatcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgaga                                                                                                                                                     >  2:220753/1‑150 (MQ=255)
tcatcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgaga                                                                                                                                                     >  1:379539/1‑150 (MQ=255)
 catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGTAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAATTTAGAGATTAGGATTGGGGAGaa                                                                                                                                                    >  1:515056/1‑150 (MQ=255)
  atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                   >  1:310378/1‑150 (MQ=255)
    caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaa                                                                                                                                                 <  2:553698/150‑1 (MQ=255)
      gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaa                                                                                                                                              <  2:486818/151‑1 (MQ=255)
        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑cc                                                                                                                                           >  2:812015/1‑151 (MQ=255)
        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑cc                                                                                                                                           <  2:350713/151‑1 (MQ=255)
            tCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑ccgccg                                                                                                                                       >  2:1078063/1‑151 (MQ=255)
               aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTtc                                                                                                                                    >  1:952161/1‑151 (MQ=255)
                  cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTtctc                                                                                                                                  >  2:471257/1‑150 (MQ=255)
                    ttATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCa                                                                                                                                 >  1:1135888/1‑149 (MQ=255)
                     tATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCg                                                                                                                              >  1:1026583/1‑151 (MQ=255)
                           aaTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAAt                                                                                                                        >  2:981405/1‑151 (MQ=255)
                              cAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAAtctc                                                                                                                     >  2:471405/1‑151 (MQ=255)
                                ggCCGGGTTTTTTTTTATGCAGCCCGGTTTTTTTATGAAAAAATTAGGAAAAAAAATGACAGGAAAAAAGAAAAAATCCTAAATAAAGGCGTTAATTAAGAGATTGGGATTGCGAAAAATAAAAC‑CCGCCGTTCTCATCGAGTAAtctc                                                                                                                     <  2:675695/149‑1 (MQ=255)
                                ggCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                   >  2:35481/1‑151 (MQ=255)
                                    gggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                >  1:493381/1‑150 (MQ=255)
                                       tttGTTTTTAGGCAGCCCGTTTTTTTAATAAAAAAATAAGGAAAAAAAAGGAGAGGAAAAAAGGAAAAATTCAAAAAAAAGGCGTTAACTAAGAGATAAGGAGGGGGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                              <  1:255009/149‑1 (MQ=255)
                                         tGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                          >  2:161945/1‑151 (MQ=255)
                                           tttttAGGCAGCCCGTTTTTTTTAAAAAAAAATTATGGAAAAAAATGAAAGGAAAAAGGGAAAAATTCAAAAAAAAGGCGGTACTTTAGAGATTGGTATGGGGAAAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                         <  2:448965/149‑1 (MQ=255)
                                             tttATGCCGCCGGGTTTTTTTATGAAAAAATTATGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTCAATAAAGGCGTTAACTTAGAGTTGAGGATGGCGGAAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                         <  1:1049834/148‑1 (MQ=255)
                                              ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                      >  2:2402/1‑150 (MQ=255)
                                              ttAGGCAGCCCGGTTTTTTTATGAAAAAATTATGAAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTAATTTAGAGATTAGGATGGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                      <  2:624561/150‑1 (MQ=255)
                                                 tGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                  >  1:271156/1‑151 (MQ=255)
                                                 tGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                  >  1:719352/1‑151 (MQ=255)
                                                   caccccGGTTTTTTTATAAAAAATTAATGAAAAAAAAGAAAAGGAAAAAAGAAGAAATTCTCAATAAGGGCGGTAATTAAGAGTTTAGGATGGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                    <  2:331406/144‑1 (MQ=255)
                                                     gCCCGGC‑TTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                              >  2:565672/1‑150 (MQ=255)
                                                      ccgggttttttaATAAAAAAATTAGGAAAAAAAATGAGAGGAAAAAAGAAAAAATTCTAAAAAAATGCGTAAACTTAGATATTAGGTTGGGGAAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                     <  1:78804/137‑1 (MQ=255)
                                                      cccGGTTTTTTTAAAAAAAAATTAGGGAAAAAAATGAAGGGAAAAAAGGAAAAATTCTAAAAAAATGCGTTAACTTAGAAATTAGGATGGGGGAAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATGACgg                                                                                                 <  2:934361/147‑1 (MQ=255)
                                                       ccGGTTTTTTTATAAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGTTAACTTAGAGATTAGGATTGGGGAAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAAt                                                                                            <  2:593097/151‑1 (MQ=255)
                                                         ggtttttttAAAAAAAATTTATGGAAAAAAATGAAAGGAAAAAAGGAAAATTTCAAAATAAAGGCGGTAATTTAGAGATTAGGTTTGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGa                                                                                          <  2:587501/148‑1 (MQ=255)
                                                           tttttttAAAAAAAAATTATGGAAAAAAAGGACAGGAAAAAAGAAAAATTTCCCAAAAAAGGCGTTAATTTAGAGATTAGGTGGGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                         <  1:42610/149‑1 (MQ=255)
                                                            ttttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                        >  2:658538/1‑150 (MQ=255)
                                                            ttttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                        >  2:6240/1‑150 (MQ=255)
                                                            ttttttATGAAAAAATTAGGAAAAAAAAAGAGAGGAAAAAAGGAAAAATTCCAAAAAAATGTGGTAACTAAGAGATTAGGTTGGCGGAAAAAACCAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                        <  1:115125/150‑1 (MQ=255)
                                                              ttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                      >  2:81418/1‑150 (MQ=255)
                                                              ttttAAAAAAAAATTAGGAAAAAAAAGGACAGGAAAAAAGGAAAAATTCTAAATAAAGCCGTTACTTTAGATTTTAGGATGGGGGAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                      <  2:1037930/150‑1 (MQ=255)
                                                                 tATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                   >  2:25212/1‑150 (MQ=255)
                                                                 tATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                   >  2:1141658/1‑150 (MQ=255)
                                                                  aaaaaaaaaTAAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAATAAATGGGTTAACTTAGAAATTGGGTTGGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                   <  1:714019/143‑1 (MQ=255)
                                                                     aaaaaaTTATGGAAAAAAATGACAGGAAAAAAGGAGAAATTCTCAATAAATGAGGTATTTTAGAGATTAGGATTGGGGAAAATAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                               <  2:1166375/147‑1 (MQ=255)
                                                                              tGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaa                                                                      >  2:686250/1‑150 (MQ=255)
                                                                               gaaaaaaaaTGACAGGAAAAAAGGAGAAATTCCAAATAAATGCGGTAAATGAGAGATTAGGAGGGGGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                        <  2:740422/143‑1 (MQ=255)
                                                                                 aaaaaaaaGAAAGGAAAAAAGGAAATATTCAAAAAAAGTGCGTTAATTTAGAAAGAAGTATGGGGAAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaa                                                                      <  1:504682/145‑1 (MQ=255)
                                                                                 aaaaaaaTGACAGGAAAAAAGGAAAAATTCTAAAAAAATGGGGTAATTTAGAAATTAGGATTGGGGAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAGTCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaaa                                                                   <  2:1083000/148‑1 (MQ=255)
                                                                                    aaaagaacgggaaaaaaGGAAAAATTCACAAAAAATGCGGTATTTAAGATATGAGGATTGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                               <  2:221405/139‑1 (MQ=255)
                                                                                    aaaaTGAAAGAAAAAAAGGAAAAATTCTAAAAAAAGGGGTAAATTTAGAGATAAGGATGGGGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                               <  1:365290/151‑1 (MQ=255)
                                                                                         gACAGGGAAAAAGGAGAAATTATCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTACGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                           >  1:74102/1‑150 (MQ=255)
                                                                                             ggaaaaaaGGAAAAATTCTAAAAAAAGGCGTAAACTAAGAGATTAGGATGGGGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                       <  2:603828/147‑1 (MQ=255)
                                                                                             ggaaaaaaGAAAAAATTCTAAATAAATGCGTTAATTTAGAGATTAGGTTTGGGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                        <  2:650778/146‑1 (MQ=255)
                                                                                               aaaaaaGAAGAATTTCTAAAAAAAGGGGGTAACTAAGAGATGAGGAGTGGGAAAAAAAATAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                         <  2:1043327/145‑1 (MQ=255)
                                                                                                aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                    >  2:246576/1‑150 (MQ=255)
                                                                                                 aaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                  >  2:697214/1‑151 (MQ=255)
                                                                                                  aaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                 >  2:90649/1‑151 (MQ=255)
                                                                                                  aaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                 >  1:108349/1‑151 (MQ=255)
                                                                                                   aagaaaaatttCTAAAAAAAGACGGTAATTTAAAATTTAGGA‑TGGTGAAAAAAACAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                                <  2:182869/142‑1 (MQ=255)
                                                                                                         aaaTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                          >  2:990149/1‑151 (MQ=255)
                                                                                                          atttCTAAAAAAATGGGATAACTTAGAGATTAGGTTTGGGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                         <  1:1098034/149‑1 (MQ=255)
                                                                                                                  aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                 >  2:1071145/1‑151 (MQ=255)
                                                                                                                  aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                 >  2:742120/1‑151 (MQ=255)
                                                                                                                        gCGGTAACTAAGAGTTTAGGATGGGGGAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                           <  2:203231/151‑1 (MQ=255)
                                                                                                                         cGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                          >  1:172876/1‑151 (MQ=255)
                                                                                                                             accTTAGAGATTAGGATTGCGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                          <  2:405960/145‑1 (MQ=255)
                                                                                                                             aaCTTAGAGATTAGGATTGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCt                        <  2:661554/149‑1 (MQ=255)
                                                                                                                             aaCTTAGAGATTAGGATGGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                      <  2:526629/151‑1 (MQ=255)
                                                                                                                                ttAGAGATTAGGATTGCGGAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCa                    <  2:1006153/150‑1 (MQ=255)
                                                                                                                                    agaTTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                >  1:654607/1‑150 (MQ=255)
                                                                                                                                       ttgggATGGCGGAAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc             <  1:471240/147‑1 (MQ=255)
                                                                                                                                       ttAGGTTGGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgca               <  2:831457/148‑1 (MQ=255)
                                                                                                                                          ggATTGCGGAGAATAACAA‑CCGCGGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac           >  1:904692/1‑149 (MQ=255)
                                                                                                                                          ggATGGGGGAAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcaca          <  2:471924/150‑1 (MQ=255)
                                                                                                                                              tGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag              <  1:1078026/142‑1 (MQ=255)
                                                                                                                                              ggCGGAGAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCt      <  2:1161605/149‑1 (MQ=255)
                                                                                                                                               ggggAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTg     <  2:599656/148‑1 (MQ=255)
                                                                                                                                                 ggAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg  >  1:8315/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  NC_000913/5011‑5306

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