Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC45.7% 84.8 / 36.8 70T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (34/4);  new base C (0/32);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 6.83e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376359‑1376652
                                                                                                                                                       |                                                                                                                                                   
tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAccc                                                                                                                                                  <  2:1319083/150‑1 (MQ=255)
    gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                             >  2:205919/1‑151 (MQ=255)
    gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                             <  1:658430/151‑1 (MQ=255)
      cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                            >  2:1115824/1‑150 (MQ=255)
         tGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                         <  2:54698/150‑1 (MQ=255)
          gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGcc                                                                                                                                       >  2:1262373/1‑151 (MQ=255)
             cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgc                                                                                                                                    >  2:917898/1‑151 (MQ=255)
             cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgc                                                                                                                                    >  1:1289286/1‑151 (MQ=255)
               cgccggGTTGCCTGTACGGCGACCCCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGACCTGTTGTGTACCGCGACCCGTTAAACCGCTATCAGGATCATATGATCAACCCCCTCGCCTATGCccgcc                                                                                                                                   <  2:862636/150‑1 (MQ=255)
                gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCCCCGCGCCTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                                   <  2:185403/149‑1 (MQ=255)
                  cgcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                                >  2:1324785/1‑150 (MQ=255)
                   gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                              >  2:461932/1‑151 (MQ=255)
                      tAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTAGCCGATGCCCGCCGTTACAtc                                                                                                                           >  2:951640/1‑151 (MQ=255)
                       aGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                                          >  1:207611/1‑151 (MQ=255)
                            ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                      >  2:632743/1‑150 (MQ=255)
                                ggCGCCCGAAAAATGGGGGGTGATCCCCTCGGCGTAGTGGAACGGGCCCCCGCGCTTCCGGGACCCGGGGTATATCTCGAACCTTTAACCCGCTAACAGAATCATAAGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                  <  1:554772/150‑1 (MQ=255)
                                 gCGACCGAAAAGTGGTGAGTGATCCCCTGCGCTTACTGGACCGGGTCGCCGCGCGCCCCGGACCCGGGGTGTATCTCGAACCGTTAACCCGCTATCAGGATCATAGGATCACCCCCCGCGCCGTTGCCCGCCGTTACATCGTCGAAAACg                                                                                                                 <  1:844279/150‑1 (MQ=255)
                                     ccGAAAATTGGTGATTGATTCCCGGCGCTTACGGAAACGGGTCGCGGCGCGTCCCGACACGGTGGTGTACCGCGACCCGTAAACCCGCTACCAGGATAATATGACCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                              <  2:655893/149‑1 (MQ=255)
                                     ccGAAAAAGGGTGAGTGATCCCCTGCGCGTACGGGAACGGGTCGCCGCGCGTACCGGACCCGTGGTGTACCTCGACCCGTTAACCCGCTATCAGGATAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                               <  1:1177619/148‑1 (MQ=255)
                                         aaaaGGGTGGGTTATCCCCTGGGCGTACGGAAACAGGTCGCCGCGCGTCCCGGACCGGTGGTGTTTCTCGAACCGTAAACCCCTTATCAGAATAATAGAACCAACCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                           <  1:848253/148‑1 (MQ=255)
                                         aaaaGGGTGAGTGATCCCCTGCGCGTCCGGGACAGGGCCCCGGCGCGTCCGGAACCGGTGTTGTTTCTCAACCTTTTAACCCGTTATCAGTATAATATGATAACCCCCCTCGCGGAGGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                              <  1:361300/145‑1 (MQ=255)
                                          aaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                        >  2:529914/1‑150 (MQ=255)
                                              ggggggtttttCCCGGCGTGTTGTGGAGCGGGCCGCCGCGCCAACGGGAGCCGTGGTTTATCTGGATCCATAAACCCGTTATAGGGATTATATAACAACCCCCCTCCCCGTCGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                                   <  2:801600/142‑1 (MQ=255)
                                                 gagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                  >  1:436769/1‑149 (MQ=255)
                                                   gtgttccccGGCGCGTACGGGAACAGGTCGCCGCGCGTACCGGACCGGGGGTGTATCTCGAACCGTAAACCCGTTATCAGGATCATAGAACCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                               <  1:662569/144‑1 (MQ=255)
                                                      aTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                           >  2:869087/1‑151 (MQ=255)
                                                       tccccTGCGCGTACTGAAACGGGTCCCCGCGCGTACCGGAACGGGGGTTTTTCTCAACCCGTTAACCCGCTATCAGGATAATATGACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                             <  1:412576/146‑1 (MQ=255)
                                                         cccTGCGCGTCCGGAAACGGTGCCCCGCGCGTCCGGAACCCGGGGTTTTTCTCGACCCTTTAACCCGCTATAAGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACa                                                                                                 <  1:730344/142‑1 (MQ=255)
                                                          ccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                         >  1:237225/1‑149 (MQ=255)
                                                           cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                      >  2:517663/1‑151 (MQ=255)
                                                                cgTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                 >  2:803634/1‑151 (MQ=255)
                                                                   aCTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                              >  1:654530/1‑151 (MQ=255)
                                                                   aCTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                              >  2:809146/1‑151 (MQ=255)
                                                                    cgggAACGGGCCGCCGGGGTACCCGAACCGGTGGTGTTCCCCGACCCGTAAACCCCCTACAAGGATAATATAACCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                    <  2:559111/142‑1 (MQ=255)
                                                                    cTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             >  2:404540/1‑151 (MQ=255)
                                                                       gAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                           >  2:95979/1‑150 (MQ=255)
                                                                                 ccgcGCGTACCGAAACCGTGGTGTTTCTCGACCCGTTAACCCGCTATCAGAATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                    <  1:544938/147‑1 (MQ=255)
                                                                                 ccgcGCGTCCGGAACGCGTGTTTTTCCCCGACCCTTAAACCCGCTATAAGAATAAAATAATCACCCCCCCCGCGGTGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                                   <  2:830211/148‑1 (MQ=255)
                                                                                  cgcgcgTCCGGAACCTGTGGTGTATCTCAACCCTTTAACCCTCTATCAGGATCATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                               <  1:594228/151‑1 (MQ=255)
                                                                                      gcgTCCCGGAACGGTGTTGTTTCCCAACCCTTTAACCCGCTATAAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                            <  1:454609/151‑1 (MQ=255)
                                                                                      gcgTCCCGGAACCGTGGTTTATCTCGAACCGTAAACCCGCTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                               <  1:390389/148‑1 (MQ=255)
                                                                                      gcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  1:265979/1‑150 (MQ=255)
                                                                                        gtcccGGACCGGGGGTGTACCTCGACCTTTAAACCCGCTATCAGAAAAAAAAGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                              <  2:174093/144‑1 (MQ=255)
                                                                                        gccccGGACCGGTGTTGTACCCCGAACCGTAAACCCGCAATAAGAATAAAATAAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                           <  2:409340/147‑1 (MQ=255)
                                                                                        gaaCCGAACCCGTGGTTTATCTCAACCCGTAAACCCGCTACAAGGATAATATAAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                               <  1:1298208/144‑1 (MQ=255)
                                                                                        gTACCGGACCGGGGTTGTATCCCAACCCTTAAACCCGCTATCAGAACAATATGCCCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             <  2:1185154/148‑1 (MQ=255)
                                                                                         caCCGAAACCGTGGTGTTTCCCGACCGTTAAACCCGCTTAGAGAACCATATGATCACCCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                    <  2:465404/139‑1 (MQ=255)
                                                                                           ccGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                        >  1:1053325/1‑150 (MQ=255)
                                                                                            cGGAACCGGTGGTTTATCCCGACCGTTTAACCCGTTATCAGGATCATATAATAAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                         <  1:1063700/149‑1 (MQ=255)
                                                                                                  ccGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                 >  2:313912/1‑151 (MQ=255)
                                                                                                    cgttgttttATCTCAACCCTTTAAACCGTTCACCGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                 <  2:597202/143‑1 (MQ=255)
                                                                                                       tggttTATCTGAACCCTTTAACCCGCTATCAGAATCTTATGATCAACCCCCCCCCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                 <  1:881941/148‑1 (MQ=255)
                                                                                                          tttaccccaacccGTTAACCCCCTATCAGGATAATAAGACCCCCCCCCTCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCa                                              <  2:203050/137‑1 (MQ=255)
                                                                                                           gtATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTACACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                            >  2:956156/1‑149 (MQ=255)
                                                                                                                tcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                                      >  1:1060573/1‑150 (MQ=255)
                                                                                                                      ctttaaagccGCTATAGGGCTCATATAAGCACCCCCCTCGTCGATTCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                 <  1:140191/142‑1 (MQ=255)
                                                                                                                      cGTTAAACCGCTATCAGGATCATCTGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGTGCTGCATGACa                                >  1:278179/1‑150 (MQ=255)
                                                                                                                         taaacccGTTATAAGGATAAAATAAACAACCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                 <  1:1116762/142‑1 (MQ=255)
                                                                                                                          aaaCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGCTTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                             >  2:464809/1‑149 (MQ=255)
                                                                                                                            aCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                          >  1:748340/1‑150 (MQ=255)
                                                                                                                              ccGATATCAGAATCATAAGACACCCCCCCCCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                             <  2:52667/146‑1 (MQ=255)
                                                                                                                                   accAGGATAATATAACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       <  2:207605/145‑1 (MQ=255)
                                                                                                                                      agaaTCATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                  <  1:18338/146‑1 (MQ=255)
                                                                                                                                          tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt            >  2:584889/1‑151 (MQ=255)
                                                                                                                                           cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg           >  1:1033755/1‑151 (MQ=255)
                                                                                                                                              ataaccAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCa             <  2:1063121/141‑1 (MQ=255)
                                                                                                                                                  tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat    >  1:153977/1‑151 (MQ=255)
                                                                                                                                                  tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat    >  2:586467/1‑151 (MQ=255)
                                                                                                                                                  tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCCCGGTCATGTGCatat    >  1:692364/1‑151 (MQ=255)
                                                                                                                                                    aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg  >  1:432741/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376359‑1376652

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