Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,130,616 C→A 8.4% A261E (GCG→GAG)  metL → Bifunctional aspartokinase/homoserine dehydrogenase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,130,6160CA 8.4% 1891.7 / 116.6 554A261E (GCG→GAG) metLBifunctional aspartokinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (265/242);  new base A (22/25);  total (287/267)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 589 total aligned reads displayed.

TATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATT  >  NC_000913/4130468‑4130764
                                                                                                                                                    |                                                                                                                                                    
tATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGgc                                                                                                                                                      >  1:5596299/1‑149 (MQ=255)
tATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGgc                                                                                                                                                      >  3:3729428/1‑149 (MQ=255)
  ttCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGgcgc                                                                                                                                                    >  1:6482676/1‑149 (MQ=255)
     cgcgACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGcc                                                                                                                                                 >  3:6477611/1‑149 (MQ=255)
     cgcgACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGcc                                                                                                                                                 <  4:1966140/149‑1 (MQ=255)
        gACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTgg                                                                                                                                              >  4:5025587/1‑149 (MQ=255)
             aaaTCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCt                                                                                                                                         >  3:1150982/1‑149 (MQ=255)
             aaaTCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCt                                                                                                                                         >  3:2622996/1‑149 (MQ=255)
               aTCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTcc                                                                                                                                       >  4:4051904/1‑149 (MQ=255)
                   gTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCg                                                                                                                                     >  2:6511154/1‑147 (MQ=255)
                       gcTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGttcttc                                                                                                                               >  2:323568/1‑149 (MQ=255)
                                tgtTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTa                                                                                                                      <  2:5372392/149‑1 (MQ=255)
                                   ttCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACtt                                                                                                                   >  1:2261414/1‑149 (MQ=255)
                                   ttCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCGCCCGTTCTTCACGCCCGTACtt                                                                                                                   <  4:6438153/149‑1 (MQ=255)
                                     ctcGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTa                                                                                                                 >  4:1341130/1‑149 (MQ=255)
                                     ctcGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTa                                                                                                                 <  3:5387335/149‑1 (MQ=255)
                                      tcGCGTAACCATCTGGAGCGACATCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTAc                                                                                                                <  3:4380682/149‑1 (MQ=255)
                                          gTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGttc                                                                                                                                  <  4:5693658/127‑1 (MQ=255)
                                          gTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGcc                                                                                                            <  1:2445267/149‑1 (MQ=255)
                                          gTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGcc                                                                                                            >  2:3844950/1‑149 (MQ=255)
                                           tAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGTTCCCGTTCTTCACGCCCGTACTTTACAGCCg                                                                                                           <  1:2594629/149‑1 (MQ=255)
                                            aaCCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCgg                                                                                                          <  1:5282339/149‑1 (MQ=255)
                                            aaCCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCgg                                                                                                          >  1:5943382/1‑149 (MQ=255)
                                            aaCCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCgg                                                                                                          >  2:736713/1‑149 (MQ=255)
                                              ccATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGtt                                                                                                        <  2:5621218/149‑1 (MQ=255)
                                                aTCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTc                                                                                                      >  4:4607549/1‑149 (MQ=255)
                                                 tCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCt                                                                                                     <  3:471924/149‑1 (MQ=255)
                                                 tCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCt                                                                                                     <  3:4049265/149‑1 (MQ=255)
                                                 tCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCt                                                                                                     <  1:1512863/149‑1 (MQ=255)
                                                 tCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCt                                                                                                     <  2:1284021/149‑1 (MQ=255)
                                                    ggAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGc                                                                                                  >  4:5855614/1‑149 (MQ=255)
                                                    ggAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGc                                                                                                  >  2:2711690/1‑149 (MQ=255)
                                                            gTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAAt                                                                                          >  1:1040738/1‑149 (MQ=255)
                                                            gTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAAt                                                                                          >  3:938666/1‑149 (MQ=255)
                                                             tCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATc                                                                                         <  3:1665632/149‑1 (MQ=255)
                                                             tCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATc                                                                                         <  4:6198997/149‑1 (MQ=255)
                                                                  ggggTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCt                                                                                    <  1:6332128/149‑1 (MQ=255)
                                                                    ggTATACAGTGCCGACCCGCGTAAAGCAACAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGACGCGCCTGGCGTCTCCCGCTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGc                                                                                  <  2:4972521/149‑1 (MQ=255)
                                                                     gTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCa                                                                                 >  2:944975/1‑149 (MQ=255)
                                                                     gTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCa                                                                                 >  2:4461330/1‑149 (MQ=255)
                                                                         aCAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTg                                                                             >  1:6891036/1‑149 (MQ=255)
                                                                          cAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTgc                                                                            >  1:2675622/1‑149 (MQ=255)
                                                                             tGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCt                                                                         <  4:170823/149‑1 (MQ=255)
                                                                               ccGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCaa                                                                                >  3:6637029/1‑140 (MQ=255)
                                                                               ccGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGt                                                                       >  3:4513015/1‑149 (MQ=255)
                                                                               ccGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGt                                                                       <  4:3729428/149‑1 (MQ=255)
                                                                               ccGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGt                                                                       >  3:6805956/1‑149 (MQ=255)
                                                                                    ccGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTa                                                                  >  2:619400/1‑149 (MQ=255)
                                                                                    ccGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTa                                                                  >  1:1833520/1‑149 (MQ=255)
                                                                                        gTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACg                                                              >  4:2984335/1‑149 (MQ=255)
                                                                                          aaaGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGcc                                                            >  2:6837232/1‑149 (MQ=255)
                                                                                              tGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGAt                                                        <  2:2605994/149‑1 (MQ=255)
                                                                                               gAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATc                                                       >  1:760936/1‑149 (MQ=255)
                                                                                                aaaGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCa                                                      >  2:1761750/1‑149 (MQ=255)
                                                                                                aaaGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCa                                                      >  4:5379569/1‑149 (MQ=255)
                                                                                                aaaGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCa                                                      >  2:2960522/1‑149 (MQ=255)
                                                                                                 acgATGCCTGACTGCTGCCGTTGCTGCGTCTGGATGATCCCAGCTGCCTTGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCaa                                                     <  2:3509719/147‑1 (MQ=255)
                                                                                                 aaGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCaa                                                     <  1:2820679/149‑1 (MQ=255)
                                                                                                 aaGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCaa                                                     <  2:1040738/149‑1 (MQ=255)
                                                                                                 aaGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCaa                                                     >  2:2481070/1‑149 (MQ=255)
                                                                                                  aGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAg                                                    <  4:3853568/149‑1 (MQ=255)
                                                                                                   gATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTg                                                                        >  1:1341440/1‑128 (MQ=255)
                                                                                                     tgcctgccTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGtt                                                 <  4:2597141/149‑1 (MQ=255)
                                                                                                     tgcctgccTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGtt                                                 <  1:400528/149‑1 (MQ=255)
                                                                                                      gcctgcctGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAg                                                                                                              >  1:3289966/1‑87 (MQ=255)
                                                                                                      gcctgcctGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTc                                                <  1:2960522/149‑1 (MQ=255)
                                                                                                        ctgcctgcTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCa                                              >  2:212876/1‑149 (MQ=255)
                                                                                                        ctgcctgcTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCa                                              >  2:3645952/1‑149 (MQ=255)
                                                                                                        ctgcctgcTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGCCGCGCCTGGCGGCTGCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACTCGCCGGATCAAGGTGCCa                                              >  2:7191159/1‑149 (MQ=255)
                                                                                                         tgcctgcTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCAc                                             <  2:411743/149‑1 (MQ=255)
                                                                                                              gctgcCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCa                                        >  3:2058955/1‑149 (MQ=255)
                                                                                                              gctgcCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCa                                        >  2:1500852/1‑149 (MQ=255)
                                                                                                               ctgcCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCg                                                                                                                                     >  4:6691435/1‑55 (MQ=255)
                                                                                                                  ccGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGa                                    >  1:6541893/1‑149 (MQ=255)
                                                                                                                  ccGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGa                                    >  1:3507202/1‑149 (MQ=255)
                                                                                                                  ccGTTGCTGCGTCTCGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGa                                    <  2:6362822/149‑1 (MQ=255)
                                                                                                                   cGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCAGGCCCGTACTTTACAGCCGGTTTCTGGCACCGAAGTCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGAGCATTGaa                                   >  4:4396007/1‑149 (MQ=255)
                                                                                                                    gTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAAc                                  <  1:4635697/149‑1 (MQ=255)
                                                                                                                     ttGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAAcg                                 >  2:3343495/1‑149 (MQ=255)
                                                                                                                     ttGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAAcg                                 >  1:4781913/1‑149 (MQ=255)
                                                                                                                        ctgcGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGt                              <  2:4781913/149‑1 (MQ=255)
                                                                                                                         tgcGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTg                             >  1:2051435/1‑149 (MQ=255)
                                                                                                                         tgcGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTg                             >  4:712817/1‑149 (MQ=255)
                                                                                                                          gcGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGc                            >  1:427650/1‑149 (MQ=255)
                                                                                                                            gTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGTTg                          <  2:574015/149‑1 (MQ=255)
                                                                                                                             tCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGAt                                                        <  3:6397457/118‑1 (MQ=255)
                                                                                                                               tGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGcc                       >  1:1706404/1‑149 (MQ=255)
                                                                                                                                ggATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGa                                                                                       >  3:2299246/1‑84 (MQ=255)
                                                                                                                                 gATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCTGTACTTTACAGCCGGGTTCTGGCAGCGAAATCGACCTTCAACTTCGCTGCAGGTAAGCGCCTGCTCAAGGTGTCAGGCGCATTGAACGCTTGCTGGCCTc                     >  4:1145033/1‑149 (MQ=255)
                                                                                                                                 gATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTc                     >  1:7108397/1‑149 (MQ=255)
                                                                                                                                    gAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAAc                                  <  1:4933753/133‑1 (MQ=255)
                                                                                                                                    gAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCgg                  <  1:6124048/149‑1 (MQ=255)
                                                                                                                                       gCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTcc                                               >  4:2269188/1‑117 (MQ=255)
                                                                                                                                       gCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTAc               >  3:1871110/1‑149 (MQ=255)
                                                                                                                                       gCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTAc               >  3:3206083/1‑149 (MQ=255)
                                                                                                                                        ccAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACt              <  1:3778886/149‑1 (MQ=255)
                                                                                                                                          aGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTgg            >  2:5993879/1‑149 (MQ=255)
                                                                                                                                          aGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTgg            >  3:6075386/1‑149 (MQ=255)
                                                                                                                                               aCTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTgcgc       <  3:5749074/149‑1 (MQ=255)
                                                                                                                                                 tGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGcnngn     <  2:4563213/149‑2 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATT  >  NC_000913/4130468‑4130764

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