Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,449,050 T→A 5.6% intergenic (+32/+26) feaB → / ← tynA phenylacetaldehyde dehydrogenase/tyramine oxidase, copper‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,449,0500TA 5.6% 240.6 / 2.9 81intergenic (+32/+26)feaB/tynAphenylacetaldehyde dehydrogenase/tyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base T (37/39);  new base A (2/3);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCAC  >  NC_000913/1448901‑1449192
                                                                                                                                                     |                                                                                                                                              
cTTAATTGACGCTAACTTACCGTTTGGTGGGATGCAGCAGTCAGGAACGGGCCGTGATGTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCt                                                                                                                                                <  2:1597742/150‑1 (MQ=255)
  tAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTgg                                                                                                                                              >  1:507336/1‑150 (MQ=255)
    aTTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTg                                                                                                                                            >  2:737052/1‑150 (MQ=255)
         cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtt                                                                                                                                       >  1:149649/1‑150 (MQ=255)
           cTAACTTCCCGTTTGGTGGGATGAAGCAGGCTGTAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGGTGGGGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtttt                                                                                                                                     <  1:1572202/150‑1 (MQ=255)
            tAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCACCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttttt                                                                                                                                    <  2:1694569/150‑1 (MQ=255)
             aaCTTTACGTTTGGTGGGATGAAGCAGTCAGGAACGGGGCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtttttt                                                                                                                                   <  2:441189/150‑1 (MQ=255)
             aaCTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGtttttt                                                                                                                                   <  1:1461403/150‑1 (MQ=255)
               cTTACCGTTTGGTGGGATGAAGCAGTCAGGTACGGGCCGTGTTTTTTGCCCCGACTGGCTGGACGGTGGGTGTGAAACTTAGTCGGTGTGTGTTCGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGttttttac                                                                                                                                 <  2:1599685/150‑3 (MQ=255)
                ttACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGttttttact                                                                                                                                <  1:697172/150‑4 (MQ=255)
                  aCCGTTTGGTGGGTTGAAGCAGTCTTGAACGGGCCGTGATTTTGGCCCCGACTGTCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCt                                                                                                                              <  2:801799/150‑1 (MQ=255)
                   ccGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTg                                                                                                                             <  2:403975/150‑1 (MQ=255)
                     gTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCAGTCTGGCTTGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCAAAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGaa                                                                                                                           <  2:1158105/150‑1 (MQ=255)
                       ttGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAAc                                                                                                                         <  1:769120/150‑1 (MQ=255)
                       ttGGGGGGATGAAGCAGTCAGGAACGGGCCGGGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCTCCTGGTGCCGTTTTTTTGTCTGAAAc                                                                                                                         <  1:968437/150‑1 (MQ=255)
                              gATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACt                                                                                                                  >  1:1575475/1‑150 (MQ=255)
                                tGAAGCAGTCAGGAACGGGCCTTGATTTTGGCCCCGACTGGCTGGACTGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTa                                                                                                                <  2:1055339/150‑1 (MQ=255)
                                tGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTa                                                                                                                <  2:1633644/150‑1 (MQ=255)
                                tGAACCCGTCAGGAACGGGCCGTGATTTTGGCCCCGTCTGGCGGGTCGGTTGGTGTTTATCTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTGATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTa                                                                                                                <  2:1519864/150‑1 (MQ=255)
                                 gAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAt                                                                                                               >  2:157563/1‑150 (MQ=255)
                                     cAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtcttt                                                                                                           <  1:1266958/150‑1 (MQ=255)
                                     cAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGAACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtcttt                                                                                                           >  1:1276886/1‑150 (MQ=255)
                                      aGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGATCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctttc                                                                                                          <  1:130025/150‑1 (MQ=255)
                                        gcAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGATGGGGGGATACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctttctt                                                                                                        <  2:1775882/149‑1 (MQ=255)
                                         cAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTttcttc                                                                                                       >  2:338860/1‑150 (MQ=255)
                                              aCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGTTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAgcgc                                                                                                  >  2:300858/1‑150 (MQ=255)
                                                gggCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGccc                                                                                                <  2:1455044/150‑1 (MQ=255)
                                                        aTTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGtt                                                                                        <  1:1439557/150‑1 (MQ=255)
                                                         ctttGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTg                                                                                       >  2:46289/2‑150 (MQ=255)
                                                         ctttGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTg                                                                                       >  2:51907/2‑150 (MQ=255)
                                                           ttGGCCCCGCCTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGc                                                                                     >  2:1102755/1‑150 (MQ=255)
                                                           ttGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGc                                                                                     <  1:625259/150‑1 (MQ=255)
                                                            tGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCg                                                                                    <  1:1578609/150‑1 (MQ=255)
                                                            tGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCg                                                                                    <  2:657541/150‑1 (MQ=255)
                                                             ggCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGt                                                                                   <  2:42657/150‑1 (MQ=255)
                                                              gCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGtt                                                                                  >  1:51758/1‑150 (MQ=255)
                                                               ccccGTCTGGCGGGACGGTTTGTGTGAAACTAAGTCGGTGGGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAACACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGttt                                                                                 <  2:229819/150‑1 (MQ=255)
                                                               ccccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGttt                                                                                 <  2:149649/150‑1 (MQ=255)
                                                                 ccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCg                                                                               <  2:507336/150‑1 (MQ=255)
                                                                 ccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCg                                                                               >  2:477492/1‑150 (MQ=255)
                                                                 ccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCg                                                                               <  1:117883/150‑1 (MQ=255)
                                                                    aCTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCa                                                                            <  2:831542/150‑1 (MQ=255)
                                                                     ctggctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaa                                                                           >  1:1016882/1‑150 (MQ=255)
                                                                      tggctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaaa                                                                          <  1:850994/150‑1 (MQ=255)
                                                                      tggctggACGGGTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaaa                                                                          <  2:802532/150‑1 (MQ=255)
                                                                       ggctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTATCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAaag                                                                         >  1:1667326/1‑150 (MQ=255)
                                                                        gctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGAGTTGTCGTTTCGTCAaaga                                                                        >  1:242025/1‑150 (MQ=255)
                                                                          tggACGGTTGGTGTGAAACTAAGTCGGGGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAaagaag                                                                      <  1:541601/150‑1 (MQ=255)
                                                                           ggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGt                                                                     >  2:1615028/1‑150 (MQ=255)
                                                                           ggACGGTTGGTGGGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGt                                                                     <  1:1810808/150‑1 (MQ=255)
                                                                                 ttGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATg                                                               <  2:534829/150‑1 (MQ=255)
                                                                                 ttGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCGCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATg                                                               <  1:937148/150‑1 (MQ=255)
                                                                                    gtgtGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGtt                                                            >  2:514854/1‑150 (MQ=255)
                                                                                         aaaCTAAGTAGGAGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcagc                                                       <  1:1723532/150‑1 (MQ=255)
                                                                                           aCTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAAAGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcagcag                                                     <  1:942853/150‑1 (MQ=255)
                                                                                            cTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCaga                                                    <  2:698779/150‑1 (MQ=255)
                                                                                            cTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTATTCAGCGCCCCTTGCGTTGGCGTTTAGTCTAAGAAGTTACATGGTTTcagcagt                                                    >  1:1806462/1‑149 (MQ=255)
                                                                                                   ggTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTa                                             <  2:1070324/150‑1 (MQ=255)
                                                                                                   ggTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTa                                             >  1:758746/1‑150 (MQ=255)
                                                                                                   ggTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCATCAGCGCCCCTAGCGTTGGCGTTTAGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTa                                             >  1:1566875/1‑150 (MQ=255)
                                                                                                          gtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGt                                                  <  2:178446/138‑1 (MQ=255)
                                                                                                          gtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGt                                                  >  1:178446/1‑138 (MQ=255)
                                                                                                                tATTAATCTGGTTCGCTCATAAGTAAAAAACGTCACCTGGTGCAGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATtcggtcgg                                >  1:642333/1‑150 (MQ=255)
                                                                                                                   tAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCAt                             >  2:66929/1‑150 (MQ=255)
                                                                                                                      tCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAAt                          >  2:692814/1‑150 (MQ=255)
                                                                                                                       cTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATc                         <  2:1413686/150‑1 (MQ=255)
                                                                                                                           ttCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGcc                     <  1:737052/150‑1 (MQ=255)
                                                                                                                              gCTCATAAGTAAAAAACGGCACCTGGTGCATTTTTTTTGTCTGAAACAATCACTTTTATTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTAGACATAATCGGCCAca                  >  2:450280/1‑149 (MQ=255)
                                                                                                                               cTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCActc                 >  2:1687247/1‑150 (MQ=255)
                                                                                                                               cTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAAACGGCCActc                 >  2:578567/1‑150 (MQ=255)
                                                                                                                                 cATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCtt               <  1:406546/150‑1 (MQ=255)
                                                                                                                                   tAAGTAAAAAACGGCACCTGGTGCAGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTAAAAGAAGTTCCATGGTTTCAGAAGAGTATGTACACATTCGGTCGCCATAAACGGCCACTCTTCg             >  2:448456/1‑150 (MQ=255)
                                                                                                                                   tAAGTAAAAAACGGCACCTGGTGACGTTTTTTTGTCTGAAACAATCAATTATCTTTCTTAAGCGCCCATAGCGTTGGCGTTTAGTCAAAGAAGTTCCCTGGTTTCCGCAGACTATGTACCCCTTCGGTCGGCATAATCTGCCACTCTTCg             >  2:259819/1‑150 (MQ=255)
                                                                                                                                       tAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTAATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTAAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCTCTCTTCGgcgc         >  2:1172748/1‑150 (MQ=255)
                                                                                                                                       aaaaaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGGTGGCGTTTCGTCAAAGTAGTTCCATGGTTTCAGCAGAGTATCTACCCATTCGGTCGGCATATTCGGCCACTCTTCGgcgc         >  1:112872/2‑150 (MQ=255)
                                                                                                                                         aaaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCgg       >  2:1225062/1‑150 (MQ=255)
                                                                                                                                         aaaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCgg       >  2:714610/1‑150 (MQ=255)
                                                                                                                                          aaaaCGGCGCCTGGTGCCGTTTTTTTGTCTGAAACAATCTCTTATCTTTCTTCAGCGCCCCTAGCGTAGGCGTTTCGTCAAAGAAGTTACATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGGATAATCGGCCACTCTTCGGCGCggg      >  2:147236/1‑150 (MQ=255)
                                                                                                                                          aaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCggg      >  2:1658652/1‑150 (MQ=255)
                                                                                                                                           aaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGc     >  1:555603/1‑150 (MQ=255)
                                                                                                                                           aaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGc     >  1:546115/1‑150 (MQ=255)
                                                                                                                                              cGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGGTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCcac  >  2:1790487/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
CTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCAC  >  NC_000913/1448901‑1449192

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: