Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,736,157 G→A *393* (TAG→TAA xylR → xylose divergent operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,736,1570GA98.8% 258.2 / NA 80*393* (TAG→TAAxylRxylose divergent operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (44/35);  minor base T (0/1);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAGCGATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAA  >  NC_000913/3736044‑3736287
                                                                                                                 |                                                                                                                                   
atatCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              >  2:1340852/1‑137 (MQ=255)
atatCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              <  1:1340852/137‑1 (MQ=255)
atatCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTca                                                                                                            <  2:1270908/139‑1 (MQ=255)
  atCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          <  1:1607867/139‑1 (MQ=255)
    cGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTc                                                                                                                   >  2:1343312/1‑128 (MQ=255)
    cGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTc                                                                                                                   <  1:1343312/128‑1 (MQ=255)
        aaTGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATTGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGt                                                                                                    >  1:1156570/1‑139 (MQ=255)
         aTGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTa                                                                                                   >  1:697744/1‑139 (MQ=255)
              cGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACga                                                                                                                                 >  2:1427335/1‑103 (MQ=255)
              cGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACga                                                                                                                                 <  1:1427335/103‑1 (MQ=255)
                gTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGc                                                                                            >  1:670512/1‑139 (MQ=255)
                 ttATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           <  2:697744/139‑1 (MQ=255)
                 ttATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           <  1:1124595/139‑1 (MQ=255)
                     ccATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGt                                                                                       >  1:1307946/1‑139 (MQ=255)
                      cATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGtt                                                                                      <  2:484301/139‑1 (MQ=255)
                        tCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGt                                                                                    >  2:496880/1‑139 (MQ=255)
                         cGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              >  1:1228238/1‑112 (MQ=255)
                         cGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              <  2:1228238/112‑1 (MQ=255)
                          gCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCg                                                                                  >  1:776725/1‑139 (MQ=255)
                             gCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGAt                                                                               <  2:776725/139‑1 (MQ=255)
                                   tttCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATgcggcg                                                                         >  1:998917/1‑139 (MQ=255)
                                      cTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGa                                                                      >  2:559281/1‑139 (MQ=255)
                                        aCTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAAc                                                                    >  2:482810/1‑139 (MQ=255)
                                         ctctGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGTTGCGTGGTTGTCGGATGCGGCGTGAACg                                                                   <  2:998917/139‑1 (MQ=255)
                                         ctctGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTg                                                                                             >  2:923967/1‑113 (MQ=255)
                                         ctctGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTg                                                                                             <  1:923967/113‑1 (MQ=255)
                                             gTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCtt                                                               >  1:1081706/1‑139 (MQ=255)
                                              tttttAAAAAAGCATATGACGCGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTCAACGCCTTa                                                              >  1:948759/1‑139 (MQ=255)
                                              tttttAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTa                                                              >  2:850204/1‑139 (MQ=255)
                                                        aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          >  1:351537/1‑85 (MQ=255)
                                                        aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          <  2:351537/85‑1 (MQ=255)
                                                        aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTAc                                                    >  1:668981/1‑139 (MQ=255)
                                                        aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTAc                                                    >  1:1162530/1‑139 (MQ=255)
                                                             atGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGcc                                               <  1:385890/139‑1 (MQ=255)
                                                                acacGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAgg                                            <  2:1214277/139‑1 (MQ=255)
                                                                  acgacgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCa                                          >  2:250547/1‑139 (MQ=255)
                                                                  acgacgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCa                                          <  2:668981/139‑1 (MQ=255)
                                                                      cgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           <  1:1574186/86‑1 (MQ=255)
                                                                      cgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           >  2:1574186/1‑86 (MQ=255)
                                                                      cgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGAt                                      >  1:117839/1‑139 (MQ=255)
                                                                         cAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATaaa                                   >  1:23842/1‑139 (MQ=255)
                                                                          aaaaGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATaaaa                                  <  2:1307946/139‑1 (MQ=255)
                                                                            aaGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACg                                >  2:895930/1‑139 (MQ=255)
                                                                                gTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCAcc                            <  1:250547/139‑1 (MQ=255)
                                                                                 tATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAGTTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACca                           >  2:1530401/1‑139 (MQ=255)
                                                                                 tATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACca                           >  2:429298/1‑139 (MQ=255)
                                                                                 tATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACca                           >  2:701825/1‑139 (MQ=255)
                                                                                  aTCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACcaa                          >  2:1577477/1‑139 (MQ=255)
                                                                                    cgcgATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGTATGATAAAACGCACcaaca                        <  1:496880/139‑1 (MQ=255)
                                                                                     gcgATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACcaacaa                       <  2:117839/139‑1 (MQ=255)
                                                                                       gATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGa                                                                      >  1:4463/1‑90 (MQ=255)
                                                                                       gATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGa                                                                      <  2:4463/90‑1 (MQ=255)
                                                                                            aaaTAGCGAGGTCATGTTGTATCGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCa                <  2:948759/139‑1 (MQ=255)
                                                                                            aaaTAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCa                <  1:1577477/139‑1 (MQ=255)
                                                                                                aGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCt            <  2:1162530/139‑1 (MQ=255)
                                                                                                  cGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGcc                                                                 >  2:879661/1‑84 (MQ=39)
                                                                                                  cGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGcc                                                                 <  1:879661/84‑1 (MQ=39)
                                                                                                   gAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          <  2:483632/42‑1 (MQ=255)
                                                                                                   gAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          >  1:483632/1‑42 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               >  1:1080443/1‑35 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               <  2:1080443/35‑1 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              <  2:1256864/36‑1 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              >  1:1256864/1‑36 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           >  2:1216956/1‑55 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           <  1:594293/55‑1 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           >  2:594293/1‑55 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           <  1:1216956/55‑1 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGcc                                               >  2:680226/1‑99 (MQ=38)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGcc                                               <  1:680226/99‑1 (MQ=38)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAgg                                            >  2:1267946/1‑102 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAgg                                            <  1:1267946/102‑1 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCt            >  1:1170829/1‑134 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCt            <  2:1170829/134‑1 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAg       >  2:239604/1‑139 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAg       >  2:757820/1‑139 (MQ=255)
                                                                                                     ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAg       >  1:889660/1‑139 (MQ=255)
                                                                                                       tCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCt     >  2:791318/1‑139 (MQ=255)
                                                                                                          tgttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCaa  <  2:1532820/139‑1 (MQ=255)
                                                                                                          tgttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCaa  <  1:482810/139‑1 (MQ=255)
                                                                                                          tgttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCaa  >  1:1532820/1‑139 (MQ=255)
                                                                                                                 |                                                                                                                                   
ATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAGCGATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAA  >  NC_000913/3736044‑3736287

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: