Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,449,050 T→A 6.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 6.6% 254.1 / 5.3 90intergenic (+32/+26)feaB/tynAphenylacetaldehyde dehydrogenase/tyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base T (49/35);  new base A (2/4);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

TTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCACATGTGTG  >  NC_000913/1448906‑1449199
                                                                                                                                                |                                                                                                                                                     
ttGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGc                                                                                                                                                  >  2:1499648/1‑150 (MQ=255)
  gACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTTGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCg                                                                                                                                                <  1:460250/150‑1 (MQ=255)
   aCGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGTGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGGGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGt                                                                                                                                               <  2:312212/150‑1 (MQ=255)
    cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtt                                                                                                                                              >  1:1220049/1‑150 (MQ=255)
      cTAACTTCCCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGtttt                                                                                                                                            <  1:445008/150‑1 (MQ=255)
       tAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttttt                                                                                                                                           <  1:1470413/150‑1 (MQ=255)
       tAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGttttt                                                                                                                                           <  2:940781/150‑1 (MQ=255)
        aaCTTACTGTTTGGTGGGATGAAGCAGTCAGGAACGGTCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGtttttt                                                                                                                                          <  1:237846/150‑1 (MQ=255)
        aaCTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCTGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtttttt                                                                                                                                          >  1:434097/1‑150 (MQ=255)
          cTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                                                                                                                                        <  2:1192026/150‑1 (MQ=255)
          cTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGttttttac                                                                                                                                        <  1:361725/150‑3 (MQ=255)
          cTTAACGTTTGGTGGGATGAAGCAGGCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGTGGTGAAACTAAGTATGTGTGTGTACGGTATTAATCTGGTTCGCTCTTAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                                                                                                                                        <  2:286325/150‑1 (MQ=255)
             aCCGTTTGGTGGTATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCt                                                                                                                                     <  2:522192/150‑1 (MQ=255)
              ccGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTg                                                                                                                                    <  1:1458467/150‑1 (MQ=255)
              ccGGTTGGTGGGAGGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTg                                                                                                                                    <  1:1451860/150‑1 (MQ=255)
               cGTGTGTTGGGATGAAGCTGTATGGAACGGGCCGTGATTTTGGCCCGGACTGGCTGGACGGTTGGTGTGTAACTAAGTCGGTGTGTGTTAGGTTTTAATCTGGTTCGCTCATAAGTAAAAAACGGCTCCTGGTGCCGTTTTTTTGTCTGa                                                                                                                                   <  2:987947/150‑1 (MQ=255)
                       gggATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCa                                                                                                                           <  1:444786/150‑1 (MQ=255)
                          aTGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCCCCTGGTGCCGTTTTTTTGTCTGAAACAATCACtt                                                                                                                        >  1:1557673/1‑150 (MQ=255)
                             aaGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtc                                                                                                                     >  2:662080/1‑150 (MQ=255)
                             aaGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtc                                                                                                                     >  2:1717996/1‑150 (MQ=255)
                               gCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACa                                                                                                                               >  1:1017111/1‑138 (MQ=255)
                               gCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACa                                                                                                                               <  2:1017111/138‑1 (MQ=255)
                                   gcAGGACCGGGCCGTGTTTTTGGCCCCGACTGGCTGTACGGGTTGTGTGAAACTAACTCGGAGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAACCGGCACCTGGGGCCGTTTTTTTGGCTGTAACAATCACTTAtctttctt                                                                                                               <  1:1165596/149‑1 (MQ=255)
                                       gAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTACGTCGGTGTGTGTAAGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAgc                                                                                                           <  2:1714836/150‑1 (MQ=255)
                                       gAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAgc                                                                                                           >  2:958232/1‑150 (MQ=255)
                                        aaCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAgcg                                                                                                          <  1:1499648/150‑1 (MQ=255)
                                           gggCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAgcgcac                                                                                                       >  2:473422/1‑148 (MQ=255)
                                            ggCCGTGATTTTGGCCCCGACTGGCTCGACGTTAGGTGTGAAACTACCTCGGTGTGTGTACGCTATTAATCTGGTTCGCTCATAAGTAAAAAACGGAAACTTGAGCCTTTTTATTGTCTGAAACAATCACTAATCTTTCTACAGCTcccc                                                                                                      >  1:45749/1‑150 (MQ=255)
                                             gCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCt                                                                                                     >  1:433180/1‑150 (MQ=255)
                                             gCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCt                                                                                                     >  1:1326866/1‑150 (MQ=255)
                                              ccGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTa                                                                                                    >  2:17810/1‑150 (MQ=255)
                                                gTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGc                                                                                                  >  2:964154/1‑150 (MQ=255)
                                                 tGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCg                                                                                                 <  2:785677/150‑1 (MQ=255)
                                                 tGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAAAAATCACTTATATTTCTTCAGCGCCCCTAGCg                                                                                                 >  1:1644160/1‑150 (MQ=255)
                                                   aTTTTTGCCCCGACTGGCTGGACGGTTGGGGTGAAACTAAGGCGGTGTGTGTACGGTATTAATCTGGGTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGtt                                                                                               <  2:890133/150‑1 (MQ=255)
                                                   aTTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAAACTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGtt                                                                                               >  1:1165087/1‑150 (MQ=255)
                                                    ttttGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTg                                                                                              >  2:1540687/1‑150 (MQ=255)
                                                         gCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGtt                                                                                         >  1:957383/1‑150 (MQ=255)
                                                          ccccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTTTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGttt                                                                                        <  2:900322/150‑1 (MQ=255)
                                                          ccacGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGCAACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCATCGCCCCTAGCGTTGGCGttt                                                                                        >  2:122275/4‑150 (MQ=255)
                                                              gACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTc                                                                                    <  1:1717996/150‑1 (MQ=255)
                                                               aCTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCa                                                                                   >  2:1390078/1‑150 (MQ=255)
                                                                     tggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAaagaag                                                                             >  2:1024996/1‑150 (MQ=255)
                                                                        aCGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTc                                                                          >  2:1608971/1‑150 (MQ=255)
                                                                                  tgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTca                                                                >  2:419130/1‑150 (MQ=255)
                                                                                   gAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcag                                                               >  1:529316/1‑150 (MQ=255)
                                                                                    aaaCTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcagc                                                              <  2:593999/150‑1 (MQ=255)
                                                                                      aCTAAGTCGGGGTGTGTACGGTATTAATCTGGATCTCTCATAAGTAAAAAACGGCACCTGGTTCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCTTCAAAGAAGTTCCATGGCTTcagcag                                                            <  1:49047/150‑1 (MQ=255)
                                                                                         aaGTCGGTGTGTGTACGGTATTAATCTGGTTAGATCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTTTCTGAAACAATAACTTATCTTTCTTCAGCGCCCCTAGCGTTTGAGTTTCGTCAAATCATTTACATGGTTTCAGCAGAGt                                                         >  1:1033300/1‑150 (MQ=255)
                                                                                          aGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTa                                                        <  1:662080/150‑1 (MQ=255)
                                                                                             cGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGt                                                     >  2:1351475/1‑150 (MQ=255)
                                                                                                tgtgtgtACGGTTATAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTAcc                                                  <  2:1584755/150‑1 (MQ=255)
                                                                                                  tgtgtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCa                                                >  1:975307/1‑150 (MQ=255)
                                                                                                  tgtgtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCa                                                >  2:1163194/1‑150 (MQ=255)
                                                                                                    tgtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCAtt                                              >  2:1399981/1‑150 (MQ=255)
                                                                                                       aCGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATtcgg                                           <  1:1540687/150‑1 (MQ=255)
                                                                                                       aCGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATtcgg                                           <  1:1163194/150‑1 (MQ=255)
                                                                                                         ggTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATtcggtc                                         <  1:713712/150‑1 (MQ=255)
                                                                                                          gTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATtcggtcg                                        >  1:765121/1‑150 (MQ=255)
                                                                                                             ttAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCa                                     >  2:421553/1‑150 (MQ=255)
                                                                                                              tAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCAt                                    <  2:1365708/150‑1 (MQ=255)
                                                                                                              tAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCAt                                    >  1:954351/1‑150 (MQ=255)
                                                                                                              tAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTATGAAACAATAACTTATCTTTATTCAGCGCCCATAGCGTTGGCGTTTCGTCAAACAAGTTCCTTGGTTTCATCAGAGTATGTACCCATTCGGTCGCCAt                                    >  1:346184/1‑150 (MQ=255)
                                                                                                                aTCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATaa                                  >  1:498683/1‑150 (MQ=255)
                                                                                                                  cTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATc                                >  1:1137539/1‑150 (MQ=255)
                                                                                                                      ttCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGcc                            <  1:964154/150‑1 (MQ=255)
                                                                                                                         gCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCTGTCGGCATAATCGGCCAct                         >  2:675338/1‑150 (MQ=255)
                                                                                                                         gCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCAct                         >  1:1336848/1‑150 (MQ=255)
                                                                                                                         gCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCAct                         >  1:1595281/1‑150 (MQ=255)
                                                                                                                            cATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCtt                      >  1:442799/1‑150 (MQ=255)
                                                                                                                            cATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCtt                      >  1:1397467/1‑150 (MQ=255)
                                                                                                                             aTAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATTGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTc                     <  2:1165087/150‑1 (MQ=255)
                                                                                                                             aTAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTc                     <  1:1399981/150‑1 (MQ=255)
                                                                                                                               aaGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCgg                   <  2:1326866/150‑1 (MQ=255)
                                                                                                                                aGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGgc                  >  1:653496/1‑150 (MQ=255)
                                                                                                                                 caaaaaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGGTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGgcg                 >  1:484395/3‑150 (MQ=255)
                                                                                                                                  tAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGgcgc                <  1:1637647/150‑1 (MQ=255)
                                                                                                                                  tAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATAACTTTTCTTTCTTAAGCGCACCTATCGTTGTCGTTTCGTCAAGCAAGTTCCATGGTTTAAGCATAGTATGTACACATTCTGTCGGCATAATCGGCCACTCGTCGgcgc                >  2:36411/1‑150 (MQ=255)
                                                                                                                                   aaaaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGcgcg               >  1:1440549/1‑150 (MQ=255)
                                                                                                                                    aaaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCgg              >  1:1384905/1‑150 (MQ=255)
                                                                                                                                    aaaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCgg              >  1:514120/1‑150 (MQ=255)
                                                                                                                                    aaaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCgg              >  2:1470160/1‑150 (MQ=255)
                                                                                                                                     aaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCTAAGAAGTTCCATGGTTTCAGCAGAGTATGTTCCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCggg             <  2:1604422/150‑1 (MQ=255)
                                                                                                                                     aaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCggg             >  1:200167/1‑150 (MQ=255)
                                                                                                                                       aaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATAACTTATCTTTCTTCAGCGCCCCTAGCGTTGTCGTTTCGTCAAAGAAGATACATCGTTACAGAAGAGTATGTACCCATTCGGTCGGCATTATCGGCCACTCTTCGTCGCGGGcc           >  2:348474/1‑150 (MQ=255)
                                                                                                                                        aCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCca          <  2:434097/150‑1 (MQ=255)
                                                                                                                                          ggCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTAAGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCcaca        >  2:1053852/1‑150 (MQ=255)
                                                                                                                                              ccTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCACAtgtg    >  1:417069/1‑150 (MQ=255)
                                                                                                                                               cTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCTGGCCACAtgtgt   <  2:484395/150‑1 (MQ=255)
                                                                                                                                               cTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCACAtgtgt   <  2:975307/150‑1 (MQ=255)
                                                                                                                                                tGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCACAtgtgtg  <  2:1384905/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCACATGTGTG  >  NC_000913/1448906‑1449199

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: