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

TCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  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)
        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                           >  2:812015/1‑151 (MQ=255)
        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                           <  2:350713/151‑1 (MQ=255)
            tCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgccg                                                                                                                                       >  2:1078063/1‑151 (MQ=255)
               aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtc                                                                                                                                    >  1:952161/1‑151 (MQ=255)
                  cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                  >  2:471257/1‑150 (MQ=255)
                    ttATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCa                                                                                                                                 >  1:1135888/1‑149 (MQ=255)
                     tATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                              >  1:1026583/1‑151 (MQ=255)
                           aaTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAt                                                                                                                        >  2:981405/1‑151 (MQ=255)
                              cAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                     >  2:471405/1‑151 (MQ=255)
                                ggCCGGGTTTTTTTTTATGCAGCCCGGTTTTTTTATGAAAAAATTAGGAAAAAAAATGACAGGAAAAAAGAAAAAATCCTAAATAAAGGCGTTAATTAAGAGATTGGGATTGCGAAAAATAAAACCCGCCGTTCTCATCGAGTAAtctc                                                                                                                     <  2:675695/149‑1 (MQ=255)
                                ggCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                   >  2:35481/1‑151 (MQ=255)
                                    gggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                >  1:493381/1‑150 (MQ=255)
                                       tttGTTTTTAGGCAGCCCGTTTTTTTAATAAAAAAATAAGGAAAAAAAAGGAGAGGAAAAAAGGAAAAATTCAAAAAAAAGGCGTTAACTAAGAGATAAGGAGGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                              <  1:255009/149‑1 (MQ=255)
                                         tGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                          >  2:161945/1‑151 (MQ=255)
                                           tttttAGGCAGCCCGTTTTTTTTAAAAAAAAATTATGGAAAAAAATGAAAGGAAAAAGGGAAAAATTCAAAAAAAAGGCGGTACTTTAGAGATTGGTATGGGGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                         <  2:448965/149‑1 (MQ=255)
                                             tttATGCCGCCGGGTTTTTTTATGAAAAAATTATGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTCAATAAAGGCGTTAACTTAGAGTTGAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                         <  1:1049834/148‑1 (MQ=255)
                                              ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                      >  2:2402/1‑150 (MQ=255)
                                              ttAGGCAGCCCGGTTTTTTTATGAAAAAATTATGAAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTAATTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                      <  2:624561/150‑1 (MQ=255)
                                                 tGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                  >  1:271156/1‑151 (MQ=255)
                                                 tGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                  >  1:719352/1‑151 (MQ=255)
                                                   caccccGGTTTTTTTATAAAAAATTAATGAAAAAAAAGAAAAGGAAAAAAGAAGAAATTCTCAATAAGGGCGGTAATTAAGAGTTTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                    <  2:331406/144‑1 (MQ=255)
                                                     gCCCGGCTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                              >  2:565672/1‑150 (MQ=255)
                                                      ccgggttttttaATAAAAAAATTAGGAAAAAAAATGAGAGGAAAAAAGAAAAAATTCTAAAAAAATGCGTAAACTTAGATATTAGGTTGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                     <  1:78804/137‑1 (MQ=255)
                                                      cccGGTTTTTTTAAAAAAAAATTAGGGAAAAAAATGAAGGGAAAAAAGGAAAAATTCTAAAAAAATGCGTTAACTTAGAAATTAGGATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATGACgg                                                                                                 <  2:934361/147‑1 (MQ=255)
                                                       ccGGTTTTTTTATAAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGTTAACTTAGAGATTAGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAAt                                                                                            <  2:593097/151‑1 (MQ=255)
                                                         ggtttttttAAAAAAAATTTATGGAAAAAAATGAAAGGAAAAAAGGAAAATTTCAAAATAAAGGCGGTAATTTAGAGATTAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGa                                                                                          <  2:587501/148‑1 (MQ=255)
                                                           tttttttAAAAAAAAATTATGGAAAAAAAGGACAGGAAAAAAGAAAAATTTCCCAAAAAAGGCGTTAATTTAGAGATTAGGTGGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                         <  1:42610/149‑1 (MQ=255)
                                                            ttttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                        >  2:658538/1‑150 (MQ=255)
                                                            ttttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                        >  2:6240/1‑150 (MQ=255)
                                                            ttttttATGAAAAAATTAGGAAAAAAAAAGAGAGGAAAAAAGGAAAAATTCCAAAAAAATGTGGTAACTAAGAGATTAGGTTGGCGGAAAAAACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                        <  1:115125/150‑1 (MQ=255)
                                                              ttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                      >  2:81418/1‑150 (MQ=255)
                                                              ttttAAAAAAAAATTAGGAAAAAAAAGGACAGGAAAAAAGGAAAAATTCTAAATAAAGCCGTTACTTTAGATTTTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                      <  2:1037930/150‑1 (MQ=255)
                                                                 tATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                   >  2:25212/1‑150 (MQ=255)
                                                                 tATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                   >  2:1141658/1‑150 (MQ=255)
                                                                  aaaaaaaaaTAAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAATAAATGGGTTAACTTAGAAATTGGGTTGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                   <  1:714019/143‑1 (MQ=255)
                                                                     aaaaaaTTATGGAAAAAAATGACAGGAAAAAAGGAGAAATTCTCAATAAATGAGGTATTTTAGAGATTAGGATTGGGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                               <  2:1166375/147‑1 (MQ=255)
                                                                              tGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaa                                                                      >  2:686250/1‑150 (MQ=255)
                                                                               gaaaaaaaaTGACAGGAAAAAAGGAGAAATTCCAAATAAATGCGGTAAATGAGAGATTAGGAGGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                        <  2:740422/143‑1 (MQ=255)
                                                                                 aaaaaaaaGAAAGGAAAAAAGGAAATATTCAAAAAAAGTGCGTTAATTTAGAAAGAAGTATGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaa                                                                      <  1:504682/145‑1 (MQ=255)
                                                                                 aaaaaaaTGACAGGAAAAAAGGAAAAATTCTAAAAAAATGGGGTAATTTAGAAATTAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAGTCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaaa                                                                   <  2:1083000/148‑1 (MQ=255)
                                                                                    aaaagaacgggaaaaaaGGAAAAATTCACAAAAAATGCGGTATTTAAGATATGAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                               <  2:221405/139‑1 (MQ=255)
                                                                                    aaaaTGAAAGAAAAAAAGGAAAAATTCTAAAAAAAGGGGTAAATTTAGAGATAAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                               <  1:365290/151‑1 (MQ=255)
                                                                                         gACAGGGAAAAAGGAGAAATTATCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTACGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                           >  1:74102/1‑150 (MQ=255)
                                                                                             ggaaaaaaGGAAAAATTCTAAAAAAAGGCGTAAACTAAGAGATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                       <  2:603828/147‑1 (MQ=255)
                                                                                             ggaaaaaaGAAAAAATTCTAAATAAATGCGTTAATTTAGAGATTAGGTTTGGGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                        <  2:650778/146‑1 (MQ=255)
                                                                                               aaaaaaGAAGAATTTCTAAAAAAAGGGGGTAACTAAGAGATGAGGAGTGGGAAAAAAAATAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                         <  2:1043327/145‑1 (MQ=255)
                                                                                                aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                    >  2:246576/1‑150 (MQ=255)
                                                                                                 aaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                  >  2:697214/1‑151 (MQ=255)
                                                                                                  aaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                 >  2:90649/1‑151 (MQ=255)
                                                                                                  aaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                 >  1:108349/1‑151 (MQ=255)
                                                                                                   aagaaaaatttCTAAAAAAAGACGGTAATTTAAAATTTAGGATGGTGAAAAAAACAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                                <  2:182869/142‑1 (MQ=255)
                                                                                                         aaaTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                          >  2:990149/1‑151 (MQ=255)
                                                                                                          atttCTAAAAAAATGGGATAACTTAGAGATTAGGTTTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                         <  1:1098034/149‑1 (MQ=255)
                                                                                                                  aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                 >  2:1071145/1‑151 (MQ=255)
                                                                                                                  aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                 >  2:742120/1‑151 (MQ=255)
                                                                                                                        gCGGTAACTAAGAGTTTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                           <  2:203231/151‑1 (MQ=255)
                                                                                                                         cGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                          >  1:172876/1‑151 (MQ=255)
                                                                                                                             accTTAGAGATTAGGATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                          <  2:405960/145‑1 (MQ=255)
                                                                                                                             aaCTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCt                        <  2:661554/149‑1 (MQ=255)
                                                                                                                             aaCTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                      <  2:526629/151‑1 (MQ=255)
                                                                                                                                ttAGAGATTAGGATTGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCa                    <  2:1006153/150‑1 (MQ=255)
                                                                                                                                    agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                >  1:654607/1‑150 (MQ=255)
                                                                                                                                       ttgggATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc             <  1:471240/147‑1 (MQ=255)
                                                                                                                                       ttAGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgca               <  2:831457/148‑1 (MQ=255)
                                                                                                                                          ggATTGCGGAGAATAACAACCGCGGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac           >  1:904692/1‑149 (MQ=255)
                                                                                                                                          ggATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcaca          <  2:471924/150‑1 (MQ=255)
                                                                                                                                              tGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag              <  1:1078026/142‑1 (MQ=255)
                                                                                                                                              ggCGGAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCt      <  2:1161605/149‑1 (MQ=255)
                                                                                                                                               ggggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTg     <  2:599656/148‑1 (MQ=255)
                                                                                                                                                 ggAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg  >  1:8315/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  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: