Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,130,614 G→T 8.4% L260L (CTG→CTT 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,6140GT 8.4% 1984.2 / 121.6 576L260L (CTG→CTTmetLBifunctional aspartokinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base G (275/252);  new base T (23/26);  total (298/278)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-03

Only 100 of 585 total aligned reads displayed.

ACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTA  >  NC_000913/4130466‑4130762
                                                                                                                                                    |                                                                                                                                                    
  tATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGgc                                                                                                                                                    >  1:5596299/1‑149 (MQ=255)
  tATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGgc                                                                                                                                                    <  2:5677053/149‑1 (MQ=255)
   aTTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGgcg                                                                                                                                                   <  2:7248785/149‑1 (MQ=255)
   aTTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGgcg                                                                                                                                                   >  1:4272923/1‑149 (MQ=255)
    ttCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGgcgc                                                                                                                                                  >  1:6482676/1‑149 (MQ=255)
      ccGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGgcgcgc                                                                                                                                                <  4:56979/149‑1 (MQ=255)
      ccGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGgcgcgc                                                                                                                                                <  3:3788106/149‑1 (MQ=255)
       cgcgACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGcc                                                                                                                                               >  3:6477611/1‑149 (MQ=255)
                    ggTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGt                                                                                                                                  >  2:1329510/1‑149 (MQ=255)
                     gTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGtt                                                                                                                                 >  1:5322341/1‑149 (MQ=255)
                         gcTGGCGGGTGTTTCTCGCGTAACCATCTGAAGCGACGTCGCCGGAGTATACAGTGCCGACCCGCGTAAAGTGATAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGGCTGGCGGCTCCCGttcttc                                                                                                                             <  2:2839012/149‑1 (MQ=255)
                               ggggGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGccc                                                                                                                       >  2:6602352/1‑149 (MQ=255)
                                 gtgtTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGt                                                                                                                     >  3:3305915/1‑149 (MQ=255)
                                     ttCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACtt                                                                                                                 <  4:4252335/149‑1 (MQ=255)
                                     ttCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACtt                                                                                                                 >  2:5914122/1‑149 (MQ=255)
                                     ttCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACtt                                                                                                                 >  1:2261414/1‑149 (MQ=255)
                                       ctcGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTa                                                                                                               >  4:1341130/1‑149 (MQ=255)
                                       ctcGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTa                                                                                                               <  3:5387335/149‑1 (MQ=255)
                                       ctcGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTa                                                                                                               >  2:612628/1‑149 (MQ=255)
                                        tcGCGTAACCATCTGGAGCGACATCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTAc                                                                                                              <  3:4380682/149‑1 (MQ=255)
                                        tcGCGTAACCATCTGCAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGCCTCCCGTTCTTCACGCCCGTACTTTAc                                                                                                              >  2:3719126/1‑149 (MQ=255)
                                         cgcgTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACa                                                                                                             <  3:238203/149‑1 (MQ=255)
                                              aaCCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCgg                                                                                                        <  2:6992409/149‑1 (MQ=255)
                                              aaCCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCgg                                                                                                        >  1:5943382/1‑149 (MQ=255)
                                                ccATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCgg                                                                                                        <  2:1191937/147‑1 (MQ=255)
                                                 cATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGttt                                                                                                     >  4:4470201/1‑149 (MQ=255)
                                                 cATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGttt                                                                                                     >  3:1069520/1‑149 (MQ=255)
                                                  aTCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGATCCCGTTCTCCACGCCCGTACTTTACAGCCGGTTTc                                                                                                    <  4:6738351/149‑1 (MQ=255)
                                                   tCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCt                                                                                                   <  1:1512863/149‑1 (MQ=255)
                                                   tCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCt                                                                                                   <  3:471924/149‑1 (MQ=255)
                                                     tGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTgg                                                                                                 <  1:2959904/149‑1 (MQ=255)
                                                     tGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTgg                                                                                                 <  4:3303976/149‑1 (MQ=255)
                                                      ggAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGc                                                                                                >  1:4001229/1‑149 (MQ=255)
                                                             cGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGaaa                                                                                         <  2:5078347/149‑1 (MQ=255)
                                                              gTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAAt                                                                                        >  1:1040738/1‑149 (MQ=255)
                                                              gTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAAt                                                                                        >  3:938666/1‑149 (MQ=255)
                                                               tCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATc                                                                                       <  1:3743276/149‑1 (MQ=255)
                                                                 gCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGTCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGa                                                                                     >  2:4841518/1‑149 (MQ=255)
                                                                  ccGGGGTATACAGTGCCGACCCGCGTAAAGTGAAATATGCCTGTCTGCCGGTGTTGGTACGTCTGGATGAGGCGAGCGAACTGGCGGTCCTGGGGGCACCCGTGCTACACGCTCGTACTTTATAGCCGTTTTCTGGTAGCGAAATCGAc                                                                                    <  4:3914544/149‑1 (MQ=255)
                                                                  ccGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGAc                                                                                    >  1:4939774/1‑149 (MQ=255)
                                                                    ggggTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCt                                                                                  <  1:6332128/149‑1 (MQ=255)
                                                                     gggTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTg                                                                                 <  1:2640340/149‑1 (MQ=255)
                                                                       gTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCa                                                                               >  3:4131856/1‑149 (MQ=255)
                                                                           aCAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTg                                                                           >  1:6891036/1‑149 (MQ=255)
                                                                              gTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACa                                                                                                             <  3:2150534/112‑1 (MQ=255)
                                                                                 ccGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCaa                                                                              >  3:6637029/1‑140 (MQ=255)
                                                                                 ccGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGt                                                                     >  3:4513015/1‑149 (MQ=255)
                                                                                 ccGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGt                                                                     >  3:6805956/1‑149 (MQ=255)
                                                                                  cGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTa                                                                    <  4:549995/149‑1 (MQ=255)
                                                                                   gACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAg                                                                   <  3:1736948/149‑1 (MQ=255)
                                                                                     cccGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGcc                                                                                                          <  3:2677859/108‑1 (MQ=255)
                                                                                      ccGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTa                                                                <  4:6439921/149‑1 (MQ=255)
                                                                                      ccGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTa                                                                >  1:1833520/1‑149 (MQ=255)
                                                                                      ccGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTa                                                                >  3:1355546/1‑149 (MQ=255)
                                                                                         cgTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGAc                                                                                    >  1:2564637/1‑126 (MQ=255)
                                                                                         cgTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTacac                                                             <  2:6482676/149‑1 (MQ=255)
                                                                                            aaaGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGa                                                                                     >  1:7138084/1‑122 (MQ=255)
                                                                                               gTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGa                                                       >  4:2671624/1‑149 (MQ=255)
                                                                                               gTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGa                                                       >  2:3169911/1‑149 (MQ=255)
                                                                                               gTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGa                                                       >  4:5991347/1‑149 (MQ=255)
                                                                                                tGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGAt                                                      <  2:3860838/149‑1 (MQ=255)
                                                                                                   aaGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCaa                                                   <  2:2620663/149‑1 (MQ=255)
                                                                                                   aaGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCaa                                                   <  4:1786703/149‑1 (MQ=255)
                                                                                                     gATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTg                                                                      >  1:1341440/1‑128 (MQ=255)
                                                                                                     gATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAgg                                                 <  2:4755350/149‑1 (MQ=255)
                                                                                                     gATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAgg                                                 >  1:778885/1‑149 (MQ=255)
                                                                                                        gcctgcctGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAg                                                                                                            <  2:3289966/87‑1 (MQ=255)
                                                                                                        gcctgcctGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAg                                                                                                            >  1:3289966/1‑87 (MQ=255)
                                                                                                          ctgcctgcTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCa                                            >  3:5760367/1‑149 (MQ=255)
                                                                                                              ctgctgCCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCAcgcg                                        <  2:4299465/149‑1 (MQ=255)
                                                                                                              ctgctgCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCAcgcg                                        >  2:631741/1‑149 (MQ=255)
                                                                                                               tgctgcCGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACgcgc                                       <  1:2816074/149‑1 (MQ=255)
                                                                                                                  tgcCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCAtt                                    <  4:51102/149‑1 (MQ=255)
                                                                                                                  tgcCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCAtt                                    >  4:3731520/1‑149 (MQ=255)
                                                                                                                    ccGTTGCTGCGTCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGa                                  >  1:6541893/1‑149 (MQ=255)
                                                                                                                    ccGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGa                                  >  1:3507202/1‑149 (MQ=255)
                                                                                                                      gTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAAc                                <  1:4635697/149‑1 (MQ=255)
                                                                                                                       ttGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAAcg                               >  1:4781913/1‑149 (MQ=255)
                                                                                                                        tgctgcGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAAcgc                              <  2:5548038/149‑1 (MQ=255)
                                                                                                                          ctgcGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGt                            <  2:2179559/149‑1 (MQ=255)
                                                                                                                           tgcGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACa                                                                                                             <  4:1981570/67‑1 (MQ=255)
                                                                                                                           tgcGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTg                           >  4:712817/1‑149 (MQ=255)
                                                                                                                           tgcGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTg                           >  1:2051435/1‑149 (MQ=255)
                                                                                                                              gTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTgcgc                                                                        <  4:6237426/101‑1 (MQ=255)
                                                                                                                               tCTGGATGAGGCCAGCGACCTTGAGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTgg                       >  2:828955/1‑149 (MQ=255)
                                                                                                                               tCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGAt                                                      <  3:6397457/118‑1 (MQ=255)
                                                                                                                                cTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGc                      >  2:3448306/1‑149 (MQ=255)
                                                                                                                                 tGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGaa                                 >  2:2487325/1‑137 (MQ=255)
                                                                                                                                 tGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGcc                     >  1:1706404/1‑149 (MQ=255)
                                                                                                                                  ggATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGa                                                                                     <  4:2299246/84‑1 (MQ=255)
                                                                                                                                   gATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCTGTACTTTACAGCCGGGTTCTGGCAGCGAAATCGACCTTCAACTTCGCTGCAGGTAAGCGCCTGCTCAAGGTGTCAGGCGCATTGAACGCTTGCTGGCCTc                   >  4:1145033/1‑149 (MQ=255)
                                                                                                                                      gAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCgg                <  1:6124048/149‑1 (MQ=255)
                                                                                                                                      gAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCgg                >  2:354418/1‑149 (MQ=255)
                                                                                                                                        ggCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTa              >  2:158078/1‑149 (MQ=255)
                                                                                                                                         gCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTAc             <  1:2444704/149‑1 (MQ=255)
                                                                                                                                           cAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTg           <  2:1993350/149‑1 (MQ=255)
                                                                                                                                            aGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTgg          >  3:6075386/1‑149 (MQ=255)
                                                                                                                                             gCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGt         >  2:4555064/1‑149 (MQ=255)
                                                                                                                                               gAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTgc       >  4:4746288/1‑149 (MQ=255)
                                                                                                                                                aaCTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTgcg      >  3:1890386/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTATTCCGCGACACAAATCGGTGCGCTGGCGGGTGTTTCTCGCGTAACCATCTGGAGCGACGTCGCCGGGGTATACAGTGCCGACCCGCGTAAAGTGAAAGATGCCTGCCTGCTGCCGTTGCTGCGTCTGGATGAGGCCAGCGAACTGGCGCGCCTGGCGGCTCCCGTTCTTCACGCCCGTACTTTACAGCCGGTTTCTGGCAGCGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTA  >  NC_000913/4130466‑4130762

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