Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,855 C→A 5.1% 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 5.1% 381.1 / 3.9 109intergenic (‑82/+36)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (55/48);  new base A (2/4);  total (57/52)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 110 total aligned reads displayed.

TTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAA  >  NC_000913/4366710‑4367003
                                                                                                                                                  |                                                                                                                                                    
ttCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGa                                                                                                                                                   <  2:322201/149‑1 (MQ=255)
 tCCCCCCAATCAGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGAc                                                                                                                                                  <  3:573672/149‑1 (MQ=255)
   cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGt                                                                                                                                                >  2:774710/1‑149 (MQ=255)
    ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                               >  1:72998/1‑149 (MQ=255)
      ccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACCAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGcc                                                                                                                                             <  4:62268/149‑1 (MQ=255)
       cAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCt                                                                                                                                            <  3:120341/149‑1 (MQ=255)
         aTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAATAAAAAAGGCACGTCAGATGACGTGCCtt                                                                                                                                           <  1:796502/149‑1 (MQ=255)
         aTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCttt                                                                                                                                          <  2:141867/149‑1 (MQ=255)
              gcgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                                                                                                                                     >  2:718544/1‑149 (MQ=255)
               cgcCCAAGGAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAATTTGATATTATATTGTTATTTTTAAGTTACTACTCAC‑AAGAAAAATGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                                    <  3:529344/149‑1 (MQ=255)
                gcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAttt                                                                                                                                   <  4:809698/149‑1 (MQ=255)
                gcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAttt                                                                                                                                   >  2:147591/1‑149 (MQ=255)
                  ccAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGt                                                                                                                                 >  1:848936/1‑149 (MQ=255)
                   cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                >  4:733419/1‑149 (MQ=255)
                   cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                >  3:195935/1‑149 (MQ=255)
                    aagaagATCGCCAGCAAAACTATGATTAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtac                                                                                                                               >  4:308806/1‑149 (MQ=255)
                     agaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtact                                                                                                                              >  2:121727/1‑149 (MQ=255)
                        agaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAcc                                                                                                                           <  1:669221/149‑1 (MQ=255)
                             gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTa                                                                                                                      <  1:774710/149‑1 (MQ=255)
                                aGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGa                                                                                                                   >  2:783942/1‑149 (MQ=255)
                                 gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAt                                                                                                                  <  3:273661/149‑1 (MQ=255)
                                 gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAt                                                                                                                  >  1:98941/1‑149 (MQ=255)
                                  cAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAtt                                                                                                                 >  1:466742/1‑149 (MQ=255)
                                    aaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                                    <  1:290134/128‑1 (MQ=255)
                                    aaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                                    >  2:290134/1‑128 (MQ=255)
                                          tgatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCg                                                                                                         >  3:175443/1‑149 (MQ=255)
                                                gTTCTACAACTAGCATATTATCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                                                                                                                                     <  4:80847/115‑1 (MQ=255)
                                                gTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                                                                                                                                     >  3:80847/1‑115 (MQ=255)
                                                 ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                                  >  2:781867/1‑149 (MQ=255)
                                                 ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                                  >  2:495271/1‑149 (MQ=255)
                                                   cTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                <  2:848936/149‑1 (MQ=255)
                                                        aCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtata                                                                                           <  4:195935/149‑1 (MQ=255)
                                                        aCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtata                                                                                           <  4:504959/149‑1 (MQ=255)
                                                          tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtata                                                                                           <  4:852959/147‑1 (MQ=255)
                                                          tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtata                                                                                           >  3:852959/1‑147 (MQ=255)
                                                          tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCC‑TTTTTATTTGTACTAACCTGTACGATTACTGTTCGCTTTCATCAGTATAGCg                                                                                        >  2:509442/1‑149 (MQ=255)
                                                              atatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTCCTGTTCGCTTTCATCAGTATAGCGttt                                                                                     >  4:788087/1‑149 (MQ=255)
                                                                atTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                               >  1:791576/1‑89 (MQ=255)
                                                                atTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                               <  2:791576/89‑1 (MQ=255)
                                                                   aGCCTTCCTTGTTTTTTAACAAGTTGATATTCGATTGTTATTTTTTAGTTTCTGCTCAC‑AAGAAAAAAGGCACGTCATCTAACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCtt                                                                                <  2:526293/149‑1 (MQ=255)
                                                                    gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                               <  3:38794/149‑1 (MQ=255)
                                                                    gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                               <  3:383186/149‑1 (MQ=255)
                                                                    gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                               <  1:329140/149‑1 (MQ=255)
                                                                      cttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTga                                                                             <  2:760158/149‑2 (MQ=255)
                                                                           ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                                  <  2:810635/123‑1 (MQ=255)
                                                                           ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                                  >  1:810635/1‑123 (MQ=255)
                                                                              ttttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCggg                                                                     >  1:772240/1‑149 (MQ=255)
                                                                              ttttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCggg                                                                     <  2:466742/149‑1 (MQ=255)
                                                                                         gTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                               >  2:451247/1‑128 (MQ=255)
                                                                                          ttGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTGCTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                         <  1:781867/149‑1 (MQ=255)
                                                                                          ttGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                         <  1:495271/149‑1 (MQ=255)
                                                                                            gATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTg                                                       >  3:700199/1‑149 (MQ=255)
                                                                                              tatTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTa                                                     >  2:575726/1‑149 (MQ=255)
                                                                                                ttAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                <  2:29592/104‑1 (MQ=255)
                                                                                                ttAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                >  1:29592/1‑104 (MQ=255)
                                                                                                     ttGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGaa                                              <  3:322092/149‑1 (MQ=255)
                                                                                                      tGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGaaa                                             <  4:624259/149‑1 (MQ=255)
                                                                                                        ttATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAat                                           >  2:91470/1‑149 (MQ=255)
                                                                                                                  gTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGtttt                                                                                    <  2:709481/98‑1 (MQ=255)
                                                                                                                  gTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGtttt                                                                                    >  1:709481/1‑98 (MQ=255)
                                                                                                                  gTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTCTCTTCAGTATAGCGTTTTGATTTGTAAGGCGGGTTCATCAGATTCTGTACGGAGAATATATCGTCAAGt                                 >  2:832887/1‑149 (MQ=255)
                                                                                                                   ttACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                   <  4:551238/82‑1 (MQ=255)
                                                                                                                   ttACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                   >  3:551238/1‑82 (MQ=255)
                                                                                                                    tACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTc                               >  1:623106/1‑149 (MQ=255)
                                                                                                                         cTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCtt                          <  1:718544/149‑1 (MQ=255)
                                                                                                                          tCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTc                         >  1:196324/1‑149 (MQ=255)
                                                                                                                           cac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCa                        >  1:536863/1‑149 (MQ=255)
                                                                                                                            ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACGCTGTACGATTACAGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAata                                          >  4:841352/1‑130 (MQ=255)
                                                                                                                            ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAata                                          <  3:841352/130‑1 (MQ=255)
                                                                                                                            ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAg                       >  1:557446/1‑149 (MQ=255)
                                                                                                                            ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAg                       >  1:249303/1‑149 (MQ=255)
                                                                                                                               aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                     <  1:121727/149‑1 (MQ=255)
                                                                                                                               aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                     >  2:538414/1‑149 (MQ=255)
                                                                                                                                 gAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCaa                   <  1:147591/149‑1 (MQ=255)
                                                                                                                                  aaaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATGTCGTCaa                                   <  4:312102/132‑1 (MQ=255)
                                                                                                                                  aaaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCaa                                   >  3:312102/1‑132 (MQ=255)
                                                                                                                                   aaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATCACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTTAGCCGGGTGCATCAGATTCTTTACGGAGGAAATATCGTCAAGTTCCTCTTCAGTTAACa                 >  1:454854/1‑149 (MQ=255)
                                                                                                                                    aaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCtt                                                                                <  3:237060/85‑1 (MQ=255)
                                                                                                                                    aaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCtt                                                                                >  4:237060/1‑85 (MQ=255)
                                                                                                                                    aaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAg                >  1:612595/1‑149 (MQ=255)
                                                                                                                                     aaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGa               >  2:485709/1‑149 (MQ=255)
                                                                                                                                     aaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGa               >  3:867983/1‑149 (MQ=255)
                                                                                                                                      aaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                      <  1:71210/93‑1 (MQ=255)
                                                                                                                                      aaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                      <  1:632324/93‑1 (MQ=255)
                                                                                                                                      aaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                      >  2:71210/1‑93 (MQ=255)
                                                                                                                                      aaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAc              <  1:538414/149‑1 (MQ=255)
                                                                                                                                      aaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGCGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAc              >  1:155191/1‑149 (MQ=255)
                                                                                                                                       aGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCAGCAGTATAGCGTTGTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc             >  2:206603/1‑149 (MQ=255)
                                                                                                                                         gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc           >  4:201056/1‑149 (MQ=255)
                                                                                                                                         gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc           >  2:632002/1‑149 (MQ=255)
                                                                                                                                         gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc           >  2:707515/1‑149 (MQ=255)
                                                                                                                                         gCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc           >  1:252841/1‑149 (MQ=255)
                                                                                                                                          cACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgcg          <  2:557446/149‑1 (MQ=255)
                                                                                                                                            cGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGtt        <  2:98941/149‑1 (MQ=255)
                                                                                                                                             gTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCt                           <  2:646323/129‑1 (MQ=255)
                                                                                                                                             gTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCt                           >  1:646323/1‑129 (MQ=255)
                                                                                                                                             gTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTc       >  1:318821/1‑149 (MQ=255)
                                                                                                                                                aTCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg    >  3:637865/1‑149 (MQ=255)
                                                                                                                                                 tCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATTGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGa   <  4:866773/149‑1 (MQ=255)
                                                                                                                                                  cTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGaa  <  3:476321/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAA  >  NC_000913/4366710‑4367003

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