Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,073,445 Δ1 bp 6.9% coding (247/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,4450G. 6.9% 239.1 / 11.1 72coding (247/1020 nt)epdD‑erythrose‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base . (2/3);  total (34/38)
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.03e-01

GTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCA  >  NC_000913/3073299‑3073594
                                                                                                                                                  |                                                                                                                                                      
gTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGa                                                                                                                                                     >  2:869122/1‑150 (MQ=255)
   gagaTCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACg                                                                                                                                                  <  1:172613/150‑1 (MQ=255)
          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                                                                                                                                      >  1:623403/1‑139 (MQ=255)
                        tgAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCATGAACGTTCACTGCAAtc                                                                                                                                      >  1:627967/1‑139 (MQ=255)
                        tgAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCATGAACGTTCACTGCAAtc                                                                                                                                      <  2:627967/140‑2 (MQ=255)
                        tgAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCATGAACGTTCACTGCAAtc                                                                                                                                      <  2:623403/140‑2 (MQ=255)
                         gAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCg                                                                                                                           >  1:971769/1‑150 (MQ=255)
                         gAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCg                                                                                                                           >  1:976747/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                                                                                                     <  2:793143/150‑1 (MQ=255)
                                               cggcggCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGAt                                                                                                     <  2:798109/150‑1 (MQ=255)
                                               cggcggCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGAt                                                                                                     <  1:426289/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                                                                                                                        <  2:398532/120‑1 (MQ=255)
                                                          tGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGtca                                                                                                                        >  1:398532/1‑120 (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)
                                                                                                                                                  |                                                                                                                                                      
GTCGAGATCGTTACTGCCAGGATGTGAAAAGAGCACTTTTTTGGCCCCGGCGGCAATATGCGCTTCGCCATGCTCGCGGGAGCCATATACGCCGGTGCAGTCGAGGACTACATCAACGCCAAGTTCACGCCAGGGGAGCGATTGCAGTGAACGTTCATGCAATACGCGGATGGCGTCATCACCAACAAAAAGTTGATCGCGTTCCTGTCGTACTTCCCATGCAAAACGGCCATGGCTGGTGTCATATTTCAACAAATGCGCCATGCCCGCAGCATCCGCCAGTTCGTTGATTGCCA  >  NC_000913/3073299‑3073594

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: