Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,855 C→A 8.7% intergenic (‑82/+36) 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,8550CA 8.7% 420.9 / 21.1 125intergenic (‑82/+36)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (58/56);  new base A (5/6);  total (63/62)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 128 total aligned reads displayed.

CTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAA  >  NC_000913/4366707‑4367003
                                                                                                                                                     |                                                                                                                                                    
cTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATc                                                                                                                                                      <  3:198669/149‑1 (MQ=255)
 tATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCt                                                                                                                                                     >  4:744667/1‑149 (MQ=255)
 tATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCt                                                                                                                                                     >  3:876153/1‑149 (MQ=255)
 tATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCt                                                                                                                                                     >  2:511738/1‑149 (MQ=255)
  aTTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTg                                                                                                                                                    >  1:227658/1‑149 (MQ=255)
  aTTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTg                                                                                                                                                    >  2:10273/1‑149 (MQ=255)
     ccccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAATAAAAAAGGCACGTCAGATGAc                                                                                                                                                  >  3:214915/1‑149 (MQ=255)
      cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGt                                                                                                                                                >  3:283728/1‑149 (MQ=255)
       ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                               >  4:288713/1‑149 (MQ=255)
        cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGc                                                                                                                                              >  1:188109/1‑149 (MQ=255)
        cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGc                                                                                                                                              >  3:459068/1‑149 (MQ=255)
           aaTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCtt                                                                                                                                           >  1:141386/1‑149 (MQ=255)
            aTCTGGCGCCCAAGAAGATCGCCAGCAAAGCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTTTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCttt                                                                                                                                          <  4:533039/149‑1 (MQ=255)
             tCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCtttt                                                                                                                                         <  1:276414/149‑1 (MQ=255)
              cTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCttttt                                                                                                                                        <  4:267825/149‑1 (MQ=255)
               tGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCtttttt                                                                                                                                       >  1:758929/1‑149 (MQ=255)
                ggCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTa                                                                                                                                      <  2:141386/149‑1 (MQ=255)
                  cgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                                    <  2:493985/149‑1 (MQ=255)
                  cgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCTTTTTTCtt                                                                                                                                    <  4:327311/149‑1 (MQ=255)
                   gcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAttt                                                                                                                                   >  2:227926/1‑149 (MQ=255)
                    cccAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAATAAAAAAGGCACGTCAGATGACGTGCCTTTTTTCttg                                                                                                                                   <  4:214915/149‑2 (MQ=38)
                     ccAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGt                                                                                                                                 >  1:826284/1‑149 (MQ=255)
                      cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCttttttc                                                                                                                                      <  3:757602/143‑2 (MQ=255)
                      cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCttttttc                                                                                                                                      >  4:757602/1‑142 (MQ=255)
                         gaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTTAC‑AATGAAAAAGACACGTCATCTTACTTGCCTTTTTTATTTGtacta                                                                                                                             >  2:349921/1‑149 (MQ=255)
                          aagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtactac                                                                                                                            >  1:99479/1‑149 (MQ=255)
                          aagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtactac                                                                                                                            <  2:815710/149‑1 (MQ=255)
                             aTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCt                                                                                                                         <  1:981599/149‑1 (MQ=255)
                                gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTa                                                                                                                      <  1:504532/149‑1 (MQ=255)
                                    gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAt                                                                                                                  >  1:704560/1‑149 (MQ=255)
                                     cAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAtt                                                                                                                 <  4:461750/149‑1 (MQ=255)
                                      aaaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTa                                                                                                                >  4:129235/1‑149 (MQ=255)
                                        aaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACt                                                                                                              <  4:836213/149‑1 (MQ=255)
                                              gatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGc                                                                                                        <  4:80182/149‑1 (MQ=255)
                                               atgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCt                                                                                                       <  1:511738/149‑1 (MQ=255)
                                                tgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCtt                                                                                                      >  4:254795/1‑149 (MQ=255)
                                                 gaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCttt                                                                                                     <  3:129235/149‑1 (MQ=255)
                                                  aGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtc                                                                                                    <  4:876153/149‑1 (MQ=255)
                                                   gTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                   <  1:227926/149‑1 (MQ=255)
                                                      cTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                >  4:605751/1‑149 (MQ=255)
                                                       tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                               >  2:470117/1‑149 (MQ=255)
                                                             tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGc                                                                                         >  1:574756/1‑149 (MQ=255)
                                                             tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGc                                                                                         >  2:408978/1‑149 (MQ=255)
                                                             tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGc                                                                                         >  2:1065242/1‑149 (MQ=255)
                                                               gCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGt                                                                                       >  3:627423/1‑149 (MQ=255)
                                                                     tAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCt                                                                                 <  4:111336/149‑1 (MQ=255)
                                                                       gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCttt                                                                                                     <  3:739328/127‑1 (MQ=255)
                                                                       gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                               <  1:687796/149‑1 (MQ=255)
                                                                       gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                               <  3:280615/149‑1 (MQ=255)
                                                                          ttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCATTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTa                                                                            >  2:880524/1‑149 (MQ=255)
                                                                              ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGcc                                                                        >  2:815626/1‑149 (MQ=255)
                                                                              ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGcc                                                                        >  3:482644/1‑149 (MQ=255)
                                                                                gTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                      <  4:950600/149‑1 (MQ=255)
                                                                                gTTTTTTAAAAAATTTATATCCGATTTTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                      <  2:260588/149‑1 (MQ=255)
                                                                                 ttttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCggg                                                                     <  4:459068/149‑1 (MQ=255)
                                                                                  tttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGt                                                                    <  2:269396/149‑1 (MQ=255)
                                                                                      tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCAt                                                                <  3:795507/149‑1 (MQ=255)
                                                                                       aaCAAGTTGATATTAGATTGTTATTTTTCAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATc                                                               <  1:563482/149‑1 (MQ=255)
                                                                                       aaCAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATc                                                               >  3:221557/1‑149 (MQ=255)
                                                                                          aaGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGa                                                            <  4:390924/149‑1 (MQ=255)
                                                                                            gTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGAtt                                                          >  1:266280/1‑149 (MQ=255)
                                                                                             ttGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAATGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                         >  3:20291/1‑149 (MQ=255)
                                                                                               gATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTg                                                       <  4:627423/149‑1 (MQ=255)
                                                                                               gATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTg                                                       <  4:456002/149‑1 (MQ=255)
                                                                                                        ttGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGaa                                              >  4:1043189/1‑149 (MQ=255)
                                                                                                           ttATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCt                                                        <  3:442126/136‑1 (MQ=255)
                                                                                                           ttATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAat                                           >  4:405686/1‑149 (MQ=255)
                                                                                                            tATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAata                                          <  3:836886/149‑1 (MQ=255)
                                                                                                                tttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGt                                      <  1:359684/149‑1 (MQ=255)
                                                                                                                 ttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                     <  1:880524/149‑1 (MQ=255)
                                                                                                                 ttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                     <  2:568486/149‑1 (MQ=255)
                                                                                                                   aaGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCaa                                   <  2:761141/149‑1 (MQ=255)
                                                                                                                       ttcTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTc                               <  1:863170/147‑1 (MQ=255)
                                                                                                                          tGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGTAGAAAATATCGTCAAGTTCCGc                            <  2:934204/149‑1 (MQ=255)
                                                                                                                          tGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGc                            <  2:227658/149‑1 (MQ=255)
                                                                                                                             tCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTc                         >  1:619821/1‑149 (MQ=255)
                                                                                                                               ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAg                       >  3:261329/1‑149 (MQ=255)
                                                                                                                               ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAg                       >  3:848568/1‑149 (MQ=255)
                                                                                                                                  aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                     >  1:619699/1‑149 (MQ=255)
                                                                                                                                  aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                     <  3:288713/149‑1 (MQ=255)
                                                                                                                                    gAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTc                               >  1:218402/1‑137 (MQ=255)
                                                                                                                                     aaaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAAc                  >  3:622795/1‑149 (MQ=255)
                                                                                                                                     aaaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAAc                  >  4:546198/1‑149 (MQ=255)
                                                                                                                                       aaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAg                >  4:782452/1‑149 (MQ=255)
                                                                                                                                       aaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAg                >  2:927492/1‑149 (MQ=255)
                                                                                                                                          aGGCACGTCGGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc             >  3:733410/1‑149 (MQ=255)
                                                                                                                                          aGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATGGTCAAGTTCCGCTTCa                        >  1:845678/1‑138 (MQ=255)
                                                                                                                                          aGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCa                        <  2:845678/138‑1 (MQ=255)
                                                                                                                                            gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGt                                                                    >  3:885257/1‑92 (MQ=255)
                                                                                                                                            gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGt                                                                    <  4:885257/92‑1 (MQ=255)
                                                                                                                                            gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTTTGTACGGAGAAAATATc                                        >  2:354880/1‑120 (MQ=255)
                                                                                                                                            gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc           >  3:899024/1‑149 (MQ=255)
                                                                                                                                            gCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGc                                                                                                        <  1:1048526/56‑1 (MQ=37)
                                                                                                                                              aCGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGtt                                <  3:100062/126‑1 (MQ=255)
                                                                                                                                              aCGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGtt                                >  4:100062/1‑126 (MQ=255)
                                                                                                                                                gTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTc       >  2:116245/1‑149 (MQ=255)
                                                                                                                                                gTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTc       >  3:160352/1‑149 (MQ=255)
                                                                                                                                                 tCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTcc      <  4:813828/149‑1 (MQ=255)
                                                                                                                                                   aTCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg    >  2:171953/1‑149 (MQ=255)
                                                                                                                                                     cTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGaa  <  4:829428/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAA  >  NC_000913/4366707‑4367003

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