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

ACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376360‑1376652
                                                                                                                                                    |                                                                                                                                                   
aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                                >  1:803502/1‑151 (MQ=255)
 ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                                <  1:1103824/150‑1 (MQ=255)
   gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                             <  1:111827/151‑1 (MQ=255)
   gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                             >  1:630463/1‑151 (MQ=255)
       aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                         >  2:94854/1‑151 (MQ=255)
         gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGcc                                                                                                                                       >  1:1048778/1‑151 (MQ=255)
            cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCAATGCccgc                                                                                                                                    >  1:990153/1‑151 (MQ=255)
                       gctggggcGGCGACCGAAAAAGGGTGGATATTTCCCGGTGTTTAAGGAAACAGGCCGCCGCGCGTCCCGGACCCGTGGTTTATCCGAACCATTAAAACCCTTTTAGGAATATTTTGAACAACCCCCCCGCCGATGCCCGCCGTTACAtc                                                                                                                           <  2:750687/142‑1 (MQ=255)
                           ggAGGGCGCCCGAAAAATGGTGAGTTTTCCCCTGCGCTTACGGAAACAGGTCGCCGCGCGTCCCGAAACGGTGGTGTTTCTCGACCCTTAAACCCCCTATCAGGATAATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGaa                                                                                                                     <  1:841897/151‑1 (MQ=255)
                             aCGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                                    >  2:1130973/1‑150 (MQ=255)
                              cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                  >  1:778427/1‑151 (MQ=255)
                                 cGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                               >  2:534256/1‑151 (MQ=255)
                                    ccGAAAATTGGGGGGTGATCCCCGGCGCGTACTGGAACGGGCCCCCGCGCGTCCCGAACCCGTGGTTTATCCCGACCCTTAAACCCGCTATCAGGATCATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                             <  2:247769/150‑1 (MQ=255)
                                    ccGAAAAGTGGTGGGTGATTCCCTGCGCGTACTGGAGCGGTTCGCCGCGCGTACCGAACCGGTGGTGTATCTCGAACCGTTAACCCGCTATCAGGATCATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                               <  1:1212083/148‑1 (MQ=255)
                                    ccGAAAAAGGGTGAGTGATCCCCTGCGCGTACTGGAACGGGTCGCCGCGCGTCCCGGACCCGTGGTTTATCCCGACCCTTTAACCCCCTATCAGGATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                                <  1:1092273/147‑1 (MQ=255)
                                    cagaaaagTGGTGGGTGCTCCCCGGCGCGTCCGGGAACGGGTCGCCGCGCGCCCGGGACCCGTGGTGTACCTCAACCCTTAAACCCGCTATCAGGATAATAAGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                                <  2:258645/143‑1 (MQ=255)
                                       aaaagTGGTGAGTGCTCCCCTGCGCTTCCTGAACCAGGCCCCCGCCCGTACCGAACCCGGGGTTTCCCTCAACCCTTAAACCCGCTATAAGCATCATAGAATCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                         <  1:617373/150‑1 (MQ=255)
                                          aaTGGTGGGTGATCCCCTGCGCGTACTGAAACAGGTCGCCGCGCTTACCGAACCCGTGGTTTATCCCGACCCTTTAACCCGCTATAAGGATCATATAATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                      <  1:118051/151‑1 (MQ=255)
                                             ggtggttgttccccGGCGCGTACTTGAACAGGCCCCCGCGCGTCCGGGACCCGTGTTGTACCCCGACCCGTAAACCCGCTATAAGAATAATATAACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                                   <  2:1065507/140‑1 (MQ=255)
                                                 agtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                  >  2:1010621/1‑148 (MQ=255)
                                                 agtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                  >  1:1159760/1‑148 (MQ=255)
                                                  gttatccccTGGGCGTCCTGAAACGGGTCGCGGCGCGTCCCGAACCCGGGGTTTATCCCGACCCTTTAAACCGCTACAAGAATCATAGGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                              <  1:821789/145‑1 (MQ=255)
                                                      ttCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                          >  2:74233/1‑151 (MQ=255)
                                                      ttCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                          >  2:57106/1‑151 (MQ=255)
                                                        cccTGCGCTTACGGAAACGGGTCGCCGCGCGTACCGGACCCGGGGTTTATCTCGACCCTTTAACCCGCTATCAGGATAATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                        <  1:536143/151‑1 (MQ=255)
                                                              gcgTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                  >  1:542088/1‑151 (MQ=255)
                                                                gTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                >  1:136685/1‑151 (MQ=255)
                                                                   cTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             >  1:947244/1‑151 (MQ=255)
                                                                    tGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             >  2:416624/1‑150 (MQ=255)
                                                                    tGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             >  1:1124728/1‑150 (MQ=255)
                                                                    gggaacggtccgccgcGCGTCCGGGACCGGTGGTGTTCCCCGACCGGTTAACCCGCTATCAGAATCATATGACAACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                <  2:830020/137‑1 (MQ=255)
                                                                     ggAACAGGCCGCCGGGCGTCCGGCACCCGTGGTTTATCCCGACCCGTAAACCCGCTATCAGGTTAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             <  2:934391/149‑1 (MQ=255)
                                                                       aaCAGGCCGCCGCGCGTCCCGAACCGGTGGTGTACCCCGACCCGTAAACCCGTAACCAGGATCATATGAACAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                           <  2:426984/149‑1 (MQ=255)
                                                                              cgccgcGCGTACCGGAACGGTGGTGTATCTCGAACCGTTAACCCGCTATCAGAATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                                  <  1:120246/151‑1 (MQ=255)
                                                                                     cgTCCCGGACCCGTGGTGTCTCTCGACCCTTAAACCCCTTATAAGAATAATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                 <  2:432603/145‑1 (MQ=255)
                                                                                     cgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  2:907413/1‑149 (MQ=255)
                                                                                     cgTACCGAACCCGTGGTGTTTCTCGAACCTTAAAACCGCTATCAGGATAATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                              <  1:804298/148‑1 (MQ=255)
                                                                                        cccGGAACCGTGTTGTATCTCAACCCTTTAACCCGTTACCAGGATCATATATCCAACCCCCTCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                            <  1:1145615/146‑1 (MQ=255)
                                                                                           gaacccgggttgtATCTCCACCCTTTAACCCCCTATAAGGATCATAGAATAAACCCCCTCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                           <  1:534541/135‑1 (MQ=255)
                                                                                             acccGGGGTGTATCTCGACCCTTAAAACCCCTATCGGGATCAAATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                   <  2:998090/149‑1 (MQ=255)
                                                                                             aaCCGTGGTGTACCCCGACCCGTAAAACCGCTATCAGGATCATATGATAAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                   <  1:40956/151‑1 (MQ=255)
                                                                                              aCCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                  >  1:341465/1‑151 (MQ=255)
                                                                                                cGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                >  2:101831/1‑151 (MQ=255)
                                                                                                 gtggttTATCTCGAACCGTAAACCCCCTATAAGAATAATATGATCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                   <  2:1022698/147‑1 (MQ=255)
                                                                                                 gtggttTACCCCGAACCTTAAAACCGCTATCAGTATAATATGATCAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                 <  2:902863/149‑1 (MQ=255)
                                                                                                 gtggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                               >  1:1180087/1‑151 (MQ=255)
                                                                                                 gtggtgATTCCCTAACCGTTAAACCCCAATCAGAACAATATAAACACCTCCCTCGCGCCTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                               <  2:705183/151‑1 (MQ=255)
                                                                                                   ggtgTATCCCGAACCTTTAACCCGCTATAAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCa                                              <  1:205398/150‑1 (MQ=255)
                                                                                                     tgtATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                           >  1:709763/1‑151 (MQ=255)
                                                                                                         tcccaaaCCTTAAAACCGCTATCAGGATCATATGAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                           <  2:407296/142‑1 (MQ=255)
                                                                                                          ctcGACCCTTAAACCCGTTATCGGAATCATATAATCAACCCCCCCCCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                          <  2:268959/147‑1 (MQ=255)
                                                                                                          cccaacccGTTAACCCGCTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                        <  1:349632/143‑1 (MQ=255)
                                                                                                              aaCCTTAAAACCGCTATGGGGATCATATGATCAACACCCCCCCGGACCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCa                                     <  1:597296/148‑1 (MQ=255)
                                                                                                               aCCCTTAAACCCGCTATCAGGATAATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                    <  2:872440/149‑1 (MQ=255)
                                                                                                                  gTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                              >  1:877833/1‑151 (MQ=255)
                                                                                                                     aaaCCGCTAAAAGAATCATATGAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                             <  1:781721/149‑1 (MQ=255)
                                                                                                                      accccccTACCAGGAACATATAACAACCCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                           <  2:419886/145‑1 (MQ=255)
                                                                                                                      aaCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                          >  2:331075/1‑151 (MQ=255)
                                                                                                                       aCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                          >  1:752699/1‑150 (MQ=255)
                                                                                                                          cGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       >  2:475869/1‑151 (MQ=255)
                                                                                                                          cGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       >  1:815492/1‑151 (MQ=255)
                                                                                                                          cGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       >  2:236962/1‑151 (MQ=255)
                                                                                                                                tCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                  >  1:756260/1‑151 (MQ=255)
                                                                                                                                   ggATAATAGAATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt               <  2:1085194/151‑1 (MQ=255)
                                                                                                                                    gaataatatGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                 <  2:527774/144‑1 (MQ=255)
                                                                                                                                         atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt          >  1:927432/1‑151 (MQ=255)
                                                                                                                                         atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt          >  2:937562/1‑151 (MQ=255)
                                                                                                                                           atGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg         >  1:485910/1‑150 (MQ=255)
                                                                                                                                                 aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg  >  2:37113/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376360‑1376652

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: