Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,381,7280AC46.7% 99.9 / 21.2 75G197G (GGT→GGGfrdAanaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (31/9);  new base C (0/35);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 2.52e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGAT  >  NC_000913/4381579‑4381871
                                                                                                                                                       |                                                                                                                                               
ccTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGAc                                                                                                                                                >  2:604568/1‑151 (MQ=255)
ccTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGAc                                                                                                                                                >  2:405751/1‑151 (MQ=255)
    cGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCg                                                                                                                                             <  1:707137/150‑1 (MQ=255)
         atcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGc                                                                                                                                       >  1:939832/1‑151 (MQ=255)
           caGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGcc                                                                                                                                      <  1:475283/150‑1 (MQ=255)
             ggATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGt                                                                                                                                   >  2:209773/1‑151 (MQ=255)
               aTACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGt                                                                                                                                   >  1:652396/1‑149 (MQ=255)
                    ggACCCGGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGGGGCTTAGCGCCATCCCCATACCGCCACCGGTACCGATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGcc                                                                                                                               <  1:1157939/148‑1 (MQ=255)
                    ggAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCAt                                                                                                                             >  1:857603/1‑150 (MQ=255)
                         cTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGAc                                                                                                                       >  1:1047705/1‑151 (MQ=255)
                          tGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGAcc                                                                                                                      >  2:962825/1‑151 (MQ=255)
                          tGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGAc                                                                                                                       >  1:806906/1‑150 (MQ=255)
                            gCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACcgc                                                                                                                    >  1:370015/1‑151 (MQ=255)
                                 ccGGTTGGGTGATCCTGAACAAATTCCATGCCCCGCGGCGGACCGCCGTGGCTTGGCGCCATCCCCATCCCGCCCCGGGTAACGATCCCGCCTTTGGTGTTAAACCGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGtt                                                                                                                 <  1:863188/149‑1 (MQ=255)
                                 ccGGTGGGGTGATACTGAACAATTTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATCCCCATCCCTTCCCCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                                                                                               <  2:709889/151‑1 (MQ=255)
                                   ggTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGc                                                                                                              >  2:368288/1‑150 (MQ=255)
                                     ttGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACg                                                                                                           >  2:521118/1‑151 (MQ=255)
                                      tGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                                                                                          >  1:425728/1‑151 (MQ=255)
                                         gTGATACTGAACGAATTCAATGTCACGCAGCGGAACGCCGTGGCTTAGCGACATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgca                                                                                                         >  1:1044061/1‑147 (MQ=255)
                                             tACTGAACGAATTCCATGTCCCGCGGCGGAACGCCGGGGTTTAGCGCCATCCCCACACCGCCACGGGTAACGAGCCCGCGGTTGGTGTTGTAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTg                                                                                                     <  2:728202/149‑1 (MQ=255)
                                                  accGAATTCCATTTCCCCGGGCGGAGCCCCGTGGTTTGGCGCAATCCCCATCCCGTCACGGGTAACGGTGCCGCCTTTGGTTTTGTACCGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGc                                                                                              <  2:160972/149‑1 (MQ=255)
                                                  aaCGAATCCCATTTCACGCAGCGGACCGCCGTGGCTTAGCGCCATACCCATCCCGCCACCGGTAACGAGGCCGCCGTTGGTGTTGTAACGATAACCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCa                                                                                                <  1:330767/149‑1 (MQ=255)
                                                   acaaaTTCCTTGCCGCGCGGCGCAACGCCGTGGCTGAGCAACCTCCCCATTCCATCGCGGTACACGTTCCGCCGGTGGTTATTATAACGACAAGCGCGCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCa                                                                                                <  2:111312/145‑1 (MQ=255)
                                                   aCGAATTCATTTTCCCGCGGCGGACCGCCGGGGCTTAGCGCCCTCCCCATACCGTCACCGGAAACGATGCCGCCTTTGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCAcc                                                                                                 <  2:493682/147‑1 (MQ=255)
                                                    cGAATCCCATGTCACGCAGCGGACCCCCGTGGCTTAGCCCCATCCCCATCCCGTCCCGGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGt                                                                                            <  2:19679/151‑1 (MQ=255)
                                                      aaTTCCATGTCACGCAGCGAAACGCCGTGGCTTAGCCCCATCCCCATACCTTCACCGTTAACAATCCCGCCGTTGGTTTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGt                                                                                            <  2:73067/149‑1 (MQ=255)
                                                          ccATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTc                                                                                      >  1:160246/1‑151 (MQ=255)
                                                            aTGTCACGCGGCGGACCGCGGTGCCTTGGCCCCATCCCCAACCCGCCCCCGTAAACGAGGCCGCCGTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTcc                                                                                     <  1:40923/150‑1 (MQ=255)
                                                             tgccACGCGGCGGACGGCGGTGGCTTAGCGCCATCCCCATCCCGTCACCGGTAACGATGCCGCCGTGGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTcc                                                                                     <  1:342183/146‑1 (MQ=255)
                                                                      gCGGACCGCCGGGGCTTAGCGCCAACCCCACCCCTTCCCCGTAAGCGAGGCCGCCGTTGGTGTTGTAACGATAACGGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCtt                                                                                       <  1:767168/138‑1 (MQ=255)
                                                                           aCGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCt                                                                      >  2:263729/1‑150 (MQ=255)
                                                                            cgcgGTGGCTTGGCGCCATCCCCATCCCTTCCCCGGAAACGATGCCGCCTTGGGTTTTGAAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTa                                                                     <  2:675945/150‑1 (MQ=255)
                                                                             gCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTAcc                                                                   >  1:710290/1‑151 (MQ=255)
                                                                              ccggggtttgccgcCATCCCCATCCCCGCCCGGTTAACGGTGCCGCCTTTGGTTTTGTAGCGATAACGCCCCCCCGCGCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCa                                                                  <  2:821413/140‑1 (MQ=255)
                                                                               cGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTAcc                                                                   >  1:428800/1‑149 (MQ=255)
                                                                                  ggtttggcgcCATACCCATCCGTTCCCCGGTAACGAGCCCGCGGTTGGTTTTGTAACGAAAACCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCa                                                                           <  1:909635/132‑1 (MQ=255)
                                                                                     ttgcccccatccccatccCTTCCCCGTTAACGTTCCCGCCGGTGGTGTTGAAAGAATAACCGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAg                                                                 <  1:639543/139‑1 (MQ=255)
                                                                                     ttAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcgcg                                                           >  2:586697/1‑151 (MQ=255)
                                                                                       aGCGCCATCCCCATCCCGTCACCGGTAACGATGCCGCCGTTGGTTTTGAAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGa                                                          <  1:392638/150‑1 (MQ=255)
                                                                                         cgcCATCCCCATCCCGCCCCGGGTAACTATGCGCCGGTTGGTTTTATACCGAAAAGCCCGCCCGGGACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGaa                                                         <  2:864184/149‑1 (MQ=255)
                                                                                         cgcCATACCCAACCCGTCCCCAGTAACGATCCGGCAGTTGTTTTTGCAACATTAAACGCGCCCCCCTCCGCCAGTAGCCATAACGTCCGCGTTAGCACTGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAAc                                                        <  1:323984/150‑1 (MQ=255)
                                                                                         cgcCACCCCCATCCCCCCCCCGGTAACGATGCCGCGGTTGGTTTTGAACCGATAACCGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAAc                                                        <  1:65620/150‑1 (MQ=255)
                                                                                           ccatccccatccCGTCCCCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACAt                                                      <  2:476382/150‑1 (MQ=255)
                                                                                             atccccatcccTCCACGGGAAAGGAGGCCGCCGTTGGTGTTGAAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGa                                                    <  1:763506/147‑1 (MQ=255)
                                                                                                cccatcccGCCACGGGTAACGATCCCGCGGTTGGTGTTGAAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcat                                                <  1:943935/151‑1 (MQ=255)
                                                                                                 ccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcatc                                               <  2:566706/151‑1 (MQ=255)
                                                                                                   ataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcatcat                                             >  2:396093/1‑151 (MQ=255)
                                                                                                      ccGTCACCGGTAACGAGGCCGCCGTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCaa                                          <  1:180540/151‑1 (MQ=255)
                                                                                                          cACCGGAAACGGTGCCGCGTTTGGTGTTATAACAATAAACGCAACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCa                                       <  2:656135/150‑1 (MQ=255)
                                                                                                            cgggTAACGATCCCGCCTTTGGTGTTTAAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGaa                                    <  2:1157470/149‑1 (MQ=255)
                                                                                                            ccGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGa                                     >  2:882894/1‑150 (MQ=255)
                                                                                                              ggTAACGACCCCGCCGTGGTTTTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAata                                  <  2:606483/151‑1 (MQ=255)
                                                                                                               gTAACGATGCCGCCGTTGGTGTTGTACCGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAatat                                 <  2:1045135/151‑1 (MQ=255)
                                                                                                                  aaggATGCCGCCGTTGGTGTTGAAACAATAACCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAatat                                 <  2:1072520/146‑1 (MQ=255)
                                                                                                                   aCGAGGCCGCGGTTGTTGTGGTAACAAGAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCa                              <  1:641219/151‑1 (MQ=255)
                                                                                                                    cGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATcc                               >  1:160249/1‑149 (MQ=255)
                                                                                                                     gATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGc                            >  1:345129/1‑151 (MQ=255)
                                                                                                                      aTGCCCCCGTTGGTTTTGTACCGAAAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCa                           <  2:416175/151‑1 (MQ=255)
                                                                                                                        cccgccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACg                         >  1:522250/2‑151 (MQ=255)
                                                                                                                         ccgcttttGGTGTTTAAACGATAAGCGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCa                           <  1:347718/142‑1 (MQ=255)
                                                                                                                         ccgcggTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCa                           <  1:107428/148‑1 (MQ=255)
                                                                                                                           gccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGaa                       >  2:158996/1‑150 (MQ=255)
                                                                                                                           gccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGaa                       >  2:1168568/1‑150 (MQ=255)
                                                                                                                                tGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTc                 >  1:62170/1‑151 (MQ=255)
                                                                                                                                    ttttgtAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTc                 <  1:28023/146‑1 (MQ=255)
                                                                                                                                    tgttgtAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTc              >  1:305007/1‑151 (MQ=255)
                                                                                                                                     gttgtAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTc              >  1:1000592/1‑150 (MQ=255)
                                                                                                                                       tgtAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCaa            <  2:242822/150‑1 (MQ=255)
                                                                                                                                       ggtAAAGAAAACCCCGCCCCGCAGCGCCAGTAGCCATAACGACCGCGTTAGCACGGAGCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAGCCAGAATATCCAGCACGAAATGTTCGTCaa            <  1:849444/149‑1 (MQ=255)
                                                                                                                                         tAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAAc         >  1:79101/1‑151 (MQ=255)
                                                                                                                                          aaCGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTTGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACg        >  2:512794/1‑151 (MQ=255)
                                                                                                                                           aCGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGc       <  2:1047739/151‑1 (MQ=255)
                                                                                                                                              aaaaccgcgCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTgg    <  2:939862/146‑1 (MQ=255)
                                                                                                                                              aaaaaCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTgg    <  2:697222/149‑1 (MQ=255)
                                                                                                                                                aaaCGCACCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGAt  <  1:307549/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
CCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGAT  >  NC_000913/4381579‑4381871

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: