Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,232,782 C→T 5.7% R151Q (CGA→CAA)  lysC ← lysine‑sensitive aspartokinase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,232,7820CT 5.7% 393.3 / 7.4 123R151Q (CGA→CAA) lysClysine‑sensitive aspartokinase 3
Reads supporting (aligned to +/- strand):  ref base C (43/72);  new base T (5/2);  total (49/74)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 127 total aligned reads displayed.

GCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCC  >  NC_000913/4232634‑4232931
                                                                                                                                                    |                                                                                                                                                     
gCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCg                                                                                                                                                      <  1:1004457/150‑1 (MQ=255)
 ccTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGCCCAAATCgg                                                                                                                                                     >  1:209340/1‑150 (MQ=255)
  cTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCggg                                                                                                                                                    >  2:1691467/1‑149 (MQ=255)
  cTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGt                                                                                                                                                    <  2:77874/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)
        ggCCAAGCGTCGTTGTACGACATTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTgg                                                                                                                                              <  2:1556736/150‑1 (MQ=255)
         gCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGt                                                                                                                                             <  1:953620/150‑1 (MQ=255)
          ccAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCGATTTGGTCGtgcaga                                                                                                                                            <  1:332040/150‑6 (MQ=255)
           cAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTAc                                                                                                                                           <  2:1144189/150‑1 (MQ=255)
           cAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTAc                                                                                                                                           >  2:57087/1‑150 (MQ=255)
              gCGTCGTTGTACGACCGTTATTTTCGCTACCGATAAATCCCTGGGTGAGCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCa                                                                                                                                        <  1:1215593/150‑1 (MQ=255)
                  cgtTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCAc                                                                                                                                    >  1:1511905/1‑150 (MQ=255)
                  cgtTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCAc                                                                                                                                    <  2:50974/150‑1 (MQ=255)
                   gtTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACt                                                                                                                                   <  2:779961/150‑1 (MQ=255)
                      gTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTa                                                                                                                                >  2:275228/1‑150 (MQ=255)
                       tACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTAc                                                                                                                               <  2:236108/150‑1 (MQ=255)
                       gaCGACCTTTATTTTCGCTACCGATAAATCCCTGGGAGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTGCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTAc                                                                                                                               <  1:1507161/149‑1 (MQ=255)
                         cGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGt                                                                                                                             >  1:557480/1‑150 (MQ=255)
                         cGACCGTTATTTTCGCCCCCGCTAAATCCCGCCGTGAGCCCTAAGCCTTCATTGAGACGTGGGAGCAGCTGCATCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAACTCGGGCGTTGGTACGCATCACTTTCCGt                                                                                                                             <  2:1065182/150‑1 (MQ=255)
                           aCCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTAc                                                                                                                           <  2:438289/150‑1 (MQ=255)
                             cTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACAt                                                                                                                         >  2:895286/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)
                                    ttCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCa                                                                                                                  <  2:1438300/150‑1 (MQ=255)
                                     tCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCAc                                                                                                                 <  2:376531/150‑1 (MQ=255)
                                      cGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACt                                                                                                                >  2:153084/1‑150 (MQ=255)
                                       gCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACtg                                                                                                               >  1:1107979/1‑150 (MQ=255)
                                         tACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACtgtg                                                                                                             <  2:1523467/150‑1 (MQ=255)
                                         tACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACtgtg                                                                                                             <  2:111436/150‑1 (MQ=255)
                                         tACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACtgtg                                                                                                             <  1:440554/150‑1 (MQ=255)
                                          aCCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGc                                                                                                            >  1:888846/1‑150 (MQ=255)
                                           ccGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGcc                                                                                                           <  2:909681/150‑1 (MQ=255)
                                             gATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTg                                                                                                         <  1:153084/150‑1 (MQ=255)
                                                ataTCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGaat                                                                                                      <  1:246633/148‑2 (MQ=255)
                                                    cccTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATcg                                                                                                  <  1:1131628/150‑1 (MQ=255)
                                                      cTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATcgcg                                                                                                >  2:32275/1‑150 (MQ=255)
                                                       tGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGt                                                                                               <  1:1423836/150‑1 (MQ=255)
                                                        gggTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGCCCAAATCGGTCGTTTGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGtt                                                                                              >  2:1122598/1‑150 (MQ=255)
                                                        gggTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATcg                                                                                                  >  2:501468/1‑146 (MQ=255)
                                                         ggTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTc                                                                                             >  1:705783/1‑150 (MQ=255)
                                                          gTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTcg                                                                                            <  2:936966/150‑1 (MQ=255)
                                                          gTGATCACTAAGCCTTCATAGCGACGTGGGAGCAGCTGCAGCGCGGCCAGTGCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTcg                                                                                            <  1:98588/150‑1 (MQ=255)
                                                               cACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGATCAACAACCCACTGTGCCTGAACCTCGCGTTCgcgctg                                                                                       >  2:319253/1‑148 (MQ=255)
                                                                aCTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGCCCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAgg                                                                                      >  2:1046201/1‑150 (MQ=255)
                                                                aCTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAgg                                                                                      <  1:27042/150‑1 (MQ=255)
                                                                aCTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAgg                                                                                      >  1:1231130/1‑150 (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)
                                                                       cTTCATTGAGACTTGGGAGCAGATGCAGCGCGGCCAGATCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTcaa                                                                               <  2:321246/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)
                                                                           aTTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaa                                                                           >  2:1638223/1‑150 (MQ=255)
                                                                            ttGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaac                                                                          <  1:1237107/150‑1 (MQ=255)
                                                                            ttGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaac                                                                          <  1:1713837/150‑1 (MQ=255)
                                                                            ttGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaac                                                                          <  1:895286/150‑1 (MQ=255)
                                                                             tGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCCCTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaaca                                                                         >  2:1114796/1‑150 (MQ=255)
                                                                             tGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCaacaaaca                                                                         >  2:1426524/1‑150 (MQ=255)
                                                                              gagaCGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAAcag                                                                        <  1:1372691/150‑1 (MQ=255)
                                                                              gagaCGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAAcag                                                                        <  1:737982/150‑1 (MQ=255)
                                                                                  cGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAgg                                                                    >  2:1718057/1‑150 (MQ=255)
                                                                                     gggAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTc                                                                 <  1:670069/150‑1 (MQ=255)
                                                                                      ggAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCg                                                                <  1:1328791/150‑1 (MQ=255)
                                                                                       gAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCgc                                                               >  1:1662115/1‑149 (MQ=255)
                                                                                         gcagcTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACa                                                             <  1:574494/150‑1 (MQ=255)
                                                                                            gcTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCa                                                          >  1:437380/1‑150 (MQ=255)
                                                                                              tGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGc                                                        <  2:887572/150‑1 (MQ=255)
                                                                                                cAGCGCGGCCAGGTCCGCCAGCGCGGCTATATCTGGCGCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTc                                                      <  1:1219918/150‑1 (MQ=255)
                                                                                                 aGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCg                                                     <  2:620108/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)
                                                                                                                           tatcTCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAAGATCGCGTTCGCGCAGGATCTCACCAAACACCTGGGTCGACATCAGCTCGCAGTGGCTGACGAGATCATCAGTccg                           >  1:254625/1‑148 (MQ=255)
                                                                                                                             tatCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAgcg                         >  1:53745/1‑150 (MQ=255)
                                                                                                                              atCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAgcgc                        >  1:313485/1‑150 (MQ=255)
                                                                                                                               tCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGcc                       >  2:642949/1‑150 (MQ=255)
                                                                                                                                cTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCg                      >  2:1247880/1‑150 (MQ=255)
                                                                                                                                   gCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGgag                   <  2:1028942/150‑1 (MQ=255)
                                                                                                                                    ctctGCACGACCAAATCGGTCGTTGGTACGCAGCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGgaga                  <  1:1109046/150‑1 (MQ=255)
                                                                                                                                      ctGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACg                <  2:557480/150‑1 (MQ=255)
                                                                                                                                       ccaaCGACCGATTTGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACCGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGCCGt               >  1:222713/4‑150 (MQ=255)
                                                                                                                                       ccaaCGACCGATTTGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGt               >  2:609254/4‑150 (MQ=255)
                                                                                                                                         aaCGACCGATTTGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTg             >  2:433614/2‑150 (MQ=255)
                                                                                                                                           cGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGcc           <  1:275228/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)
                                                                                                                                                 aaTCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCAGAGACGTTGCCAGcgcc     <  1:1027999/150‑1 (MQ=255)
                                                                                                                                                    cGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGcgccgcc  >  2:203668/1‑150 (MQ=255)
                                                                                                                                                    cGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGcgccgcc  >  1:464290/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCTCCACGGCCAAGCGTCGTTGTACGACCTTTATTTTCGCTACCGATAAATCCCTGGGTGATCACTAAGCCTTCATTGAGACGTGGGAGCAGCTGCAGCGCGGCCAGTTCCGCCAGCGCGGCTATATCTGGCTCTGCACGACCAAATCGGTCGTTGGTACGCATCACTTTACGTACATCAAACCACTGTGCCTGAACATCGCGTTCGCGCAGGATCTCAACAAACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCC  >  NC_000913/4232634‑4232931

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: ‑