Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC40.3% 101.6 / 31.9 67T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (38/2);  new base C (0/27);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 9.89e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376356‑1376655
                                                                                                                                                           |                                                                                                                                                      
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                                       >  2:776968/1‑151 (MQ=255)
gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCATGCGCGTACTGGAACAGGTCGCCGCGCGTACTGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                                       >  1:1090681/1‑151 (MQ=255)
gCTAACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAaca                                                                                                                                                        <  2:1006778/150‑1 (MQ=255)
 cTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAcc                                                                                                                                                      >  1:151532/1‑151 (MQ=255)
    aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                                       <  2:495256/147‑1 (MQ=255)
     ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAcc                                                                                                                                                      <  1:1031257/147‑1 (MQ=255)
      cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGc                                                                                                                                                 >  1:1153905/1‑151 (MQ=255)
      cgacgATGACTTCGCCGCGTAGCATGGAAGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGCTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCg                                                                                                                                                  >  1:548233/4‑150 (MQ=255)
       gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                                >  1:65894/1‑151 (MQ=255)
         cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                               >  2:61593/1‑150 (MQ=255)
          gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                             <  2:290625/151‑1 (MQ=255)
                       cgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                                 >  1:24913/1‑150 (MQ=255)
                            cgggggcGGCGCCCGAAAAATGGTGAGTGTTTCCCTGCCCGTACGGGAACAGGTCGCCGCGCTTCCGGGACCGGTGTTGTCTCCCGACCCTTAAACCCTCTAAAAGGATAATATGACAACCCCCCCCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                           <  1:299177/145‑1 (MQ=255)
                              tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                          >  1:1063888/1‑150 (MQ=255)
                               ggACGGGAACCGCAAAAGGGGGAGTTATTCCCTGCGGGTAGTGAAACGGGCCGCCCCGCGTCCCGAACCGGTGGTTTATCTCAACCGTTAAACCCGTTATGAGAATAAATTGATAACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGaa                                                                                                                        <  1:526956/151‑1 (MQ=255)
                                  cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                     >  2:823480/1‑151 (MQ=255)
                                   ggCGACCGCAAAATGGGGAGTTATTCCCTGGGCGTACGGGACCGGGCCCCCGCGGGTCCCGGACCGGTGGTGTATCTCGACCCGTAAACCCGCTATCAGGATCATATGACAAACCCCCCCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                           <  1:371312/144‑1 (MQ=255)
                                   gccagccGAAAATTGGTGAGTTATTCCCGGCGCTTCCAGAAACGGGGCGCCGCGCGGCCCGAATCCTTGTTTTATCTCGACCCTTAAACCCACTATCAGGATAATATAACCAACCCCCTCGCCAATGCCCGCCGTTACATCGTCGaaaa                                                                                                                      <  2:1216423/144‑1 (MQ=255)
                                    gCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                    >  1:83896/1‑150 (MQ=255)
                                     cGACCGAAAAAGGGGGAGTTATCCCCCGCGCGTACTGAAAAAGTTCCCCGCGCGCACCGGACCCGGGGTGTTTCTCGACCCGTTAACCCGCTATCAGGATCATATAATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                                   <  2:433555/150‑1 (MQ=255)
                                        cgaaaaaaTGGTGAGTTTTCCCCTCCGTTTAGTGAACCGGGTCGCCGCGCGTACCGGACCAGTGTTGTATCTCAACCCTTAAAACTCTTACCGGCAAAATAGGATCACCCCCCCCGCGGATCGCCGCCGTTACATCGTCGAAAACGa                                                                                                                   <  2:65899/143‑1 (MQ=255)
                                        ccGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                               >  1:605015/1‑151 (MQ=255)
                                        ccGAAAAATGGTGAGTGATCCCCTGCGCTTCCTGGAACAGGCGCCCGCGCGTCCCGAAACCGTGGTGTATCTCGACCCTTTAACCCGCTATCAGGATCATATGATCACCCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                               <  2:90377/151‑1 (MQ=255)
                                        ccGAAAAAGGGGGAGTGATCCCCGGCGCTTACGGAAAAAGGTCGCCGCGCGTCCGGGACCCGGGGTTTCTCCCGACCGTTAAACCCGCTATAAGGATCATATGAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                                  <  2:122678/148‑1 (MQ=255)
                                          gCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                              >  1:1058964/1‑150 (MQ=255)
                                          gCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                              >  1:796330/1‑150 (MQ=255)
                                               aTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                        >  1:656792/1‑151 (MQ=255)
                                                gggggagattctCCCTGCGTGTCTGGGAACGGTCGCCCGCGCTCCCGGACACGTGGTTTTTCCCGGACCCTTAAAACCTCTAGAGGAATAATATAAAAAACCCCCTGGCCGAGACCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                          <  1:23311/137‑1 (MQ=255)
                                                 ggtgggtgTTCCCCGGCGCGTACGGGAACGGGCCGCCGCGCTTCCGGGACCGGTGTTGTACCCCGACCCGTAAACCCGCTACCAGGTTAATATTAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                       <  1:838978/150‑1 (MQ=255)
                                                      gtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     >  2:636055/1‑147 (MQ=255)
                                                         ttccccTGCGCTTCTGGAAAGGGCCCCCGGCGCGTCCGGGACCGGTGTTTTTTCCCAAACCGTAAACCCTCTATAGGGAATAAATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                                 <  1:816275/145‑1 (MQ=255)
                                                           tCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                            >  2:1157013/1‑151 (MQ=255)
                                                                 cgggtatgggAACAGGGCGCCGCGTGTCCCGGAACCGGGGTGTATCTCGAACTTTTAACCCGATATCAGGAACAAAGGATCAGCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                          <  2:197026/139‑1 (MQ=255)
                                                                   cgTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                    >  2:32300/1‑151 (MQ=255)
                                                                     tACGGGAACGGGCCGCCGCGCTTACCGGAACGGTGGTGTATCTCGACCCTTAAAACCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                   <  2:1035057/150‑1 (MQ=255)
                                                                         gaaaCGGGCCGCCGCGCGCCCCGGACCCGTGGTGTACCGCAACCCGTTAACCCGCTAGAAGGATCAAAGGAACACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                                  <  2:732614/145‑1 (MQ=255)
                                                                           aaCGGGCCGCCGCGCGTCCCGGAACCGTGGTGTATCTCGACCCTTAAACCCGCTATCAGAATAATATGAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                            <  1:897183/151‑1 (MQ=255)
                                                                               ggccgccgcGCGTACGGAACCGGGGTTTTACCCCGACCCGTAAACCCGCTTTAAGAATAATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                          <  1:332441/146‑1 (MQ=255)
                                                                                     cgcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                   >  2:747696/1‑150 (MQ=255)
                                                                                            gtcccGGACCGGTGGTTTATCCCAACCCGTTAACCCGCTATCAGAATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                                  <  2:832293/143‑1 (MQ=255)
                                                                                               ccGGACCCGGGTTGTATCTCGAACCGTTAACCCGCTACAAGGAACATATAATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  1:556052/149‑1 (MQ=255)
                                                                                                    aCCGTGGGGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATTAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGACGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                     >  2:182531/1‑151 (MQ=255)
                                                                                                     ccGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                    >  1:782336/1‑151 (MQ=255)
                                                                                                       gtggttTATCTCAAACCGTAAAACCTCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                   <  2:1032090/150‑1 (MQ=255)
                                                                                                       gtggtgTACCTCGACCCTTAAACCCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  2:832152/141‑1 (MQ=255)
                                                                                                         ggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                  >  2:629793/1‑149 (MQ=255)
                                                                                                         ggtgTATCCCGACCCGTAAAACCGTTATCAGGATCATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                                <  2:548602/151‑1 (MQ=255)
                                                                                                            tttATCTCGAACCATAAACCCGCTATAAGAATAAAATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                               <  1:840895/148‑1 (MQ=255)
                                                                                                                tctcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                          >  2:308157/1‑151 (MQ=255)
                                                                                                                 ctcGACCCTTAAACCCGCTATCAGAATCATATGACCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                           <  2:135359/149‑1 (MQ=255)
                                                                                                                 ctcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                                         >  1:641747/1‑151 (MQ=255)
                                                                                                                    gAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAGGATAACCTGGCGCAGGCGCTGCATg                                      >  2:779667/1‑151 (MQ=255)
                                                                                                                     aaCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                      >  2:770529/1‑150 (MQ=255)
                                                                                                                           tAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                               >  1:115787/1‑151 (MQ=255)
                                                                                                                              aaCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                              >  1:1121956/1‑150 (MQ=255)
                                                                                                                               aaccctttataAGGATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACaa                                  <  1:1003933/140‑1 (MQ=255)
                                                                                                                                 ccGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                           >  2:5030/1‑151 (MQ=255)
                                                                                                                                    gttaccAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                          <  1:506568/145‑1 (MQ=255)
                                                                                                                                    cTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                        >  2:671020/1‑151 (MQ=255)
                                                                                                                                        tCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                     >  2:825658/1‑151 (MQ=255)
                                                                                                                                         cAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGAAGGCGCTGCATGACAACCGCGACCTGCTCg                    >  2:663037/1‑151 (MQ=255)
                                                                                                                                          aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg                   >  1:819004/1‑151 (MQ=255)
                                                                                                                                            gatattattaacAACCCCCCCCCCGACCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                    <  1:508612/137‑1 (MQ=255)
                                                                                                                                            gATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATCACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                 >  1:393270/1‑151 (MQ=255)
                                                                                                                                                 tatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg            >  1:259400/1‑151 (MQ=255)
                                                                                                                                                  atGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc           >  1:1081852/1‑151 (MQ=255)
                                                                                                                                                   tGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                 <  2:966706/144‑1 (MQ=255)
                                                                                                                                                      tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat       >  1:516238/1‑151 (MQ=255)
                                                                                                                                                        aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg     >  1:467471/1‑151 (MQ=255)
                                                                                                                                                        aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg     >  2:470363/1‑151 (MQ=255)
                                                                                                                                                        aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg     >  2:400032/1‑151 (MQ=255)
                                                                                                                                                           acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  1:792602/1‑151 (MQ=255)
                                                                                                                                                           acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  1:1072186/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 < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: