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

ACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376360‑1376652
                                                                                                                                                    |                                                                                                                                                   
aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTc                                                                                                                                                >  1:803502/1‑151 (MQ=255)
 ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTc                                                                                                                                                <  1:1103824/150‑1 (MQ=255)
   gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGcc                                                                                                                                             <  1:111827/151‑1 (MQ=255)
   gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGcc                                                                                                                                             >  1:630463/1‑151 (MQ=255)
       aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATg                                                                                                                                         >  2:94854/1‑151 (MQ=255)
         gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGcc                                                                                                                                       >  1:1048778/1‑151 (MQ=255)
            cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCAATGCccgc                                                                                                                                    >  1:990153/1‑151 (MQ=255)
                       gctggggcGGCGACCGAAAAAGGGTGGATATTTCCCGGTGTTTAAGGAAACAGGCCGCCGCGCGTCCCGGACCCGTGGTTTATC‑CGAACCATTAAA‑ACCC‑TTTTAGGAATATTTTGAACAACCCCCCCGCCGATGCCCGCCGTTACAtc                                                                                                                           <  2:750687/142‑1 (MQ=255)
                           ggAGGGCGCCCGAAAAATGGTGAGTTTTCCCCTGCGCTTACGGAAACAGGTCGCCGCGCGTCCCGAAACGGTGGTGTTTCTCGACCCTTAAAC‑CCCC‑TATCAGG‑ATAATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGaa                                                                                                                     <  1:841897/151‑1 (MQ=255)
                             aCGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                                    >  2:1130973/1‑150 (MQ=255)
                              cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                  >  1:778427/1‑151 (MQ=255)
                                 cGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                               >  2:534256/1‑151 (MQ=255)
                                    ccGAAAATTGGGGGGTGATCCCCGGCGCGTACTGGAACGGGCCCCCGCGCGTCCCGAACCCGTGGTTTATCCCGACCCTTAAAC‑CCGC‑TATCAGG‑ATCATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                             <  2:247769/150‑1 (MQ=255)
                                    ccGAAAAGTGGTGGGTGATTCCCTGCGCGTACTGGAGCGGTTCGCCGCGCGTACCGAACCGGTGGTGTATCTCGAACCGTTAAC‑CCGC‑TATCAGG‑ATCATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                               <  1:1212083/148‑1 (MQ=255)
                                    ccGAAAAAGGGTGAGTGATCCCCTGCGCGTACTGGAACGGGTCGCCGCGCGTCCCGGACCCGTGGTTTATCCCGACCCTTTAAC‑CCCC‑TATCAGG‑ATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                                <  1:1092273/147‑1 (MQ=255)
                                    cagaaaagTGGTGGGTGCTCCCCGGCGCGTCCGGGAACGGGTCGCCGCGCGCCCGGGACCCGTGGTGTACCTC‑AACCCTTAAACCCGC‑TATCAGG‑ATAATAAGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                                <  2:258645/143‑1 (MQ=255)
                                       aaaagTGGTGAGTGCTCCCCTGCGCTTCCTGAACCAGGCCCCCGCCCGTACCGAACCCGGGGTTTCCCTC‑AACCCTTAAACCCGC‑TATAAGC‑ATCATAGAATCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                         <  1:617373/150‑1 (MQ=255)
                                          aaTGGTGGGTGATCCCCTGCGCGTACTGAAACAGGTCGCCGCGCTTACCGAACCCGTGGTTTATCCCGACCCTTTAAC‑CCGC‑TATAAGG‑ATCATATAATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                      <  1:118051/151‑1 (MQ=255)
                                             ggtggttgttccccGGCGCGTACTTGAACAGGCCCCCGCGCGTCCGGGACCCGTGTTGTACCCCGACCCGTAAAC‑CCGC‑TATAAGA‑ATAATATAACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                                   <  2:1065507/140‑1 (MQ=255)
                                                 agtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                  >  2:1010621/1‑148 (MQ=255)
                                                 agtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                  >  1:1159760/1‑148 (MQ=255)
                                                  gttatccccTGGGCGTCCTGAAACGGGTCGCGGCGCGTCCCGAACCCGGGGTTTATCCCGACCCTTTAAA‑CCGC‑TACAAGA‑ATCATAGGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                              <  1:821789/145‑1 (MQ=255)
                                                      ttCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                          >  2:74233/1‑151 (MQ=255)
                                                      ttCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                          >  2:57106/1‑151 (MQ=255)
                                                        cccTGCGCTTACGGAAACGGGTCGCCGCGCGTACCGGACCCGGGGTTTATCTCGACCCTTTAAC‑CCGC‑TATCAGG‑ATAATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                        <  1:536143/151‑1 (MQ=255)
                                                              gcgTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                  >  1:542088/1‑151 (MQ=255)
                                                                gTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                >  1:136685/1‑151 (MQ=255)
                                                                   cTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             >  1:947244/1‑151 (MQ=255)
                                                                    tGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             >  2:416624/1‑150 (MQ=255)
                                                                    tGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             >  1:1124728/1‑150 (MQ=255)
                                                                    gggaacggtccgccgcGCGTCCGGGACCGGTGGTGTTCCCCGACCGGTTAAC‑CCGC‑TATCAGA‑ATCATATGACAACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                                <  2:830020/137‑1 (MQ=255)
                                                                     ggAACAGGCCGCCGGGCGTCCGGCACCCGTGGTTTATCCCGACCCGTAAAC‑CCGC‑TATCAGG‑TTAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             <  2:934391/149‑1 (MQ=255)
                                                                       aaCAGGCCGCCGCGCGTCCCGAACCGGTGGTGTACCCCGACCCGTAAAC‑CCGT‑AACCAGG‑ATCATATGAACAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                           <  2:426984/149‑1 (MQ=255)
                                                                              cgccgcGCGTACCGGAACGGTGGTGTATCTCGAACCGTTAAC‑CCGC‑TATCAGA‑ATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                                  <  1:120246/151‑1 (MQ=255)
                                                                                     cgTCCCGGACCCGTGGTGTCTCTCG‑ACCCTTAAA‑CCCCTTATAAGA‑ATAATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                 <  2:432603/145‑1 (MQ=255)
                                                                                     cgTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  2:907413/1‑149 (MQ=255)
                                                                                     cgTACCGAACCCGTGGTGTTTCTCGAACCTTAAAA‑CCGC‑TATCAGG‑ATAATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                              <  1:804298/148‑1 (MQ=255)
                                                                                        cccGGAACCGTGTTGTATCTCAACCCTTTAAC‑CCGT‑TACCAGG‑ATCATATATCCAACCCCCTCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                            <  1:1145615/146‑1 (MQ=255)
                                                                                           gaacccgggttgtATCTCCACCCTTTAAC‑CCCC‑TATAAGG‑ATCATAGAATAAACCCCCTCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                           <  1:534541/135‑1 (MQ=255)
                                                                                             acccGGGGTGTATCTCGACCCTTAAAA‑CCCC‑TATCGGG‑ATCAAATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                   <  2:998090/149‑1 (MQ=255)
                                                                                             aaCCGTGGTGTACCCCGACCCGTAAAA‑CCGC‑TATCAGG‑ATCATATGATAAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                   <  1:40956/151‑1 (MQ=255)
                                                                                              aCCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                  >  1:341465/1‑151 (MQ=255)
                                                                                                cGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                >  2:101831/1‑151 (MQ=255)
                                                                                                 gtggttTATCTCGAACCGTAAAC‑CCCC‑TATAAGA‑ATAATATGATCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                   <  2:1022698/147‑1 (MQ=255)
                                                                                                 gtggttTACCCCGAACCTTAAAA‑CCGC‑TATCAGT‑ATAATATGATCAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                 <  2:902863/149‑1 (MQ=255)
                                                                                                 gtggtgTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                               >  1:1180087/1‑151 (MQ=255)
                                                                                                 gtggtgATTCCCTAACCGTTAAA‑CCCC‑AATCAGA‑ACAATATAAACACCTCCCTCGCGCCTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                               <  2:705183/151‑1 (MQ=255)
                                                                                                   ggtgTATCCCGAACCTTTAAC‑CCGC‑TATAAGG‑ATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCa                                              <  1:205398/150‑1 (MQ=255)
                                                                                                     tgtATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                           >  1:709763/1‑151 (MQ=255)
                                                                                                         tcccaaaCCTTAAAA‑CCGC‑TATCAGG‑ATCATATGAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                           <  2:407296/142‑1 (MQ=255)
                                                                                                          ctcGACCCTTAAAC‑CCGT‑TATCGGA‑ATCATATAATCAACCCCCCCCCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                          <  2:268959/147‑1 (MQ=255)
                                                                                                          cccaacccGTTAAC‑CCGC‑TATCAGG‑ATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                        <  1:349632/143‑1 (MQ=255)
                                                                                                              aaCCTTAAAA‑CCGC‑TATGGGG‑ATCATATGATCAACACCCCCCCGGACCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCa                                     <  1:597296/148‑1 (MQ=255)
                                                                                                               aCCCTTAAACCCGC‑TATCAGG‑ATAATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                    <  2:872440/149‑1 (MQ=255)
                                                                                                                  gTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                              >  1:877833/1‑151 (MQ=255)
                                                                                                                     aaa‑CCGC‑TAAAAGA‑ATCATATGAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                             <  1:781721/149‑1 (MQ=255)
                                                                                                                      acccccc‑TACCAGG‑AACATATAACAACCCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                           <  2:419886/145‑1 (MQ=255)
                                                                                                                      aa‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                          >  2:331075/1‑151 (MQ=255)
                                                                                                                       a‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                          >  1:752699/1‑150 (MQ=255)
                                                                                                                          cGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       >  2:475869/1‑151 (MQ=255)
                                                                                                                          cGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       >  1:815492/1‑151 (MQ=255)
                                                                                                                          cGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       >  2:236962/1‑151 (MQ=255)
                                                                                                                                tCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                  >  1:756260/1‑151 (MQ=255)
                                                                                                                                   gg‑ATAATAGAATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt               <  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)
                                                                                                                                                    |                                                                                                                                                   
ACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGG‑ATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  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: ‑