Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC41.8% 184.1 / 21.2 98T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (46/11);  new base C (0/41);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 8.76e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376356‑1376655
                                                                                                                                                        |                                                                                                                                                      
gctcgccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAaca                                                                                                                                                        <  2:456620/145‑1 (MQ=255)
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                                       >  4:1583995/1‑151 (MQ=255)
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                                       >  4:106632/1‑151 (MQ=255)
   tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                                    >  2:512615/1‑151 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                                    <  4:1384861/150‑1 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                                   <  3:166333/151‑1 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAcc                                                                                                                                                      <  2:1403119/147‑1 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTca                                                                                                                                                  <  2:1281425/151‑2 (MQ=255)
       gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                                <  2:1354119/151‑1 (MQ=255)
         cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGa                                                                                                                                              <  4:993069/151‑1 (MQ=255)
           aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                             <  2:1672766/150‑1 (MQ=255)
            tGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGc                                                                                                                                           >  4:1673665/1‑151 (MQ=255)
            acgctcttccgatctgCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGc                                                                                                                                           >  2:1467651/16‑151 (MQ=255)
             gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGc                                                                                                                                           <  4:401992/150‑1 (MQ=255)
                   gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGt                                                                                                                                    >  4:1536185/1‑151 (MQ=255)
                     cgcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                                   >  4:1000019/1‑150 (MQ=255)
                            ccTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                           >  1:12725/1‑151 (MQ=255)
                             cTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                          >  3:212784/1‑151 (MQ=255)
                               ggAGGGCGCCCGAAAAAGGGGGAGTTTTCCCCTGCGCTTACGGGAACGGGTCCCCGCGCGCCCCGGACCGGTGGTGTATCGCGACCCTTTAACCCGTTATCAGGATCATATGACAACCCCCCTCGCCGATGCCCGCCGTTACAtcgt                                                                                                                            <  2:273206/147‑1 (MQ=255)
                               ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                         >  4:766618/1‑150 (MQ=255)
                                   ggCGCCCGAAAAAGGGTGAGTGTTCCCCGGCGCGTACGGAAACAGGCCGCCGCGCGTCCCGGACCGGGGGTGTATCTCAACCCGTAAACCCGCTATCAGGATAATAGGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                    <  2:998722/151‑1 (MQ=255)
                                   ggCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                    >  4:1314221/1‑151 (MQ=255)
                                   ggCGACCGAAAAATGGTGAGTGATCCCCTGCGCTTACGGGAACGGGTCGCCGCGCGTACCGGACCCGTGGTTTTTCTCGAACCTTTAACCCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                    <  3:352022/151‑1 (MQ=255)
                                     cGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                                  >  2:117531/1‑151 (MQ=255)
                                     cGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                                  >  4:389940/1‑151 (MQ=255)
                                      gcccGAAAAATGGGGGGTGATTCCCGGCGCGTACGGGAACAGGTCGCCGCGCGTCCCGGAACGGTGTTTTTTCTCGAACGTTTAAACCGCTATAAGGATCATATGATCACCCCCCGCGCCGAGGCCCGCCGTTACATCGTCGAAAACGa                                                                                                                   <  4:28941/147‑1 (MQ=255)
                                       aCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                                >  1:379926/1‑151 (MQ=255)
                                        ccGAAAAATGGTGGGGGTTCCCGTGCGCTTACGGGACGGGGTCCCCGCGCGTCCCGAACCGGTGGTGTATCGCGAACGTTAAACCCGATACCAGGATCATATGACAAACCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                                 <  4:1243717/149‑1 (MQ=255)
                                            aaaaTGGTGAGTTTTCCCCGGCGCGTACTGAAACAGGCCCCCGCGCGTCCCGGACCGGGGGTTTATCCCAACCCTTAAACCCGCTACCAGGATCATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                             <  2:632749/149‑1 (MQ=255)
                                            aaaaTGGTGAGTTATTCCCTGCGCGTCCGGAAACAGTCCGCCGCGCCTACCGGACCCGGGGTTTATCTCGACCCGTTACCCCGCTATCAGGATCAAATAACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                           <  2:912597/151‑1 (MQ=255)
                                            aaaaTGGTGAGTGATTCCCTGCGCGTACGGGAACAGGTCGCCGCGCGTCCCGGACCCGTGGTTTATCTCGACCCTTTAAACCGCTATAAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                            <  4:241469/150‑1 (MQ=255)
                                            aaaaGGGGGAGTGTTCCCCTGCGCTTACGGGACCAGGTCGCCCCGCTTCCCGGAACCGTGTTGTATCTCGAACATTAAAACCCCTATAAGGATAAAATAATCACCCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                           <  3:92970/151‑1 (MQ=255)
                                                 ggggagtttttCCCTGGGCGTACGGGAACAGGTCGCCGCGCGTCCCGAACCGGTGGTGTATCGCGACCCGTTAACCCGTTACCAGGATAATAAGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                         <  4:1402464/139‑1 (MQ=255)
                                                  gtgagtgaTTCCCGGCGCGTACTGAAACGGGTCCCCGCGCGTCCCGGACCGGTGGTTTATCTCGACCCGTTAACCCGCTACAAGAATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                        <  2:453757/148‑1 (MQ=255)
                                                  gggagtgtcccccGGCGCGTACTGAAACGGGTCGCCGCGCGCACCGGACCCGTGTTGTTTCTCGACCTTTAAACCCGCTATCAGGAAAATTGGATCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     <  3:404371/141‑1 (MQ=255)
                                                    gagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     >  4:1291335/1‑149 (MQ=255)
                                                       tgtttCCCTGCGCTTACTGGAACGGGTCGCCGCGCGTACCGGACCGGTGGTGTATCTCGAACCGTAAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                                 <  4:744524/147‑1 (MQ=255)
                                                       tgttccccGGCGCGTACTGAAAAGGGCCGCCGCCCGTACGGAAACCGTGGTTTATCGCAACCCTTAAACCCTATATCAGAATAAATTGACCACCCCCCACGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACa                                                                                                    <  2:1447374/142‑1 (MQ=255)
                                                       tgaTTCCCTGCGCGTACGGGAACGGGTCGCCGCGCGTACCGGACCGGTGGTTTATCTCAACCCGTTAACCCGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                                   <  4:882519/148‑1 (MQ=255)
                                                         aTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                              >  2:1212503/1‑151 (MQ=255)
                                                          tccccTGCGCGTACGGGAACAGGCCCCGGCGCGTCCGGAACCCGTGGTTTACCTCAAACCGTAAACCCGTTATAAGGATCATATAAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                                 <  3:246668/145‑1 (MQ=255)
                                                            cccGGCGCGTACGGAAACGGGTCGCCGCGCGTCCCGGACCGGTGGTTTATCTCAACCCGTAAACCCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                           <  4:1455157/151‑1 (MQ=255)
                                                             ccTGCGCGTGCTGAAACGGGTCGCCGCGCGTACGGGAACCGTGGTGTATCCCGAACCTTAAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                          <  4:1583673/151‑1 (MQ=255)
                                                                    gTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                   >  4:1235783/1‑151 (MQ=255)
                                                                    gTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                   >  4:822187/1‑151 (MQ=255)
                                                                      aCTGGAACAGGCCGCCGCGCGTCCGGAACCGGGGGTGTATCGCGACCCGTAAACCCGTTATAAGGTTAATATGATCAACCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                     <  1:444141/147‑1 (MQ=255)
                                                                         ggAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                              >  2:1603610/1‑151 (MQ=255)
                                                                               ggTCGCCGCGCGTCCCGGACCCGTGGTGTATCTCGACCCGTAAACCCGTTACAAGGATCATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                        <  4:1640464/151‑1 (MQ=255)
                                                                                   gccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    >  3:422671/1‑151 (MQ=255)
                                                                                   gccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    >  4:902239/1‑151 (MQ=255)
                                                                                    ccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                    >  2:95711/1‑150 (MQ=255)
                                                                                        gcgtccggaaaCGGTGGTGTATCTCAACCCTTAAACCCGCTATCAGGATCATATGATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                <  4:171363/141‑1 (MQ=255)
                                                                                        gcgTGCCGGACCCGGGGTGTATCTCGAACCTTTAAACCGTTAGCAGGATAATATAATCAACACCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                               <  1:110149/151‑1 (MQ=255)
                                                                                         cgcaCCGGAACCGTGGTTTATCTCGACCCGTTAACCCGCTATAAGGATCATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                               <  4:695294/147‑1 (MQ=255)
                                                                                         cgTACCGAACCCGGGGGGTATCCCAACCCTTAAACCCGCTATCAGAATAATATAATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                                 <  2:1389552/148‑1 (MQ=255)
                                                                                          gTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  4:249997/1‑148 (MQ=255)
                                                                                           tACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                                >  3:282798/1‑147 (MQ=255)
                                                                                             ccGGACCCGTGGTGTATCTCGACCCGTTAACCCGCTATCAGGATCATATGATCAACACCCTCACCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                               <  3:465819/146‑1 (MQ=255)
                                                                                             ccGGAACCGTGTTGTATCTCGAACCGTAAACCCGTTATCAGGATAATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  3:269732/149‑1 (MQ=255)
                                                                                               ggAACGGTGGTGTCTCTCGACCCTTAAACCCGCTACCAGGATAATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                           <  2:1434003/148‑1 (MQ=255)
                                                                                               ggAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                        >  2:358661/1‑151 (MQ=255)
                                                                                               cgggaccgggggtgTATCTCGACCGGTAAACCCGCTATCGGGATCATATAATAACCCCCCGCGCTGTTTCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                       <  3:285666/128‑1 (MQ=255)
                                                                                                gAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                       >  3:442939/1‑151 (MQ=255)
                                                                                                  acccGTGGTGTATCCCAAACCTTAAACCCGCTATCAGGATAATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                       <  2:1543627/148‑1 (MQ=255)
                                                                                                    ccGTGGTTTACCCCGACCCGTAAACCCCCAATCAGAATCATAGGACAACCCCCCCCGCGGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                        <  2:1671637/147‑1 (MQ=255)
                                                                                                    ccGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                    >  3:420470/1‑151 (MQ=255)
                                                                                                     cGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                   >  3:297080/1‑151 (MQ=255)
                                                                                                     cGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                   >  4:92921/1‑151 (MQ=255)
                                                                                                        tggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                  >  3:232514/1‑150 (MQ=255)
                                                                                                           tgtATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                              >  4:892369/1‑151 (MQ=255)
                                                                                                           tgtATCCCGAAGCTTAAAACCGTTAGCAGAATAATATGATAACCCCCCCCCCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                               <  1:9949/150‑1 (MQ=255)
                                                                                                            gtATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                             >  2:1324602/1‑151 (MQ=255)
                                                                                                             tATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcgc                                            >  3:2958/1‑151 (MQ=255)
                                                                                                               tctcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                          >  2:4477/1‑151 (MQ=255)
                                                                                                               tctcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                          >  4:594671/1‑151 (MQ=255)
                                                                                                                ctcGAACCGTTAACCCGCTATCAGGATAATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                              <  4:320326/146‑1 (MQ=255)
                                                                                                                  cGACCCGTAAACCCGATAACAGGATAATATAATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCa                                        <  2:1234228/150‑1 (MQ=255)
                                                                                                                   gAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                          >  3:167445/1‑115 (MQ=255)
                                                                                                                    acccGTTAACCCGCTACCAGAATCATATAAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                             <  4:1410291/141‑1 (MQ=255)
                                                                                                                    aaCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                      >  3:334141/1‑150 (MQ=255)
                                                                                                                       ccttAAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                          <  4:622134/142‑1 (MQ=255)
                                                                                                                           aaaCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                              >  4:971735/1‑151 (MQ=255)
                                                                                                                            aaccctttATCAGTATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                                <  4:551952/142‑1 (MQ=255)
                                                                                                                                     tCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                     >  4:1568699/1‑151 (MQ=255)
                                                                                                                                        gattCATATGACAACCCCCCCCGCCGTGCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                     <  1:238366/145‑1 (MQ=255)
                                                                                                                                        gaaTCATATAACCAACCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                     <  4:95873/146‑1 (MQ=255)
                                                                                                                                         gATAATATAATCAACCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg                   <  1:252472/149‑1 (MQ=255)
                                                                                                                                          ataatatGATCTACACCCTCGCCGATGCCCGCCGTTACATCTTCGAGAACGATCTTAACCATGTACAGATTATCGGCGATTGCTATCACATGAATAACGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCCCGGTc                 >  2:300254/4‑150 (MQ=255)
                                                                                                                                             atagGACCAACCCCCTCGCCTATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCa                <  2:1395606/148‑1 (MQ=255)
                                                                                                                                              tatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGTGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg            >  2:1272145/1‑151 (MQ=255)
                                                                                                                                              tatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg            >  2:812740/1‑151 (MQ=255)
                                                                                                                                               atGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc           >  4:640098/1‑151 (MQ=255)
                                                                                                                                               atGAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCa                <  2:1173254/146‑1 (MQ=255)
                                                                                                                                                  aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat         >  1:347177/1‑150 (MQ=255)
                                                                                                                                                   tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat       >  3:137012/1‑151 (MQ=255)
                                                                                                                                                    cAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAtt      >  4:1610541/1‑151 (MQ=255)
                                                                                                                                                     aaCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAtt      <  2:821099/150‑1 (MQ=255)
                                                                                                                                                     aaCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAtt      <  3:342101/150‑1 (MQ=255)
                                                                                                                                                       cacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCg   >  4:1154519/1‑151 (MQ=255)
                                                                                                                                                        acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  3:78239/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                      
GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376356‑1376655

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