Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,381,7280AC42.2% 150.5 / 21.4 91G197G (GGT→GGGfrdAanaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (44/8);  new base C (0/38);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 7.61e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGG  >  NC_000913/4381578‑4381878
                                                                                                                                                          |                                                                                                                                                      
aCCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGa                                                                                                                                                        <  1:410355/151‑1 (MQ=255)
aCCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGa                                                                                                                                                        <  1:119631/151‑1 (MQ=255)
 ccTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGAc                                                                                                                                                       <  1:60842/151‑1 (MQ=255)
    tCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCg                                                                                                                                                    >  1:320260/1‑151 (MQ=255)
         catcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCg                                                                                                                                               >  1:10797/1‑151 (MQ=255)
         catcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCg                                                                                                                                               >  2:1033756/1‑151 (MQ=255)
            caGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGcc                                                                                                                                             <  1:152682/150‑1 (MQ=255)
            caGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGcc                                                                                                                                             <  2:1120435/150‑1 (MQ=255)
                      gAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATa                                                                                                                                   >  2:129502/1‑150 (MQ=255)
                          cTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGAc                                                                                                                              >  1:1272141/1‑151 (MQ=255)
                            gggaggcgggTTGGGTGACACTAAACAATTTCCGTGTCCCGGGGCGGAACGGGGTGGTTTAGCGCCAACCCCATCCCTGCGGCGGTAAGGGTCCGCCGGTTGTTGTTGTACAGAAAAACGCGCCCGGCTCCGCCAGTAGCCATAACGAcc                                                                                                                             <  2:1527456/142‑1 (MQ=255)
                            ggaaGACCGGTGGGTGATTACTAAAGGATTTCTTTGTCAGGCAGGGGAACGCGGGGCCTTGGCGCCACCCCCATCCGTCCCCGGGAAACGACGCCGCCGTTGGTGTTGTAGGGAAAACCGCGACCCGCCCCGCCAGTAGCCATAACGACcg                                                                                                                            <  2:675060/148‑1 (MQ=255)
                               aGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGt                                                                                                                         >  1:473167/1‑151 (MQ=255)
                                    ggTTGGGTGATACTGAACGAATCCCATGTCACGCAGCGGAACGCCGTGGTTTAGCGCCATCCCCATCCCGTCCCCGGTAACGAGGCCGCCGTTGGTGTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                                                                                                    <  1:1195565/151‑1 (MQ=255)
                                      ttGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACg                                                                                                                  >  2:864044/1‑151 (MQ=255)
                                      ttGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACg                                                                                                                  >  1:877237/1‑151 (MQ=255)
                                        ggttGATATTAAACGAATTCCATCCCCCACACGGGAACCAGGTGTTTTCGCCCCATCCCCATCCGGTCACCGTTAAGGTGCCCCCCGTTGGTGTTGTACCAATAAAGGCGGCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGc                                                                                                                     <  1:231054/143‑1 (MQ=255)
                                        gggTGATACTGAACAAATTCCATGTCCCGCAGCGGACCGCCGTGGCTAAGCGCCATCCCCATACCGTCCCCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                                                                                                 <  1:1346633/150‑1 (MQ=255)
                                        gggTGATACTAAACTAATTCCATGTCACGCAGCGGACCGCCGTGGTTTAGCGCCATCCCCACCCCGTCCCCGGTAACGACGCCGCGGTTGGTGTTGTAACAATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                                                                                                 <  1:664371/150‑1 (MQ=255)
                                        gggTGATACTAAACGAATTCCATCTCACGCAGCGGACCGCGGTGGCTTAGCGCCATCCCCACCCCGTCCCCGATAACGATGCCCCCGTTGGTGTTGTAACGATAAACGCGACCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                                                                                                 <  2:65121/150‑1 (MQ=255)
                                         ggataTACTAAACAAATTCCATGTCACGCAGAGCAGCGCGTGGGCGCAGCCCCATACCCACACCGCCCGGGGAAACGACGCCGCCGGGGGTGTTGTAAGGAAAAGCGGGCCCGGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGAt                                                                                                               <  2:746214/147‑1 (MQ=255)
                                           tGATCCTAAACAATTCCCATGCCCCGCAGCGGACCGCGGTGGTTTAGCCCCATCCCCATCCCTCCGCGGGTAACGAGGCCGCCTTTGGTTTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATc                                                                                                              <  2:279318/150‑1 (MQ=255)
                                            gATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTg                                                                                                            >  1:107447/1‑151 (MQ=255)
                                             gttACTAAACGATTCCCAGGTCCCGCAGGGGAAGGCGGTGGTTGAGCCCCATCCCCATCCCGCCCCGGTTACGGAGGCCGCGGTTGGTTTTGTAACAAAAACGGGGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCt                                                                                                             <  2:302959/148‑1 (MQ=255)
                                                acgaaaCAAATTCCTTGCCCCGCGGGGGACCGCGGTGGTTTAGCCACTTCCCCAGCCCGTCCCGGGAAAGGTTGCGGCTGTGGGTGTTGTACCAAAAAGCGAGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCa                                                                                                          <  2:1037875/146‑1 (MQ=255)
                                                  tGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCAcc                                                                                                        >  1:774214/1‑150 (MQ=255)
                                                   gAACGAATTCCATGTCCCGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTTTAACGATAACCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCa                                                                                                       <  1:353699/150‑1 (MQ=255)
                                                   gAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCAcc                                                                                                        >  1:86254/1‑149 (MQ=255)
                                                    accaaaTTCCAGGTCACGCGGCGGAACGCTGGGGCTTCGCCCCATCCCCATCCCGCCCCGGGAAACGATGCCGCTTTTGGTGTTGTACCAAAAACCGCACCCCGACCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAg                                                                                                      <  2:1094691/146‑1 (MQ=255)
                                                      cGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGt                                                                                                   >  2:597842/1‑151 (MQ=255)
                                                        aaTTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGc                                                                                                 >  1:512843/1‑151 (MQ=255)
                                                        aaTCCCATGCCACCAAGCGGACCGCCGGGGCTTAGCCCCATCCCCATACCGTCACCGGTAACGATGCCGCCGTTGTTTTTGTAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGc                                                                                                 <  1:975000/151‑1 (MQ=255)
                                                         aTTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGcc                                                                                                >  2:813158/1‑151 (MQ=255)
                                                           tCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCtt                                                                                              >  1:245720/1‑151 (MQ=255)
                                                                gTCACGCGGCGGACCGCGGTGGTTTAGCCCCATCCCCATCCCGTCCCGGGTAACGATGCCGCGGTGGGTTTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTc                                                                                             <  2:670365/147‑1 (MQ=255)
                                                                  cccgcgGCGGACCGCCGTGGTTTAGCCCCACCCCCATCCCGTCCCGGGTAAGGATCCCGCCGTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat                                                                                       <  2:795429/149‑1 (MQ=255)
                                                                  cACGCAGCGGAGCGCCGTGGCTGGTCCCCATCCCCACCCCGTCGCCGGTAACGGTGCCGCTGTTGGTGTTGAAACAAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatca                                                                                        <  1:1023267/150‑1 (MQ=255)
                                                                  cACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat                                                                                       >  2:1323897/1‑151 (MQ=255)
                                                                        gCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCAt                                                                                 >  2:1305790/1‑151 (MQ=255)
                                                                         cGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCAtt                                                                                >  1:1377079/1‑151 (MQ=255)
                                                                          ggaccccggTGGTTTAGCGCCATCCCCACCCCGTCCCGGGTAAGGATGCCGCCGTTGGTTTTGTAAGAATAACGGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTc                                                                                   <  2:233768/139‑1 (MQ=255)
                                                                             aCGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTg                                                                               >  1:813317/1‑148 (MQ=255)
                                                                                ccGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTAcc                                                                          >  2:69449/1‑150 (MQ=255)
                                                                                 cGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTAcc                                                                          >  1:636004/1‑149 (MQ=255)
                                                                                  gTGGTTTAGCTCCACCCCCAATCCATCCCCGGTAGCGATCCCGCCTTTGTTTTTATATCAAAAAGCGCGCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCa                                                                         <  2:240571/149‑1 (MQ=255)
                                                                                  gTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAgg                                                                       >  1:1198705/1‑151 (MQ=255)
                                                                                    ggtttAGCGCCATCCCCACCCCGTCCCCGGTAACGAGGCCGCCGTTGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGcc                                                                     <  1:1158608/148‑1 (MQ=255)
                                                                                    ggCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGcc                                                                     >  2:772006/1‑151 (MQ=255)
                                                                                       ttAGCCCCATCCCCATCCCTTCCCGTGAAACGTTGCCGCGGTTTGTTTTGAACCAAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcg                                                                    <  2:1212339/149‑1 (MQ=255)
                                                                                        tAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcgcg                                                                  >  1:1015279/1‑150 (MQ=255)
                                                                                         aGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGa                                                                 >  2:1447668/1‑150 (MQ=255)
                                                                                         aGCCCCATCCCCAAACCGCCCCCGGTAACGATGCCGCGGTTGGTGTTGTAACAAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGcc                                                                     <  2:119142/146‑1 (MQ=255)
                                                                                            gcCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACAt                                                             >  1:955054/1‑151 (MQ=255)
                                                                                             ccatccccatccCGTCCCGGGAAACGATGCCGCCGTTGGTGTTGTAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcgc                                                                   <  1:1527904/144‑1 (MQ=255)
                                                                                             caatccccatcccGCCCCCGGAACCGGTGCCCCCGTGGGTGTTGTAACGAAAAACGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATg                                                            <  2:473148/146‑1 (MQ=255)
                                                                                               atacccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGAc                                                          >  2:1499918/1‑151 (MQ=255)
                                                                                               atacccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGAc                                                          >  2:676418/1‑151 (MQ=255)
                                                                                                tacccataccGTCACCGGTAACGATGCCGTCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGAc                                                          >  2:839118/1‑150 (MQ=255)
                                                                                                   ccatcccGTCACCGGAAACGATCCCGCGGTTGTTGTTGTAACGATAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcatc                                                      <  1:738300/151‑1 (MQ=255)
                                                                                                   ccatcccGCCACGGGTAACGATGCCGCCGTTGGTGTTGAAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcatc                                                      <  1:1124910/151‑1 (MQ=255)
                                                                                                       cccgcccccGGAAACGATGCCGCCTTTGGTTTTGTAACGAAAACCGCGCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACcatc                                                      <  1:1472888/140‑1 (MQ=255)
                                                                                                             aCCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAAcc                                               >  2:598551/1‑148 (MQ=255)
                                                                                                              ccGTTAACGTTGCCCCCGTGGGTGTTGAAACAAAAACCGCCCCCCCCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAAc                                                <  2:1474273/146‑1 (MQ=255)
                                                                                                              ccGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGa                                            >  2:1212210/1‑150 (MQ=255)
                                                                                                               cGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAat                                          >  1:1506231/1‑151 (MQ=255)
                                                                                                                 gtaaaCGATGCCTCCGTTGGTGTTAAAACGACAACCGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAat                                          <  2:4500/147‑1 (MQ=255)
                                                                                                                 ggaaaCGATGCCGCGTTGGGTGTTGTAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAat                                          <  1:1382154/147‑1 (MQ=255)
                                                                                                                  gaaccGAGGCCGCTGTTGTTGTTTTAACAATAACCCCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAatat                                        <  1:1326788/147‑1 (MQ=255)
                                                                                                                            ccgcggTTGGTTTGGTAACGAAAACCGCACCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTACTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCa                                  <  2:291089/148‑1 (MQ=255)
                                                                                                                            ccgccgTTGGTGTTGTAACGATAAACGGGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCa                                  <  2:1187156/148‑1 (MQ=255)
                                                                                                                             cgccgTTGGTGTTGTAACGATAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGa                               <  1:386649/150‑1 (MQ=255)
                                                                                                                             ccgcgtttGGTTTTGTAACGATAACCGCCCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCAc                                 <  2:9816/143‑1 (MQ=255)
                                                                                                                              gccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGaa                              >  1:1370026/1‑150 (MQ=255)
                                                                                                                               ccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAAt                            >  2:537041/1‑151 (MQ=255)
                                                                                                                                cgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAAt                            >  1:1406662/1‑150 (MQ=255)
                                                                                                                                   ttGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGt                          >  2:124864/1‑150 (MQ=255)
                                                                                                                                       tgttttAACGAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGt                      <  2:774183/150‑1 (MQ=255)
                                                                                                                                         ttgtAACGAAAAACGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCa                    <  1:129508/150‑1 (MQ=255)
                                                                                                                                             aaCGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACg               >  2:992542/1‑151 (MQ=255)
                                                                                                                                             aaCGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACg               >  2:1432700/1‑151 (MQ=255)
                                                                                                                                             aaCGATAAACGAGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACg               <  1:64760/151‑1 (MQ=255)
                                                                                                                                             aaCGAAAACCGCGCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACg               <  1:662282/151‑1 (MQ=255)
                                                                                                                                             aaCGAAAAACGCCCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCaaaa                 <  2:195918/149‑1 (MQ=255)
                                                                                                                                              aCGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGc              >  1:493986/1‑151 (MQ=255)
                                                                                                                                                gATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTg            >  1:212814/1‑151 (MQ=255)
                                                                                                                                                 aaaaccgcgCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTg            <  2:1391947/145‑1 (MQ=255)
                                                                                                                                                 aaaaaCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCaaa                  <  2:613053/142‑1 (MQ=255)
                                                                                                                                                 aTAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTgg           >  2:297353/1‑151 (MQ=255)
                                                                                                                                                  tAAACGCGACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTgg           >  2:696900/1‑150 (MQ=255)
                                                                                                                                                  tAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTgg           >  2:968346/1‑150 (MQ=255)
                                                                                                                                                   aaaCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGAt         <  1:1471153/151‑1 (MQ=255)
                                                                                                                                                   aaaCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGAt         <  2:1485738/151‑1 (MQ=255)
                                                                                                                                                         gACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCCCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCgg   >  2:1438882/1‑151 (MQ=255)
                                                                                                                                                          aCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCggg  >  1:460773/1‑151 (MQ=255)
                                                                                                                                                          |                                                                                                                                                      
ACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCGCGAACATGACCATCATCAACCAGAATATCCAGCACGAAATGTTCGTCAAAACGCTGGATCTGCGGG  >  NC_000913/4381578‑4381878

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: