Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,232,778 A→G 5.7% F152F (TTT→TTC) ‡ lysC ← lysine‑sensitive aspartokinase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,232,7780AG 5.7% 368.8 / 9.0 123F152F (TTT→TTC) ‡lysClysine‑sensitive aspartokinase 3
Reads supporting (aligned to +/- strand):  ref base A (43/73);  new base G (5/2);  total (48/75)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

Only 100 of 125 total aligned reads displayed.

CGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGC  >  NC_000913/4232630‑4232927
                                                                                                                                                    |                                                                                                                                                     
  cTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAAt                                                                                                                                                    <  1:458953/150‑1 (MQ=255)
    gCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCg                                                                                                                                                  <  1:1004457/150‑1 (MQ=255)
    gCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCg                                                                                                                                                  >  1:990065/1‑150 (MQ=255)
     ccTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGCCCAAATCgg                                                                                                                                                 >  1:209340/1‑150 (MQ=255)
      cTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCggg                                                                                                                                                >  2:1691467/1‑149 (MQ=255)
      cTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGt                                                                                                                                                <  2:77874/150‑1 (MQ=255)
        ccACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCg                                                                                                                                              <  2:1097583/150‑1 (MQ=255)
        ccACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTACCGCCAGCGCGGCTATAACTGGCTCTGCACGACCAAATCGCTCg                                                                                                                                              >  2:888966/1‑150 (MQ=255)
            ggCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCGATTTGGTCGtgca                                                                                                                                          <  2:745904/150‑4 (MQ=255)
            ggCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTgg                                                                                                                                          <  2:837893/150‑1 (MQ=255)
            ggCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTgg                                                                                                                                          >  1:1523467/1‑150 (MQ=255)
            ggCCAAGCGTCGTTGTACGACATTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTgg                                                                                                                                          <  2:1556736/150‑1 (MQ=255)
             gCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGt                                                                                                                                         <  1:953620/150‑1 (MQ=255)
              ccAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCGATTTGGTCGtgcaga                                                                                                                                        <  1:332040/150‑6 (MQ=255)
               cAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTAc                                                                                                                                       >  2:57087/1‑150 (MQ=255)
                  gCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCa                                                                                                                                    <  1:1215592/150‑1 (MQ=255)
                      cgtTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCAc                                                                                                                                <  2:50974/150‑1 (MQ=255)
                      cgtTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCAc                                                                                                                                <  1:570444/150‑1 (MQ=255)
                      cgtTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCAc                                                                                                                                >  1:1511905/1‑150 (MQ=255)
                          gTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTa                                                                                                                            >  2:275228/1‑150 (MQ=255)
                          gTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGGGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTa                                                                                                                            <  1:180645/150‑1 (MQ=255)
                           gaCGACCTTTATTTTCGCTACCGATAAATCCCTGGGAGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTGCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTAc                                                                                                                           <  1:1507161/149‑1 (MQ=255)
                            aCGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACg                                                                                                                          <  1:707203/150‑1 (MQ=255)
                             cGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGt                                                                                                                         >  1:557480/1‑150 (MQ=255)
                               aCCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTAc                                                                                                                       <  2:438289/150‑1 (MQ=255)
                                ccTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACa                                                                                                                      >  1:384964/1‑150 (MQ=255)
                                 cTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACAt                                                                                                                     >  1:790449/1‑150 (MQ=255)
                                  tttATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATc                                                                                                                    <  1:1606967/150‑1 (MQ=255)
                                   ttATTTTCGCTACCGATAAATCCTTGGGTGATCACTATGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCa                                                                                                                   <  1:623976/150‑1 (MQ=255)
                                   ttATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCa                                                                                                                   <  1:919889/150‑1 (MQ=255)
                                        ttCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCa                                                                                                              <  2:1438300/150‑1 (MQ=255)
                                         tCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCAc                                                                                                             <  2:376531/150‑1 (MQ=255)
                                          cGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACt                                                                                                            >  2:153084/1‑150 (MQ=255)
                                             tACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACtgtg                                                                                                         <  2:1523467/150‑1 (MQ=255)
                                             tACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACtgtg                                                                                                         <  1:440554/150‑1 (MQ=255)
                                             tACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACtgtg                                                                                                         <  2:111436/150‑1 (MQ=255)
                                               ccGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGcc                                                                                                       <  2:909681/150‑1 (MQ=255)
                                                    ataTCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGaat                                                                                                  <  1:246633/148‑2 (MQ=255)
                                                      aTCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGCCCAAATAGCTCGTTGGTACGCATCACTTTACGAAAATCAAACCACTGTGCCTGAACAt                                                                                                >  1:1642352/1‑150 (MQ=255)
                                                        cccTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATcg                                                                                              <  1:1131628/150‑1 (MQ=255)
                                                          cTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATcgcg                                                                                            >  2:32275/1‑150 (MQ=255)
                                                            gggTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGCCCAAATCGGTCGTTTGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGtt                                                                                          >  2:1122598/1‑150 (MQ=255)
                                                            gggTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATcg                                                                                              <  1:501468/146‑1 (MQ=255)
                                                            gggTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATcg                                                                                              >  2:501468/1‑146 (MQ=255)
                                                             ggTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTc                                                                                         >  1:705783/1‑150 (MQ=255)
                                                              gTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTcg                                                                                        <  2:936966/150‑1 (MQ=255)
                                                                   cACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGATCAACAACCCACTGTGCCTGAACCTCGCGTTCgcgctg                                                                                   >  2:319253/1‑148 (MQ=255)
                                                                    aCTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGCCCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAgg                                                                                  >  2:1046201/1‑150 (MQ=255)
                                                                    aCTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAgg                                                                                  >  1:1231130/1‑150 (MQ=255)
                                                                    aCTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAgg                                                                                  <  1:27042/150‑1 (MQ=255)
                                                                       aaGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGAtc                                                                               <  2:384964/150‑1 (MQ=255)
                                                                       aaGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGGTCCGCCAGCGCGGCTATAGCTGGCTCTGCACGACCAAATCGGTCGGTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGAtc                                                                               <  1:956865/150‑1 (MQ=255)
                                                                         gCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGAtctc                                                                             <  2:268268/150‑1 (MQ=255)
                                                                          ccTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTca                                                                            <  1:1026935/150‑1 (MQ=255)
                                                                          ccTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTca                                                                            >  2:104579/1‑150 (MQ=255)
                                                                          ccTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTca                                                                            >  1:520032/1‑150 (MQ=255)
                                                                           cTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCATATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTcaa                                                                           <  2:443526/150‑1 (MQ=255)
                                                                            ttCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTcaac                                                                          <  2:595324/150‑1 (MQ=255)
                                                                             tCATTGAGACGTGGGAGCAGCGGCAGCGCGGCCAGGTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGGCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTcaaca                                                                         <  1:1610794/150‑1 (MQ=255)
                                                                              cATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCGAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTcaacaa                                                                        <  2:1568024/150‑1 (MQ=255)
                                                                              cATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCCAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTcaacaa                                                                        >  1:1569233/1‑150 (MQ=255)
                                                                               aTTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaa                                                                       >  2:1682304/1‑150 (MQ=255)
                                                                                ttGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaac                                                                      <  1:895286/150‑1 (MQ=255)
                                                                                ttGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaac                                                                      <  1:1237107/150‑1 (MQ=255)
                                                                                ttGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaac                                                                      <  1:1713837/150‑1 (MQ=255)
                                                                                 tGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCCCTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaaca                                                                     >  2:1114796/1‑150 (MQ=255)
                                                                                 tGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaaca                                                                     >  2:1426524/1‑150 (MQ=255)
                                                                                  gagaCGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAAcag                                                                    <  1:737982/150‑1 (MQ=255)
                                                                                  gagaCGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAAcag                                                                    <  1:1372691/150‑1 (MQ=255)
                                                                                      cGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCCTCACTGTTCGTACATCAACCGACTGTGCCTGAACATCGCGTTCGCGGAGGATCTCAACAAACAGCAgg                                                                >  1:263738/1‑150 (MQ=255)
                                                                                      cGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAgg                                                                >  2:1718057/1‑150 (MQ=255)
                                                                                          ggAGCAGCTGCTGCGCGGCCAGGTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCg                                                            <  2:285175/150‑1 (MQ=255)
                                                                                          ggAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCg                                                            <  1:1328791/150‑1 (MQ=255)
                                                                                             gcagcTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACa                                                         <  1:574494/150‑1 (MQ=255)
                                                                                               agcTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATc                                                       <  2:594574/150‑1 (MQ=255)
                                                                                                gcTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCa                                                      <  2:790449/150‑1 (MQ=255)
                                                                                                  tGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGc                                                    <  2:887572/150‑1 (MQ=255)
                                                                                                  tGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGc                                                    <  2:624182/150‑1 (MQ=255)
                                                                                                    cAGCGCGGCCAGGTCCGCCAGCGCGGCTATATCTGGCGCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTc                                                  <  1:1219918/150‑1 (MQ=255)
                                                                                                         cgGCCAGTTCCGACAGCGCGGCTATATCTGGCTCTGCACGACAAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGTATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGt                                             >  1:134279/1‑150 (MQ=255)
                                                                                                            ccAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGc                                          <  1:1680401/150‑1 (MQ=255)
                                                                                                               gTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGa                                       >  1:1030027/1‑150 (MQ=255)
                                                                                                                 tCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGAcc                                     <  1:1431549/150‑1 (MQ=255)
                                                                                                                   cGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAg                                   <  2:11541/150‑1 (MQ=255)
                                                                                                                               tatcTCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAAGATCGCGTTCGCGCAGGATCTCACCAAACACCTGGGTCGACATCAGCTCGCAGTGGCTGACGAGATCATCAGTccg                       >  1:254625/1‑148 (MQ=255)
                                                                                                                                  atCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAgcgc                    >  1:313485/1‑150 (MQ=255)
                                                                                                                                   tCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGcc                   >  2:642949/1‑150 (MQ=255)
                                                                                                                                    cTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCg                  >  2:1247880/1‑150 (MQ=255)
                                                                                                                                       gCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGgag               <  1:108836/150‑1 (MQ=255)
                                                                                                                                       gCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGgag               <  2:1028942/150‑1 (MQ=255)
                                                                                                                                       gCTCTGCACGACCAAATCGGTAGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGgag               <  2:1729531/150‑1 (MQ=255)
                                                                                                                                        ctctGCACGACCAAATCGGTCGTTGGTACGCAGCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGgaga              <  1:1109046/150‑1 (MQ=255)
                                                                                                                                          accaaCGACCGATTTGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACg            >  2:899172/5‑150 (MQ=255)
                                                                                                                                           ccaaCGACCGATTTGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCCGCTCATCTGTCAGCGCCGGAGACGt           >  1:1245722/4‑150 (MQ=255)
                                                                                                                                           ccaaCGACCGATTTGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGt           >  2:609254/4‑150 (MQ=255)
                                                                                                                                            gCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGtt          >  1:1411044/1‑150 (MQ=255)
                                                                                                                                             cACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCACGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTg         <  2:359262/150‑1 (MQ=255)
                                                                                                                                               cGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATAAAACCACTGTTCCTGAACATAGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGcc       <  2:98428/150‑1 (MQ=255)
                                                                                                                                                  ccAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAgc    >  1:406912/1‑150 (MQ=255)
                                                                                                                                                    aaaTCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAgcgc  >  1:1123576/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCTGCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGC  >  NC_000913/4232630‑4232927

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