Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,854 T→G 8.7% intergenic (‑81/+37) dcuA ← / ← aspA C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,366,8540TG 8.7% 326.3 / 22.3 103intergenic (‑81/+37)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base T (52/42);  new base G (3/6);  total (55/48)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 103 total aligned reads displayed.

TATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCA  >  NC_000913/4366708‑4367000
                                                                                                                                                   |                                                                                                                                                  
tATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATc                                                                                                                                                   >  1:921137/1‑148 (MQ=255)
 aTTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCt                                                                                                                                                  >  1:636236/1‑148 (MQ=255)
   tCCCCCCAATCTGGCGCCCAAGAAGATCGCTTGCAAAACTATGATGAGTCCTACATCTAGCATATTAGCCTTCCTTTTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGa                                                                                                                                                <  1:305308/148‑1 (MQ=255)
      ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACg                                                                                                                                              <  1:543555/147‑1 (MQ=255)
      ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                            >  3:116168/1‑149 (MQ=255)
      ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                            >  3:202518/1‑149 (MQ=255)
       cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                            >  1:330309/1‑148 (MQ=255)
       cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGc                                                                                                                                           >  3:177156/1‑149 (MQ=255)
           aTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAATAAAAAAGGCACGTCAGATGACGTGCCtt                                                                                                                                        <  4:85552/149‑1 (MQ=255)
             cTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCtttt                                                                                                                                      <  2:256805/148‑1 (MQ=255)
             cTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCtttt                                                                                                                                      <  2:20344/148‑1 (MQ=255)
                 cgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                                 <  4:177156/149‑1 (MQ=255)
                  gcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                                 <  2:271209/148‑1 (MQ=255)
                    ccAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AATAAAAAAGGCACGTCAGATGACGTGCCTTTTTTCttgt                                                                                                                               <  2:711862/148‑3 (MQ=39)
                       agaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGATTTTTAACAAGTTGATAGTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTCTTTTATTTGtaca                                                                                                                           <  4:620310/149‑2 (MQ=255)
                          agaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtactac                                                                                                                         >  1:875261/1‑148 (MQ=255)
                           gaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAccc                                                                                                                       >  3:247289/1‑149 (MQ=255)
                               gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTa                                                                                                                   >  4:455936/1‑149 (MQ=255)
                                   gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTCTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAt                                                                                                               >  4:240858/1‑149 (MQ=255)
                                      aaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTa                                                                                                             <  1:536319/148‑1 (MQ=255)
                                          tATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGt                                                                                                         <  2:172611/148‑1 (MQ=255)
                                              atgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTCTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCt                                                                                                    <  3:240858/149‑1 (MQ=255)
                                                  gTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                <  3:455936/149‑1 (MQ=255)
                                                    tCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                               >  2:707350/1‑148 (MQ=255)
                                                     cTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatc                                                                                              <  2:273621/148‑1 (MQ=255)
                                                      tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                            >  3:120175/1‑149 (MQ=255)
                                                       aCAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                            <  1:673699/148‑1 (MQ=255)
                                                       aCAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGt                                                                                           <  3:422621/149‑1 (MQ=255)
                                                            tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAg                                                                                       <  2:921137/148‑1 (MQ=255)
                                                             aGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCg                                                                                     >  3:299195/1‑149 (MQ=255)
                                                                   ttAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGc                                                                               >  4:185882/1‑149 (MQ=255)
                                                                   ttAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGc                                                                               >  4:865385/1‑149 (MQ=255)
                                                                    tAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCt                                                                              <  3:512087/149‑1 (MQ=255)
                                                                       ccttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                            >  1:343141/1‑148 (MQ=255)
                                                                       ccttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTg                                                                           <  4:831066/149‑1 (MQ=255)
                                                                        cttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGt                                                                          <  4:299195/149‑1 (MQ=255)
                                                                          tccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGTCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTaa                                                                        >  4:648675/1‑149 (MQ=255)
                                                                          tccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTa                                                                         <  2:636236/148‑1 (MQ=255)
                                                                             ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAACAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGTTTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGc                                                                      >  1:791697/1‑148 (MQ=255)
                                                                             ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGcc                                                                     >  3:282896/1‑149 (MQ=255)
                                                                              tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGcc                                                                     >  2:668445/1‑148 (MQ=255)
                                                                              tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCg                                                                    <  3:554729/149‑1 (MQ=255)
                                                                               gTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                   >  4:224924/1‑149 (MQ=255)
                                                                                 tttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCggg                                                                  <  1:707350/148‑1 (MQ=255)
                                                                                  ttttAACAAGTTGCTATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTg                                                                >  4:776130/1‑149 (MQ=255)
                                                                                     tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCa                                                              <  2:791697/148‑1 (MQ=255)
                                                                                     tAACAAGTTGATATTAGATTGTTATTTATAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCTTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCa                                                              <  1:266202/148‑1 (MQ=255)
                                                                                       aCAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCa                                                           <  3:224924/149‑1 (MQ=255)
                                                                                         aaGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAg                                                          <  2:343141/148‑1 (MQ=255)
                                                                                           gTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGAtt                                                       >  3:133883/1‑149 (MQ=255)
                                                                                            ttGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGAtt                                                       >  2:501440/1‑148 (MQ=255)
                                                                                             tGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                      >  2:260953/1‑148 (MQ=255)
                                                                                                tatTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGt                                                                                           <  2:819964/108‑1 (MQ=255)
                                                                                                tatTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGt                                                                                           >  1:819964/1‑108 (MQ=255)
                                                                                                tatTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTa                                                  <  4:120175/149‑1 (MQ=255)
                                                                                                     gATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGgag                                             <  4:877978/149‑1 (MQ=255)
                                                                                                      aTTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGgaga                                            >  3:239620/1‑149 (MQ=255)
                                                                                                         gTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGaaa                                          <  2:875261/148‑1 (MQ=255)
                                                                                                           tATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAata                                       <  4:675619/149‑1 (MQ=255)
                                                                                                           tATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAat                                        <  2:590257/148‑1 (MQ=255)
                                                                                                             tttttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAatat                                      <  2:350671/148‑1 (MQ=255)
                                                                                                             tttttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATc                                     <  4:133883/149‑1 (MQ=255)
                                                                                                              ttttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATc                                     >  2:257889/1‑148 (MQ=255)
                                                                                                              ttttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATc                                     >  1:635309/1‑148 (MQ=255)
                                                                                                               tttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCg                                    <  2:36947/148‑1 (MQ=255)
                                                                                                                ttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTTTGTACGGAGAAAATATCGt                                   >  2:280067/1‑148 (MQ=255)
                                                                                                                   aGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAg                               <  4:892216/149‑1 (MQ=255)
                                                                                                                       aCTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTcc                           >  3:25700/1‑149 (MQ=255)
                                                                                                                         tGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCg                          <  2:134492/148‑1 (MQ=255)
                                                                                                                         tGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGc                         <  3:312107/149‑1 (MQ=255)
                                                                                                                            tCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatc                                                                                              <  3:632391/77‑1 (MQ=255)
                                                                                                                            tCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatc                                                                                              >  4:632391/1‑77 (MQ=255)
                                                                                                                            tCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTc                      >  4:462134/1‑149 (MQ=255)
                                                                                                                             cac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTc                      >  1:453773/1‑148 (MQ=255)
                                                                                                                             cac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTc                      >  2:785833/1‑148 (MQ=255)
                                                                                                                              ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAg                    >  4:533244/1‑149 (MQ=255)
                                                                                                                               c‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAg                    >  1:206829/1‑148 (MQ=255)
                                                                                                                               c‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAg                    >  1:381619/1‑148 (MQ=255)
                                                                                                                               c‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGt                   <  3:95954/149‑1 (MQ=255)
                                                                                                                                 aataaaaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCATATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                  >  4:356486/4‑149 (MQ=255)
                                                                                                                                 aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                   <  3:683067/100‑1 (MQ=255)
                                                                                                                                 aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                   >  4:683067/1‑100 (MQ=255)
                                                                                                                                   gAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGAt                                                        >  4:263532/1‑109 (MQ=255)
                                                                                                                                   gAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGAt                                                        <  3:263532/109‑1 (MQ=255)
                                                                                                                                    aaaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTTGCTTTCATCAGTATAGCGTTTTCCTTTGTAAGCAAGGTGCATCAGATTCTGTACGGAGAAAATATCTTCGAGTTCCGCTTCAGTCaa                >  2:436207/1‑148 (MQ=255)
                                                                                                                                    aaaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAAc               >  3:414491/1‑149 (MQ=255)
                                                                                                                                    aaaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAAc               >  4:878774/1‑149 (MQ=255)
                                                                                                                                     aaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACa              <  3:185882/149‑1 (MQ=255)
                                                                                                                                      aaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCttt                                                                                                  >  1:324612/1‑64 (MQ=255)
                                                                                                                                      aaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCttt                                                                                                  <  2:324612/64‑1 (MQ=25)
                                                                                                                                       aaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAg             <  2:545010/148‑1 (MQ=255)
                                                                                                                                        aaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCa                                                              >  4:15634/1‑98 (MQ=255)
                                                                                                                                        aaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCa                                                              <  3:15634/98‑1 (MQ=255)
                                                                                                                                           gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGGTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc        >  3:650292/1‑149 (MQ=255)
                                                                                                                                           gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc        >  3:326014/1‑149 (MQ=255)
                                                                                                                                           gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc        >  3:123323/1‑149 (MQ=255)
                                                                                                                                             aCGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGt      >  4:510878/1‑149 (MQ=255)
                                                                                                                                             aCGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGt      >  3:189601/1‑149 (MQ=255)
                                                                                                                                                tCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTc    >  1:139043/1‑148 (MQ=255)
                                                                                                                                                 cATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCa  >  4:819611/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCA  >  NC_000913/4366708‑4367000

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