Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,450,962 T→G 5.4% intergenic (+44/‑96) ytfQ → / → ytfR galactofuranose ABC transporter periplasmic binding protein/putative sugar ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,450,9620TG 5.4% 595.3 / 21.7 166intergenic (+44/‑96)ytfQ/ytfRgalactofuranose ABC transporter periplasmic binding protein/putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (80/77);  new base G (4/5);  total (84/82)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

Only 100 of 168 total aligned reads displayed.

AAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATT  >  NC_000913/4450814‑4451110
                                                                                                                                                    |                                                                                                                                                    
    aCGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgc                                                                                                                                                 <  4:1113123/149‑1 (MQ=255)
     cGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgc                                                                                                                                                 <  2:1095522/149‑1 (MQ=255)
      ggCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                                                                                                                                <  3:814622/149‑1 (MQ=255)
      ggCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                                                                                                                                >  1:375661/1‑149 (MQ=255)
         aCCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCa                                                                                                                                             >  2:97339/1‑149 (MQ=255)
         aCCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCa                                                                                                                                             >  1:119918/1‑149 (MQ=255)
             tGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCg                                                                                                                                         >  2:356131/1‑149 (MQ=255)
              gCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCAGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCgg                                                                                                                                        >  4:115717/1‑149 (MQ=255)
                cTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCg                                                                                                                                      >  3:253748/1‑149 (MQ=255)
                cTGAAAAGCTGACGTTAACCAAAGCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                                      <  3:656572/149‑2 (MQ=255)
                 tGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                                                                                                     >  3:922215/1‑149 (MQ=255)
                  gAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGAt                                                                                                                                    <  1:1125309/149‑1 (MQ=255)
                     aaGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGt                                                                                                                                 <  3:1343820/149‑1 (MQ=255)
                      aGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTg                                                                                                                                <  4:922215/149‑1 (MQ=255)
                          gACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACt                                                                                                                            <  3:472461/149‑1 (MQ=255)
                             gTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                                                                                                         >  2:376273/1‑149 (MQ=255)
                             gTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                                                                                                         <  3:1066291/149‑1 (MQ=255)
                             gTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                                                                                                         <  3:41838/149‑1 (MQ=255)
                             gTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                                                                                                         >  1:1380083/1‑149 (MQ=255)
                               tAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGc                                                                                                                       >  2:728733/1‑149 (MQ=255)
                                 aCCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGTTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                                                                                                     >  4:840055/1‑149 (MQ=255)
                                 aCCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                                                                                                     >  1:1480452/1‑149 (MQ=255)
                                  ccAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                                                                                    >  4:323205/1‑149 (MQ=255)
                                   cAAATCCACCCTATACCTGCCTGATACCGCAAAAGATGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCGTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGAGGGTGCACTGAAGCCTGa                                                                                                                   >  4:1232862/1‑149 (MQ=255)
                                    aaaTCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                                                                                  <  1:8939/149‑1 (MQ=255)
                                    aaaTCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                                                                                  <  3:323205/149‑1 (MQ=255)
                                      aTCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                                       <  2:338785/126‑1 (MQ=255)
                                          aCCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGcgccttat                                                                                                                                            >  2:1267667/1‑113 (MQ=255)
                                          aCCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACg                                                                                                            >  4:1449783/1‑149 (MQ=255)
                                          aCCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACg                                                                                                            >  3:191082/1‑149 (MQ=255)
                                          aCCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCGAGcgcctt                                                                                                                                              <  2:1277385/115‑3 (MQ=255)
                                          aCCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCGAGcgcctt                                                                                                                                              >  1:1277385/1‑113 (MQ=255)
                                             cTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                                                                                  <  1:188088/140‑1 (MQ=255)
                                             cTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                                                                                  >  2:188088/1‑140 (MQ=255)
                                                tACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                                                                                      <  1:97798/149‑1 (MQ=255)
                                                  cctgcctgATACCGCAAAAGAACAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGTATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                                                                                                    <  3:893379/149‑1 (MQ=255)
                                                    tgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCa                                                                                                                              <  2:1383962/121‑1 (MQ=255)
                                                    tgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                                                                                  >  1:99517/1‑149 (MQ=255)
                                                     gcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                                                                                 <  2:738421/149‑1 (MQ=255)
                                                       ctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                                                                                               >  1:728856/1‑149 (MQ=255)
                                                        tgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                                                                              <  2:794271/149‑1 (MQ=255)
                                                        tgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                                                                              >  4:1393039/1‑149 (MQ=255)
                                                         gATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                                                                                             <  1:639694/149‑1 (MQ=255)
                                                          aTACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                                                                                            <  2:89379/149‑1 (MQ=255)
                                                          aTACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                                                                                            <  2:1195039/149‑1 (MQ=255)
                                                          aTACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                                                                                            <  2:822099/149‑1 (MQ=255)
                                                           tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCGGGCGATTGTGCACCGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                                                                                           >  4:980424/1‑149 (MQ=255)
                                                              cGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                                                                                        >  2:250601/1‑149 (MQ=255)
                                                                 aaaaGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                                                                                     >  1:1238542/1‑149 (MQ=255)
                                                                      aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                                                                                <  1:1062216/149‑1 (MQ=255)
                                                                       agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                                                                               <  2:1238542/149‑1 (MQ=255)
                                                                        gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                                                                              >  1:1376057/1‑149 (MQ=255)
                                                                         aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTACTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc                                                                             >  1:180873/1‑149 (MQ=255)
                                                                          agtAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa                                                                            >  1:683615/3‑149 (MQ=255)
                                                                          aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa                                                                            >  4:1205417/1‑149 (MQ=255)
                                                                           ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                                                                              >  2:986742/1‑130 (MQ=255)
                                                                           ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc                                                                           <  3:302450/149‑1 (MQ=255)
                                                                             agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGc                                                                         >  4:297023/1‑149 (MQ=255)
                                                                              gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt                                                                        <  3:1089155/149‑1 (MQ=255)
                                                                              gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt                                                                        <  1:409501/149‑1 (MQ=255)
                                                                              gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt                                                                        >  1:332044/1‑149 (MQ=255)
                                                                               agaagaagaAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCtt                                                                       >  4:617122/1‑149 (MQ=255)
                                                                               agaagaagaAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCtt                                                                       <  4:419390/149‑1 (MQ=255)
                                                                                 aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTAGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTa                                                                     >  1:957436/1‑149 (MQ=255)
                                                                                      gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTATGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg                                                                >  3:686950/1‑149 (MQ=255)
                                                                                      gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg                                                                >  3:1041800/1‑149 (MQ=255)
                                                                                      gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg                                                                >  4:812152/1‑149 (MQ=255)
                                                                                       aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCgg                                                               <  3:812152/149‑1 (MQ=255)
                                                                                        aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCggg                                                              <  1:970063/149‑1 (MQ=255)
                                                                                        aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCggg                                                              <  3:1146044/149‑1 (MQ=255)
                                                                                        aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCggg                                                              <  3:663491/149‑1 (MQ=255)
                                                                                           atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGa                                                           >  1:833528/1‑149 (MQ=255)
                                                                                             atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGaaa                                                         <  3:402231/149‑1 (MQ=255)
                                                                                              tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAAc                                                        <  4:188416/149‑1 (MQ=255)
                                                                                               ggggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAAcc                                                       >  4:429946/1‑149 (MQ=255)
                                                                                                gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCa                                                      >  1:520397/1‑149 (MQ=255)
                                                                                                gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCa                                                      >  1:354484/1‑149 (MQ=255)
                                                                                                   tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                                       <  2:624573/65‑1 (MQ=17)
                                                                                                   tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                                       >  1:624573/1‑65 (MQ=18)
                                                                                                   tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATga                                                   <  3:121613/149‑1 (MQ=255)
                                                                                                     ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                                                                                                         >  1:993061/1‑77 (MQ=255)
                                                                                                      tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATgacga                                                <  3:1316079/149‑1 (MQ=255)
                                                                                                        aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGAcgacc                                              <  4:157563/149‑1 (MQ=255)
                                                                                                          ggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGAcgaccga                                            <  4:1286688/149‑1 (MQ=255)
                                                                                                           gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGAcgaccgac                                           <  2:99517/149‑1 (MQ=255)
                                                                                                            ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACgaccgacc                                          <  3:511081/149‑1 (MQ=255)
                                                                                                              gCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCaa                                        <  1:415957/149‑1 (MQ=255)
                                                                                                              gCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCaa                                        >  4:48290/1‑149 (MQ=255)
                                                                                                               cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAac                                       <  4:1409705/149‑1 (MQ=255)
                                                                                                                tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAaca                                      >  4:249381/1‑149 (MQ=255)
                                                                                                                    ccTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAg                                  <  4:714640/149‑1 (MQ=255)
                                                                                                                      tGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGga                                <  4:1075188/149‑1 (MQ=255)
                                                                                                                        aTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGgaga                              >  1:758302/1‑149 (MQ=255)
                                                                                                                          gCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGAtc                            >  3:1403290/1‑149 (MQ=255)
                                                                                                                           ccGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATgacgac                                               <  3:58034/129‑1 (MQ=255)
                                                                                                                           ccGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATgacgac                                               >  4:58034/1‑129 (MQ=255)
                                                                                                                                    aGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCg                  <  2:119918/149‑1 (MQ=255)
                                                                                                                                    aGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCg                  >  1:23533/1‑149 (MQ=255)
                                                                                                                                         ggATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGa             >  4:168467/1‑149 (MQ=255)
                                                                                                                                              aGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAg        <  1:1260044/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAATT  >  NC_000913/4450814‑4451110

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: