Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,381,7280AC45.7% 108.1 / 23.2 81G197G (GGT→GGGfrdAanaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (34/10);  new base C (0/37);  total (34/47)
Fisher's exact test for biased strand distribution p-value = 3.82e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGG  >  NC_000913/4381590‑4381877
                                                                                                                                           |                                                                                                                                                     
caGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGcc                                                                                                                                            <  2:638434/150‑1 (MQ=255)
caGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCa                                                                                                                                           <  1:1354826/151‑1 (MQ=255)
    aTACCGGAACCTGGCAGACCGGTTGGGCGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCa                                                                                                                                           >  2:1051362/1‑147 (MQ=255)
     tACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGc                                                                                                                                      >  2:1073268/1‑151 (MQ=255)
          gAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATa                                                                                                                                  >  2:1080101/1‑150 (MQ=255)
           aaCCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAAc                                                                                                                                >  1:530263/1‑151 (MQ=255)
             ccTGGCAGCCCGGTGGGGTGATACTAAACAAATTCCATGCCACGCGGCGGACCGCCGGGGCTTAGCGCAATCCCCATCCCGTCCCCGGTAACGAGGCCGCCGTTGGTGTTGAACCGATAAACGCGCCCCGCACCGCCAGTAGCCATaa                                                                                                                                 <  1:828660/148‑1 (MQ=255)
              cTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGAc                                                                                                                             >  2:464108/1‑151 (MQ=255)
               tGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGAcc                                                                                                                            >  1:1397860/1‑151 (MQ=255)
                gggagcccGGTTGGGTGATACTGAACAAATCCCATGTCACGCGGGGGACCGCCGTGCCTTACCGCCATACCCATCCCGCCCCGGTTAACGATGCCGCCGTTGGTGTTTTACAGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGAc                                                                                                                             <  2:1190934/143‑1 (MQ=255)
                      cgggTTGGGTGATTCTGAACGAATCCCTTGTCACGCAGCGGAACGCCGTGGCTTCGCCCCATCCCCACCCCCTCCCCGGTACCGATGCCGCCGTTGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                                                                                                     <  1:1063026/149‑1 (MQ=255)
                      ccGGTTGGGTAATCCTAAACGATTCCCCTGCCCCGCAGGGAACCGCCGGGGTTTACCCCCACCCCCATCCCTCCCCCGTTACCGATCCCGCCTTTGGTGTTGTAACAAAAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGtt                                                                                                                       <  2:649080/149‑1 (MQ=255)
                        ggTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGc                                                                                                                    >  2:717467/1‑150 (MQ=255)
                        ggTTGGGTGATACTGAACAAATTCCATGTCCCGCGGCGGACCGCCGTGGCTTAGCGCCATCCCCATCCCGTCCCCGGAACCGATGCCGCCGTTGGTGTTGTAACGATAACCGCGCCCCGCCCCGCCAGTAGCCATAACGACcgc                                                                                                                          <  2:151261/144‑1 (MQ=255)
                         gTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCAc                                                                                                                  >  1:124041/1‑151 (MQ=255)
                            gggTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                                                                                                >  2:466010/1‑150 (MQ=255)
                            gggTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                                                                                                >  2:1538913/1‑150 (MQ=255)
                            gggTGATACTGAACGAATCCCATGCCCCGCAGCGGACCGCCGTGGTTTAGCCCAATCCCCATCCCGCCCCCGTAAACGATGCCGCGGTTGGTGTTGTAACAATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACg                                                                                                                 <  1:86070/149‑1 (MQ=255)
                            gggTGATACTAAACGAATCCCATGTCCCGCAGCGGAACGCGGTGGTTTAGCGCCATCCCCATCCCGTCCCCGTTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGa                                                                                                               <  1:444650/151‑1 (MQ=255)
                            gggTAATCCTAAACAAATTCCATGTCCCGCAGGGGACCGCCGTGGTTTAGCGCCATCCCCACCCCGCCCCGGGTAACGAGCCCGCCTTTGGTGTTGTAACAATAACCGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGa                                                                                                               <  1:627431/151‑1 (MQ=255)
                              gTGATACTAACCGATTTCCAGTCCCCGCGGCGGACCGCCGTGGCTTAGCGCCATCCCCACCCCGCCACCGGTAACGATGCCCCCGTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATc                                                                                                             <  2:375963/151‑1 (MQ=255)
                               tGAAACGGACCGAATTCCATGCCACGCAGGGGAACGCCGTGGTTTAGCGCAATCCCCATCCCGTCCCCGGAAACGAGGCCGCCGTTGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATc                                                                                                             <  2:1265071/150‑1 (MQ=255)
                                  tACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGc                                                                                                          >  1:1387988/1‑150 (MQ=255)
                                   aCTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCAc                                                                                                        >  2:1490390/1‑151 (MQ=255)
                                       accGAATCCCATGTCCCGCAGCGGACCGCCGGGGCTTAGCGCCATCCCCATCCCGTCCCCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCAcc                                                                                                       <  1:1463036/146‑1 (MQ=255)
                                           aaTTCCATGTCCCGCAGCGGAACGCGGTGGTTTAGCGCCATCCCCATACCGTCACGGGAAACGATGCCGCCGTTGGTTTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGt                                                                                                  <  1:1441136/149‑1 (MQ=255)
                                            aTTCCATGTCCCGCGGCGGACCGCCGGGGTTAAGCGCCATCCCCACCCCGTCCCGGGTACCGATGCCGCGTTTGGTGTTGAAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAg                                                                                                     <  2:805775/145‑1 (MQ=255)
                                             ttCCATGTCCCGCAGGGGAACGCGGTGGCTTAGCGCCATCCCCATCCCGCCCCCGGAAACGATGCCGCCGTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGc                                                                                                <  2:1009096/149‑1 (MQ=255)
                                             ttCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCt                                                                                              >  2:897487/1‑151 (MQ=255)
                                                cATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTcc                                                                                           >  2:365502/1‑151 (MQ=255)
                                                 aTGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTcc                                                                                           >  1:1039274/1‑150 (MQ=255)
                                                     cACGCAGCGGAACGCGGTGGTTTAGCCCCATCCCCATCCCGTCACCGGTAACGATGCCGCCGTTGGTTTTGTAACAATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat                                                                                      <  1:1219226/151‑1 (MQ=255)
                                                       cGCAGGGGACCGCCGTGGTTTAGCGCCATCCCCATCCCGTCACCGTTAACGAGGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATg                                                                                     <  2:344277/150‑1 (MQ=255)
                                                          aGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTc                                                                                  >  2:1176783/1‑150 (MQ=255)
                                                            cGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCAtt                                                                               >  2:39524/1‑151 (MQ=255)
                                                                aCGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCt                                                                            >  1:477742/1‑150 (MQ=255)
                                                                   ccGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTAcc                                                                         >  1:775779/1‑150 (MQ=255)
                                                                     gTGGGTTAGCGCATTACCCAACCCGTCACCGGAACCGTTGCCGCGGTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCa                                                                        <  2:1459566/149‑1 (MQ=255)
                                                                     gTGGCTTGGCGCCATCCCCAACCCGCCACCGGTAACGATGCGGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAg                                                                       <  2:59853/150‑1 (MQ=255)
                                                                      tGGATTAGCGCCATCCCAATACCGCGACCGGTACCGGTCGGGCGGTGGTGTGTGTAACAAAAAGCCGGACCGGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt                                                                                   <  2:931967/137‑1 (MQ=255)
                                                                       ggCTTAGCGCCATCCCCATACCGTCACCGTAAGCGATGCCGCGGTGGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGcc                                                                    <  1:1189988/151‑1 (MQ=255)
                                                                          ttcgccccacccccatcccGCCACCGGTAACGATCCCGCCGTTGGTTTTTTAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCt                                                                            <  2:694919/129‑1 (MQ=255)
                                                                           tAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcgcg                                                                 >  2:286570/1‑150 (MQ=255)
                                                                                ccatccccatccCGCCCCCGGTAACGATGCCGCCTTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGaa                                                               <  1:867409/147‑1 (MQ=255)
                                                                                caatccccatcccGTCCCCATTAACGTTGCCGCGGTTGGTGTTGTAACGAAAACCGCGACCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACAt                                                            <  1:552174/145‑1 (MQ=255)
                                                                                  atccccatcccGTCCCCGGTAACGATGCCGCCGTTGGTGTTGTAACAAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGa                                                          <  1:1201956/147‑1 (MQ=255)
                                                                                     cccatcccGTCCCGGGAACCGAGGCCGCCGTTGGTGTTGAAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGa                                                          <  2:1536346/147‑1 (MQ=255)
                                                                                     cccatcccGTCCCCGGTAACGATGCCGCCGTTGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGAc                                                         <  1:103563/148‑1 (MQ=255)
                                                                                     cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGAc                                                         >  2:787232/1‑148 (MQ=255)
                                                                                     cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcat                                                      >  1:798906/1‑151 (MQ=255)
                                                                                      ccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcatc                                                     >  2:136137/1‑151 (MQ=255)
                                                                                          cccGTCCCCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCAtcatca                                                 <  1:1178300/150‑1 (MQ=255)
                                                                                            cGTCCCCGGTAACGATGCCGCCGTGGGTTTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcatcat                                                   <  2:164322/147‑1 (MQ=255)
                                                                                              tCACCGGTAACGATTCCGCCGTTGGAGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAAcc                                              >  2:733178/1‑150 (MQ=255)
                                                                                               ccccGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAg                                            <  1:175247/149‑1 (MQ=255)
                                                                                               cACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAAcc                                              >  1:802867/1‑149 (MQ=255)
                                                                                                 ctgataaggtggccgcgGTTGTTTTTGAAACGAAAAACGCGCCCCGCTCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGa                                           <  2:1039520/139‑1 (MQ=255)
                                                                                                  cGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAat                                         >  1:623574/1‑151 (MQ=255)
                                                                                                   gaaaaCGAGGCCGCGGTTGGTGTTGTAAGAAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGa                                           <  2:919593/145‑1 (MQ=255)
                                                                                                       aCGATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCa                                    <  1:533518/151‑1 (MQ=255)
                                                                                                         gctGCGGCCGTTGGTTTTGAAAGGATAACGGCGCCCGGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCa                                    <  2:1001255/147‑1 (MQ=255)
                                                                                                         gATGCCGCCGTTGGTGTTGTAACGATAACCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGc                                  <  1:116153/151‑1 (MQ=255)
                                                                                                           gtccgcctTTGTTGTTGTAACAAAAAGCGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGc                                  <  2:1166114/147‑1 (MQ=255)
                                                                                                             ccgccgTTGGTTTTGTACCAATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAg                                   <  2:1506679/146‑1 (MQ=255)
                                                                                                             ccgccgTTGGTGTTTTAACAAAAACCGCGCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCAc                                <  2:1186457/149‑1 (MQ=255)
                                                                                                                ccgTTGGTGTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGaa                             <  1:497086/149‑1 (MQ=255)
                                                                                                                   ttGGTGTTGTAACGAAAACCGCGCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATg                          <  1:732632/149‑1 (MQ=255)
                                                                                                                       ttttgtAACAATAAACGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCg                      <  1:1135117/147‑1 (MQ=255)
                                                                                                                       tgttgtAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTc                    >  2:623947/1‑151 (MQ=255)
                                                                                                                          ttaaaaCGAAAACCGCGCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCa                   <  2:475821/146‑1 (MQ=255)
                                                                                                                              aaCGATAAAGGCGCCCCCCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAAc               <  2:1028801/150‑1 (MQ=255)
                                                                                                                              aaCGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCaaa                 <  1:849800/148‑1 (MQ=255)
                                                                                                                              aaCGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACg              <  1:1245732/151‑1 (MQ=255)
                                                                                                                              aaCGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGCACCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATCACCATCATCCACCACAACATCGAGCACGAAATGTCAGTCAAACCg              >  2:106250/1‑151 (MQ=255)
                                                                                                                              aaCGAAAACCCCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCaaaa                <  1:331924/149‑1 (MQ=255)
                                                                                                                                 gATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGACGCGAACATGACCATCATCAACCAGAATATCCCGCACGAAATGTTCGTCAAAACGCTg           >  2:1295350/1‑151 (MQ=255)
                                                                                                                                  aTAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTgg          >  2:744263/1‑151 (MQ=255)
                                                                                                                                    aaaCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGAt        >  1:60992/1‑151 (MQ=255)
                                                                                                                                    aaaCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGAt        >  2:962909/1‑151 (MQ=255)
                                                                                                                                         cgACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCg   >  2:1335663/1‑151 (MQ=255)
                                                                                                                                          gACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCgg  >  1:157111/1‑151 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
CAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGG  >  NC_000913/4381590‑4381877

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: