Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 809,217 2 bp→AG coding (40‑41/1290 nt) bioA ← adenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913809,2170CA97.6% 128.5 / NA 43W14M (TGG→CTG) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (16/25);  total (17/26)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913809,2180AG100.0% 148.0 / NA 43W14M (TGG→CTG) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/26);  total (17/26)

CAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTC  >  NC_000913/808990‑809455
                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                             
cAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAgg                                                                                                                                                                                                                           <  1:173570/249‑1 (MQ=255)
    ggctggcGCATGGGTGATACCGCCAAGCATCACATGCGAAATGTCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGGCAGACAAAAG‑CAGCTCGCAACCTTCGGCGCTCACCGCCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGAt                                                                                                                                                                                                                     <  2:99258/251‑1 (MQ=255)
        ggcGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTc                                                                                                                                                                                                                 >  1:96757/1‑251 (MQ=255)
               gggTGATTCCGCCAAACATCCCAGGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCCCCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTc                                                                                                                                                                                                          <  2:297807/251‑1 (MQ=255)
                 gTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGGGTATGGGTGCAGGATATGGCGTTGGTCACAGGCAAGATCGTCCGTTGTCAt                                                                                                                                                                                                        >  1:240341/1‑251 (MQ=255)
                   gATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATa                                                                                                                                                                                                       <  1:135553/250‑1 (MQ=255)
                      aCCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATc                                                                                                                                                                                                    <  1:79918/250‑1 (MQ=255)
                      aCCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCg                                                                                                                                                                                                   <  1:146226/251‑1 (MQ=255)
                          ccAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGCGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACtt                                                                                                                                                                                               >  1:330426/1‑251 (MQ=255)
                                      tGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTACGTTGTCATAATCGACTTGTAAACCAAAtt                                                                                                                                                                                   >  1:175017/1‑251 (MQ=255)
                                       gCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAAtt                                                                                                                                                                                   <  1:63406/250‑1 (MQ=255)
                                           cATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGCTCGCCGCCCACCAGGACGACATACCGTCACCCAGGCGTCGGCCGTCAGACAAAAG‑CAGCTCGCAACCTTCGGCGCTCCCCACCGGCTAAACCGGCAGAGGGGCGGTCATGGGTGTGTATGGGTGCCGGCTATGGCGTTGGTCAACGGCAAGATCGTCCTTTGTCGTAATCGACTTGTAAACCAAATTggaa                                                                                                                                                                               >  2:1348/1‑247 (MQ=255)
                                                    aaTTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAACCCAAATTGGAAAGATTTAGGt                                                                                                                                                                     >  1:243781/1‑251 (MQ=255)
                                                                   cgccgcATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAGGATTTAGGTTTACAAGTCTACAcc                                                                                                                                                      <  1:149453/251‑1 (MQ=255)
                                                                        cATTAAGCTGCGGTTGATTGTAGCCGGGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCGACACCGAAtt                                                                                                                                                 <  2:147823/251‑1 (MQ=255)
                                                                              gctccGGGTGATTGTAACCGTGGACCGCCGCCCACCAGGACGACATACCGTCAACCACGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTCCGGCGCTCACCACCGGATAAACCGGAAGCGGGGAGGTCATGGATGCGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGCCTTGTAAACCAAATTGAAAAGATTTAGGTTAACAAGTCTACACCGAATTcacaac                                                                                                                                           >  2:57359/5‑251 (MQ=255)
                                                                                 gCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACaaa                                                                                                                                        <  1:32305/251‑1 (MQ=255)
                                                                                          tGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGtt                                                                                                                               >  1:232542/1‑251 (MQ=255)
                                                                                           gTAGCCGTGGATCGCCGCCCACCAGGACGTCATACCCTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCTCAACCTTCGGCGCTCACCACCGGAGAATCCGGCAGAGGGTTGGTCATGGATGTGTATGGGTGCAGGATATTGCGTTTGTCAAAGGCATGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGACTTAGGTTTACAAGTCTACACCGAATTAACAACTAAAAACACGttt                                                                                                                              <  1:205972/251‑1 (MQ=255)
                                                                                             ggCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTg                                                                                                                            <  2:330426/250‑1 (MQ=255)
                                                                                                       cgccgcCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAG‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGccc                                                                                                                   <  1:256676/250‑1 (MQ=255)
                                                                                                           gccCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAaca                                                                                                                                   <  1:310826/230‑1 (MQ=255)
                                                                                                           gaccaccaGGCCGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACAGGATAAACCGGCAGAGGAGAGGACCTGGATGTGGATGGGTGCAGGGTATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTAATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTCAAAGTCAACACCGCATTAACAACAAAAAacaagaccggaagaacgtcgtga                                                                                                               >  2:310826/3‑230 (MQ=255)
                                                                                                                    gacgacATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGcc                                                                                                     >  1:234501/1‑251 (MQ=255)
                                                                                                                      cgacATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGcc                                                                                                                    <  1:202487/234‑1 (MQ=255)
                                                                                                                      cgacATACAGTCAGCCAGGAGTCTGCCGTCAGACAAAAT‑CAGCCCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGCGGGGCGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATACTCGCCTTGTAAACCAACTTGAAACGATTTCGGTTTACAAGTATACACCGCATTAACAACAAAAAACACGTTTTGGAGAAGcc                                                                                                                    >  2:202487/1‑234 (MQ=255)
                                                                                                                            aCCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTgg                                                                                             <  1:92267/251‑1 (MQ=255)
                                                                                                                                ccAACCAGGCGGCAGCCGTCAGACAAAAG‑CAGCTCCCAGCCTTCGGCGCGCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATGGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTGTTGGAGAAGCCCCATGGCTCACCGCCCCCGCTGGACag                                                                                         <  2:74479/250‑2 (MQ=255)
                                                                                                                                 cAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAAACTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGGCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACAt                                                                                         <  2:175017/250‑1 (MQ=255)
                                                                                                                                           acTGCAGTCAGACAAAATACA‑CTCGCAACCTCCGGTGCTCACCCCCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGAACGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGCCATTGTCGCAAGTc                                                                              <  2:191184/250‑1 (MQ=255)
                                                                                                                                            cTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTg                                                                                                                            >  1:169676/1‑204 (MQ=255)
                                                                                                                                            cTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTg                                                                                                                            <  2:169676/204‑1 (MQ=255)
                                                                                                                                                    aGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAatt                                                                      >  1:270636/1‑250 (MQ=255)
                                                                                                                                                     gtgaaaaT‑CTGCTCGCAACCGTCGGCGCGTACCCCATGATTATCCGGCAGTGTGGAGGTCATGGATGTGTATGGGGTCTGGTTATTGCGTAGGTCATAGGTAAGATCGTACGTTTTCATAATCGACTTGTAAACCAAAGTGATAAGATTGAGGTTTACAAGTCTACACCGAATTAACAAAAAAAAACACTTTTTTGATAATCCCCATGGCTCACCGCCCACGTTGGACATTGTCTCAAGTCACAGAattag                                                                    <  2:308300/248‑2 (MQ=255)
                                                                                                                                                                     cAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGc                                                     >  1:189021/1‑251 (MQ=255)
                                                                                                                                                                                   caccacCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGttt                                       <  1:346431/251‑1 (MQ=255)
                                                                                                                                                                                      ccccGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGaa                                    <  2:96757/249‑1 (MQ=255)
                                                                                                                                                                                                        agGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGc                  >  1:314068/1‑251 (MQ=255)
                                                                                                                                                                                                             cggacATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATAGTCCGTTGACATAATCGACTTGTAAACCAAATGGGAAAGATTTAGGTTTACAAGTCTACAACGAATTAAAAACAAAAAACACGTTTTGGAGAAGCCCCATGGATCACGGCCCACGCGGGAAATGGTTGCAAGTCACAGAATTATTTGAAAAACCGTTGCAGGATCTGCTGTTTGACGCACAGAAGGTGAATCGcacgc              >  2:243539/5‑246 (MQ=255)
                                                                                                                                                                                                              ggTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCAt            <  1:344376/251‑1 (MQ=255)
                                                                                                                                                                                                                        gtgtATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGGTGTCAAAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGGGCATCGCCAGCATTTCGATCCTc  <  2:140898/251‑1 (MQ=255)
                                                                                                                                                                                                                            aTGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGAt      <  1:74711/243‑1 (MQ=255)
                                                                                                                                                                                                                            aTGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGAAGCAGGTGCATCGCAAGCATTTCGAt      >  2:74711/1‑243 (MQ=255)
                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                             
CAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAT‑CAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTC  >  NC_000913/808990‑809455

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑