Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,073,455:1 +C 6.8% coding (237/1020 nt) epd ← D‑erythrose‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,073,4551.C 6.8% 293.6 / 17.4 74coding (237/1020 nt)epdD‑erythrose‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (33/36);  new base C (2/3);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCACCACGGTAAT  >  NC_000913/3073309‑3073604
                                                                                                                                                   |                                                                                                                                                     
ttACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGc                                                                                                                                                    <  1:508410/150‑1 (MQ=255)
   cTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAAt                                                                                                                                                 <  2:689469/150‑1 (MQ=255)
    tGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATa                                                                                                                                                <  2:738555/150‑1 (MQ=255)
       cAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCATGAACGTTCACTGCAAtcgct                                                                                                                                             <  2:1020099/150‑5 (MQ=255)
              tgAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCATGAACGTTCACTGCAAtc                                                                                                                                                <  2:623403/140‑2 (MQ=255)
              tgAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCATGAACGTTCACTGCAAtc                                                                                                                                                <  2:627967/140‑2 (MQ=255)
              tgAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCATGAACGTTCACTGCAAtc                                                                                                                                                >  1:623403/1‑139 (MQ=255)
              tgAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCATGAACGTTCACTGCAAtc                                                                                                                                                >  1:627967/1‑139 (MQ=255)
               gAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCg                                                                                                                                     >  1:976747/1‑150 (MQ=255)
               gAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCg                                                                                                                                     >  1:971769/1‑150 (MQ=255)
                aaaaGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGt                                                                                                                                    <  1:685063/150‑1 (MQ=255)
                 aaaGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGtc                                                                                                                                   >  2:637409/1‑150 (MQ=255)
                   agagCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGtcat                                                                                                                                 >  2:1327654/1‑150 (MQ=255)
                       cACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCAcc                                                                                                                             <  2:450993/150‑1 (MQ=255)
                          ttttttGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGTGATTGCAGTGAACGTTCTTGCAATAGGCGGATGGCGTCATCTCcaac                                                                                                                          <  2:86757/150‑1 (MQ=255)
                           tttttGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACcaaca                                                                                                                         >  2:908747/1‑150 (MQ=255)
                            ttttGGCCCCGGCGGCAATATGCGCGTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACcaacaa                                                                                                                        <  2:768367/150‑1 (MQ=255)
                                  ccccGGCGGAAATTGGCGCTTCGCCATGCTCGCGGTATCAATATCCGCCGGTGCTGTCGAGCAATACATAAACGCCATTTTCACGCCCGGGTCGCTTTTTCTGTTAAGGTTCATCAAATACGGGTTTGGCTTAATCACAAAAAAAAAgct                                                                                                                  <  2:971769/150‑3 (MQ=255)
                                     cggcggCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGAt                                                                                                               <  1:426289/150‑1 (MQ=255)
                                     cggcggCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGAt                                                                                                               <  2:798109/150‑1 (MQ=255)
                                     cggcggCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGAt                                                                                                               <  2:793143/150‑1 (MQ=255)
                                        cggcAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATcgc                                                                                                            <  2:222821/150‑1 (MQ=255)
                                            aaTATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTc                                                                                                        >  1:1053090/1‑150 (MQ=255)
                                            aaTATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTc                                                                                                        >  1:1424993/1‑150 (MQ=255)
                                             atatGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTcc                                                                                                       <  1:637409/150‑1 (MQ=255)
                                              tatGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCt                                                                                                      >  1:1168601/1‑150 (MQ=255)
                                                tGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGtca                                                                                                                                  >  1:398532/1‑120 (MQ=255)
                                                tGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGtca                                                                                                                                  <  2:398532/120‑1 (MQ=255)
                                                 gcgcTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTc                                                                                                   >  2:1277875/1‑150 (MQ=255)
                                                    cTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTa                                                                                                <  2:1031530/150‑1 (MQ=255)
                                                      tCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACt                                                                                              >  2:408841/1‑150 (MQ=255)
                                                              cTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGc                                                                                      <  1:937071/150‑1 (MQ=255)
                                                                 gcgGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCaaa                                                                                   >  1:1097351/1‑150 (MQ=255)
                                                                    ggAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCATTACGGGGATGGGGTCATCACCAACTAAAAGTTGATCGCGTTCCTGTCGTACTTGCCATGAAAAACg                                                                                <  2:221781/150‑1 (MQ=255)
                                                                            ataCGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGc                                                                        <  2:1053090/150‑1 (MQ=255)
                                                                             taCGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCt                                                                       <  1:816857/150‑1 (MQ=255)
                                                                              aCGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTg                                                                      <  1:784948/150‑1 (MQ=255)
                                                                                gCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGgt                                                                    >  1:331703/1‑150 (MQ=255)
                                                                                   ggTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTc                                                                 >  1:148179/1‑150 (MQ=255)
                                                                                      gCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCata                                                              <  1:1277875/150‑1 (MQ=255)
                                                                                      gCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCata                                                              >  1:524688/1‑150 (MQ=255)
                                                                                         gTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAGCGGCCATGGCTGGTGTCATAttt                                                           >  2:16981/1‑150 (MQ=255)
                                                                                           cGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTca                                                         <  2:148179/150‑1 (MQ=255)
                                                                                           cGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTca                                                         <  1:528682/150‑1 (MQ=255)
                                                                                            gAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTcaa                                                        <  1:1327654/150‑1 (MQ=255)
                                                                                            gAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTcaa                                                        <  1:908747/150‑1 (MQ=255)
                                                                                              ggACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTcaaca                                                      >  1:289185/1‑150 (MQ=255)
                                                                                                aCTACATCAACGCCAAGTTCACGCCAGGGCAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACaaa                                                    >  1:264330/1‑150 (MQ=255)
                                                                                                 cTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAAt                                                   >  1:237277/1‑150 (MQ=255)
                                                                                                 cTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAAt                                                   >  1:606392/1‑150 (MQ=255)
                                                                                                   aCATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATgc                                                 >  2:485861/1‑150 (MQ=255)
                                                                                                     aTCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATgcgc                                               >  1:1356232/1‑150 (MQ=255)
                                                                                                       cAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCa                                             >  1:1292657/1‑150 (MQ=255)
                                                                                                         aCGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATg                                           >  1:1056637/1‑150 (MQ=255)
                                                                                                           gCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGcc                                         >  1:981160/1‑150 (MQ=255)
                                                                                                              aaGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCgc                                      <  2:690055/150‑1 (MQ=255)
                                                                                                               aGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCgca                                     <  2:1168601/150‑1 (MQ=255)
                                                                                                                gTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCgcag                                    >  1:1182161/1‑150 (MQ=255)
                                                                                                                  tCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCgcagca                                  >  1:166377/1‑150 (MQ=255)
                                                                                                                        cAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATTCCCGCAGCATCCGcc                            <  2:524688/150‑1 (MQ=255)
                                                                                                                             gAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTc                       >  1:771468/1‑150 (MQ=255)
                                                                                                                               gCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGt                     >  2:1402854/1‑150 (MQ=255)
                                                                                                                                ggATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCTCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCGCGCAGCATCCGCCAGTTCGtt                    <  2:1248864/149‑1 (MQ=255)
                                                                                                                                cGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGtt                    <  1:869122/150‑1 (MQ=255)
                                                                                                                                 gATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGt                     >  1:1181469/1‑148 (MQ=255)
                                                                                                                                 gATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGt                     <  2:1181469/148‑1 (MQ=255)
                                                                                                                                  aTTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGa                  <  2:981160/150‑1 (MQ=255)
                                                                                                                                  aTTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGa                  <  2:391717/150‑1 (MQ=255)
                                                                                                                                     gCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTg               <  2:1424993/150‑1 (MQ=255)
                                                                                                                                     gCAGTGAACGTTCATGCAATACACGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTg               >  1:1337797/1‑150 (MQ=255)
                                                                                                                                        gTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATGGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGcca            <  2:264330/150‑1 (MQ=255)
                                                                                                                                               ttCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCACCACGGt     <  1:1402854/150‑1 (MQ=255)
                                                                                                                                                 cATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCACCACGGTaa   >  1:598347/1‑150 (MQ=255)
                                                                                                                                                 cATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCACCACGGTaa   >  2:1362433/1‑150 (MQ=255)
                                                                                                                                                  aTGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCACCACGGTAAt  >  1:434413/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCACCACGGTAAT  >  NC_000913/3073309‑3073604

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: