Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,463,240 T→A 6.0% S372S (TCT→TCA paaK → phenylacetate‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,463,2400TA 6.0% 316.8 / 6.2 99S372S (TCT→TCApaaKphenylacetate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (46/46);  new base A (2/4);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

CGCATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAG  >  NC_000913/1463091‑1463388
                                                                                                                                                     |                                                                                                                                                    
cGCATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCt                                                                                                                                                      <  2:455269/150‑1 (MQ=255)
 gCATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtg                                                                                                                                                     >  2:625781/1‑150 (MQ=255)
  cATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtgt                                                                                                                                                    <  1:53318/150‑1 (MQ=255)
  cATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtgt                                                                                                                                                    >  2:53318/1‑150 (MQ=255)
   aTCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtgtg                                                                                                                                                   <  1:1046832/150‑1 (MQ=255)
   aTCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtgtg                                                                                                                                                   <  1:998259/150‑1 (MQ=255)
   aTCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCtgtg                                                                                                                                                   <  2:145858/150‑1 (MQ=255)
    tCAGCTGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGa                                                                                                                                                  <  1:926735/150‑1 (MQ=255)
    tCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGa                                                                                                                                                  <  2:715352/150‑1 (MQ=255)
      aGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGaaa                                                                                                                                                <  2:78264/150‑1 (MQ=255)
       gCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCGCAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAg                                                                                                                                               >  1:518882/1‑150 (MQ=255)
        cGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGt                                                                                                                                              <  1:673056/150‑1 (MQ=255)
         ggACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTg                                                                                                                                             <  2:655246/150‑1 (MQ=255)
          gACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTgg                                                                                                                                            >  1:757223/1‑150 (MQ=255)
           aCGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGa                                                                                                                                           <  1:1616634/150‑1 (MQ=255)
            cGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAg                                                                                                                                          <  1:1587447/150‑1 (MQ=255)
            cGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAg                                                                                                                                          <  1:1239733/150‑1 (MQ=255)
             gCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCACAGAAAGTgaatc                                                                                                                                         <  1:1183882/150‑5 (MQ=255)
              cAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGATGTGAACCGCCGGGGGCATCTTGATTCACTTTCACAGAAAGTgaatca                                                                                                                                        <  1:1437828/150‑6 (MQ=255)
               aGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTg                                                                                                                                       <  1:651686/150‑1 (MQ=255)
               aGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTg                                                                                                                                       <  2:1526479/150‑1 (MQ=255)
               aGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTg                                                                                                                                       <  2:1173481/150‑1 (MQ=255)
               aGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCACAGAAAGTgaatcaa                                                                                                                                       <  1:487225/150‑7 (MQ=255)
                  gATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTgaaa                                                                                                                                    <  2:377299/150‑1 (MQ=255)
                   aTGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTgaaag                                                                                                                                   <  2:801863/150‑1 (MQ=255)
                      aCATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTg                                                                                                                                       >  1:706114/1‑143 (MQ=255)
                      aCATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTg                                                                                                                                       <  2:706114/143‑1 (MQ=255)
                       cATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGaaagaaag                                                                                                                               <  1:344581/150‑1 (MQ=255)
                        aTGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGt                                                                                                                              <  2:350470/150‑1 (MQ=255)
                         tGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTa                                                                                                                             >  2:1372642/1‑150 (MQ=255)
                                 attCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACa                                                                                                                     >  2:1313945/1‑150 (MQ=255)
                                  ttCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACAt                                                                                                                    <  2:1236579/150‑1 (MQ=255)
                                     gcgGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGa                                                                                                                 >  1:1529596/1‑150 (MQ=255)
                                     gcgGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGa                                                                                                                 >  1:244838/1‑150 (MQ=255)
                                         tgtAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCa                                                                                                             >  1:1161980/1‑150 (MQ=255)
                                           tAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATg                                                                                                           <  1:641015/150‑1 (MQ=255)
                                              aTGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGagc                                                                                                        <  1:1372642/150‑1 (MQ=255)
                                                gTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGagcag                                                                                                      <  1:317581/150‑1 (MQ=255)
                                                  cTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGCGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCg                                                                                                    <  1:454197/150‑1 (MQ=255)
                                                  cTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCg                                                                                                    >  2:1066285/1‑150 (MQ=255)
                                                  cTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCg                                                                                                    >  1:1174750/1‑150 (MQ=255)
                                                   tttCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGt                                                                                                   <  1:107287/150‑1 (MQ=255)
                                                     tCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTg                                                                                                 <  1:499765/150‑1 (MQ=255)
                                                      ccGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGc                                                                                                >  2:1441335/1‑150 (MQ=255)
                                                       cGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCACCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGcc                                                                                               >  2:1188937/1‑150 (MQ=255)
                                                        gTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCa                                                                                              >  2:546234/1‑150 (MQ=255)
                                                              aCTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATg                                                                                        <  2:931266/150‑1 (MQ=255)
                                                              aCTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATg                                                                                        <  2:1188546/150‑1 (MQ=255)
                                                               cTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGc                                                                                       >  2:577656/1‑150 (MQ=255)
                                                                   aagaagAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGccacc                                                                                   >  1:1374321/1‑150 (MQ=255)
                                                                    agaagaGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGccacca                                                                                  >  1:1110387/1‑150 (MQ=255)
                                                                    agaagaGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGccacca                                                                                  <  1:557327/150‑1 (MQ=255)
                                                                       agagATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCt                                                                               >  2:1582976/1‑150 (MQ=255)
                                                                       agagATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCt                                                                               >  2:1528443/1‑150 (MQ=255)
                                                                              gTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCAt                                                                        <  2:393082/150‑1 (MQ=255)
                                                                               tCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATc                                                                       >  2:1390051/1‑150 (MQ=255)
                                                                                   aTTCGAACATTTATCGCCGCATTACCAACTGGATGTGAATAGCCGGGGTCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAATATAGTAGCCTGACATTGTCGCTTGTGCAGCGTTGCCAGGTATTCCACCAGCTGCTACATCGGAt                                                                   <  2:1434811/150‑1 (MQ=255)
                                                                                      cGAACATTTATCGCCGCATTTCCAACTGGAGGTGAACCGCCGGGGCCTCATTATTGACTTTCGGTGAAAGTGGTGTTGAAATGAAGGAGCCGGCCATTTCCCCAAGAGCAGCGTTGGCCGGTATGCCAACAGCGTCGGCCTCGGATTTAg                                                               >  2:1434909/1‑150 (MQ=255)
                                                                                      cGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTaa                                                                >  2:66247/1‑150 (MQ=255)
                                                                                      cGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTaa                                                                >  2:1375849/1‑150 (MQ=255)
                                                                                        aaCATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGt                                                              >  1:725496/1‑150 (MQ=255)
                                                                                          cATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCCGCTGCGCCATCGGATTAAGTCg                                                            >  2:377051/1‑150 (MQ=255)
                                                                                                tCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGAtggtg                                                      <  2:1239652/150‑1 (MQ=255)
                                                                                                tCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGAtggtg                                                      >  1:660445/1‑150 (MQ=255)
                                                                                                  gccgcATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTggg                                                    <  1:1441335/150‑1 (MQ=255)
                                                                                                    cgcATTACCAACTGGATGTGAACCGCCGGCGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGa                                                  <  2:808021/150‑1 (MQ=255)
                                                                                                    cgcATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGa                                                  >  1:1243571/1‑150 (MQ=255)
                                                                                                      cATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGTCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGAtc                                                >  2:298137/1‑150 (MQ=255)
                                                                                                      cATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGAtc                                                <  2:1366717/150‑1 (MQ=255)
                                                                                                         tACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATctct                                             <  2:48424/150‑1 (MQ=255)
                                                                                                         tACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATctct                                             <  2:42401/150‑1 (MQ=255)
                                                                                                          aCCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTa                                            >  1:487513/1‑150 (MQ=255)
                                                                                                          aCCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTa                                            >  2:764262/1‑150 (MQ=255)
                                                                                                           ccAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTAc                                           >  1:470277/1‑150 (MQ=255)
                                                                                                             aaCTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCg                                         <  2:641124/150‑1 (MQ=255)
                                                                                                                tGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGAtg                                      <  2:1186017/150‑1 (MQ=255)
                                                                                                                 ggAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGAtgt                                     >  2:1024046/1‑150 (MQ=255)
                                                                                                                   aGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGtga                                   >  2:453456/1‑150 (MQ=255)
                                                                                                                   aGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGtga                                   <  1:546234/150‑1 (MQ=255)
                                                                                                                   aGGTGAACCGCCGGGGGCATCTTGATTCACGTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGtga                                   >  2:1583642/1‑150 (MQ=255)
                                                                                                                    ggTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGtgat                                  >  1:920029/1‑150 (MQ=255)
                                                                                                                    ggTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGtgat                                  <  1:940095/150‑1 (MQ=255)
                                                                                                                    ggTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGtgat                                  >  1:413103/1‑150 (MQ=255)
                                                                                                                      tGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGtgatga                                >  2:503901/1‑150 (MQ=255)
                                                                                                                      tGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGtgatga                                >  2:1242361/1‑150 (MQ=255)
                                                                                                                         aCCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCg                             <  2:413103/150‑1 (MQ=255)
                                                                                                                         aCCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCg                             >  2:1614257/1‑150 (MQ=255)
                                                                                                                         aCCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCg                             <  2:518882/150‑1 (MQ=255)
                                                                                                                            gccgGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTa                          <  1:1390051/150‑1 (MQ=255)
                                                                                                                              cgGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAAc                        <  1:577656/150‑1 (MQ=255)
                                                                                                                                ggggCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACtg                      <  1:338997/150‑1 (MQ=255)
                                                                                                                                    cATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGc                  >  1:17837/1‑150 (MQ=255)
                                                                                                                                    cATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGc                  <  2:1529596/150‑1 (MQ=255)
                                                                                                                                      tCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAg                <  1:625781/150‑1 (MQ=255)
                                                                                                                                          actccACTTTCACAGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAg                <  2:63024/142‑1 (MQ=255)
                                                                                                                                          actccACTTTCACAGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAg                >  1:63024/5‑146 (MQ=255)
                                                                                                                                          actccACTTTCACAGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATc            >  2:690485/5‑150 (MQ=255)
                                                                                                                                           aTTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATcc           >  1:1302697/1‑150 (MQ=255)
                                                                                                                                              cACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCcgc        <  1:453456/150‑1 (MQ=255)
                                                                                                                                              cACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCcgc        >  2:1073793/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCATCAGCGGACGCAGCGATGACATGTTGATTATTCGCGGTGTAAATGTCTTTCCGTCACAACTGGAAGAAGAGATTGTCAAATTCGAACATTTATCGCCGCATTACCAACTGGAGGTGAACCGCCGGGGGCATCTTGATTCACTTTCTGTGAAAGTGGAGTTGAAAGAAAGTAGCCTGACATTGACGCATGAGCAGCGTTGCCAGGTATGCCACCAGCTGCGCCATCGGATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAG  >  NC_000913/1463091‑1463388

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