Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,999,231 T→A 6.9% T84S (ACT→TCT)  dcyD ← D‑cysteine desulfhydrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,999,2310TA 6.9% 264.6 / 6.6 86T84S (ACT→TCT) dcyDD‑cysteine desulfhydrase
Reads supporting (aligned to +/- strand):  ref base T (39/41);  new base A (4/2);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTT  >  NC_000913/1999088‑1999380
                                                                                                                                               |                                                                                                                                                     
tGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGTAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCg                                                                                                                                                 <  1:184402/150‑1 (MQ=255)
     cGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACAt                                                                                                                                            >  1:1331952/1‑150 (MQ=255)
      gCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtcg                                                                                                                                           <  1:516734/150‑5 (MQ=255)
       cacaTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATgg                                                                                                                                          >  1:867413/1‑150 (MQ=255)
         caTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGtt                                                                                                                                        <  1:361494/150‑1 (MQ=255)
         caTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGtt                                                                                                                                        <  2:890572/150‑1 (MQ=255)
         caTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtcgctg                                                                                                                                        <  1:405183/150‑8 (MQ=255)
          aTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTa                                                                                                                                       <  2:664547/150‑1 (MQ=255)
              cAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACt                                                                                                                                   <  2:588193/150‑1 (MQ=255)
               aaTCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTg                                                                                                                                  >  1:587292/1‑150 (MQ=255)
                  cTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGAt                                                                                                                               <  2:20554/150‑1 (MQ=255)
                  cTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGGTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGAGGCAGTCTGCCGCACATGGTTAGACTGGAt                                                                                                                               <  2:693568/150‑1 (MQ=255)
                         ttGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCg                                                                                                                        <  1:178673/150‑1 (MQ=255)
                         ttGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCg                                                                                                                        <  2:1376130/150‑1 (MQ=255)
                          tGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCgg                                                                                                                       >  1:549983/1‑150 (MQ=255)
                           gAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGcccc                                                                                                                         <  2:207055/147‑1 (MQ=255)
                           gAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGcccc                                                                                                                         >  1:207055/1‑147 (MQ=255)
                           gAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGc                                                                                                                      <  2:1190879/150‑1 (MQ=255)
                           gAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGc                                                                                                                      >  2:579195/1‑150 (MQ=255)
                             aCAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAg                                                                                                                    <  2:1411071/150‑1 (MQ=255)
                              cAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGt                                                                                                                   >  1:1006363/1‑150 (MQ=255)
                                gATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTaa                                                                                                                 <  1:986949/150‑1 (MQ=255)
                                  tCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATc                                                                                                               <  2:867413/150‑1 (MQ=255)
                                   ccAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCa                                                                                                              >  1:1378827/1‑150 (MQ=255)
                                     aGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGc                                                                                                            <  1:987344/150‑1 (MQ=255)
                                          caaGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATc                                                                                                       >  1:33390/1‑150 (MQ=255)
                                            aGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCgg                                                                                                     >  2:13613/1‑150 (MQ=255)
                                                      gTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACg                                                                                           >  2:1123537/1‑150 (MQ=255)
                                                       ttGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGc                                                                                          <  2:992310/150‑1 (MQ=255)
                                                        tGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCa                                                                                         >  2:1436773/1‑150 (MQ=255)
                                                          gTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGt                                                                                                                   <  2:1056605/122‑1 (MQ=255)
                                                          gTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGt                                                                                                                   >  1:1056605/1‑122 (MQ=255)
                                                           ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCagag                                                                                      <  1:352184/150‑1 (MQ=255)
                                                              aaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCagag                                                                                      <  2:1076828/147‑1 (MQ=255)
                                                              aaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCagag                                                                                      >  1:1076828/1‑147 (MQ=255)
                                                              aaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCAt                                                                                   >  2:1520824/1‑150 (MQ=255)
                                                              aaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCAt                                                                                   >  1:918837/1‑150 (MQ=255)
                                                               aTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATc                                                                                  >  1:367864/1‑150 (MQ=255)
                                                                tAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCt                                                                                 <  2:6705/150‑1 (MQ=255)
                                                                tAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCACCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCt                                                                                 >  1:204446/1‑150 (MQ=255)
                                                                      tCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCcgcg                                                                           >  1:1108926/1‑150 (MQ=255)
                                                                       cTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCga                                                                          >  2:1099176/1‑150 (MQ=255)
                                                                        tGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCgag                                                                         <  2:261821/150‑1 (MQ=255)
                                                                         gCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCgaga                                                                        >  2:1030439/1‑150 (MQ=255)
                                                                            gTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAAt                                                                     <  2:1378827/150‑1 (MQ=255)
                                                                              tgtgCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTc                                                                   <  1:907001/150‑1 (MQ=255)
                                                                               gtgCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTcc                                                                  <  2:587292/150‑1 (MQ=255)
                                                                                 gCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTATTCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg                                                                >  1:229753/1‑150 (MQ=255)
                                                                                 gCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg                                                                >  2:1159916/1‑150 (MQ=255)
                                                                                   cAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCt                                                              <  1:701283/150‑1 (MQ=255)
                                                                                   cAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCt                                                              <  2:1006363/150‑1 (MQ=255)
                                                                                    aaTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCtt                                                             <  1:1361178/150‑1 (MQ=255)
                                                                                      tAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTAc                                                           <  1:1030439/150‑1 (MQ=255)
                                                                                        ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCcgc                                                                            >  1:975122/1‑131 (MQ=255)
                                                                                        ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCcgc                                                                            <  2:975122/131‑1 (MQ=255)
                                                                                           ttttCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAAt                                                      <  1:579195/150‑1 (MQ=255)
                                                                                             ttCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAAttt                                                    >  1:793154/1‑150 (MQ=255)
                                                                                              tCCAGCAGCGACACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTa                                                   >  2:957294/1‑150 (MQ=255)
                                                                                               ccAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAt                                                  >  1:393140/1‑150 (MQ=255)
                                                                                                 agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTg                                                >  1:336142/1‑150 (MQ=255)
                                                                                                           cGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTg                                      <  1:342846/150‑1 (MQ=255)
                                                                                                             cAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGcc                                    <  1:733738/150‑1 (MQ=255)
                                                                                                               gTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCAt                                  >  2:19787/1‑150 (MQ=255)
                                                                                                                 gCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATgg                                >  2:1132445/1‑150 (MQ=255)
                                                                                                                    gACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGTATCTCCCCGGCAGTAATCAGCGTATGGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTTCCGCCCATTGCCATggggg                             <  2:681626/150‑1 (MQ=255)
                                                                                                                    gACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATggggg                             >  1:1185867/1‑150 (MQ=255)
                                                                                                                       cGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGa                          <  1:19787/150‑1 (MQ=255)
                                                                                                                              ggCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATc                   >  1:1458158/1‑150 (MQ=255)
                                                                                                                                cAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATccc                 >  1:1162906/1‑150 (MQ=255)
                                                                                                                                cAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATccc                 >  1:448802/1‑150 (MQ=255)
                                                                                                                                  ccatgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCCCGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGt               >  1:1014861/6‑150 (MQ=255)
                                                                                                                                  ccatgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGt               >  1:840197/6‑150 (MQ=255)
                                                                                                                                   cGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGtt              <  2:744410/150‑1 (MQ=255)
                                                                                                                                    gACTGCGGCAGTCTGCCGCACATGGTTAGTCTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTTCGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGttt             <  2:518885/150‑1 (MQ=255)
                                                                                                                                    atgtGCGCCAGACTGCCGCACATGGTTAGACTGGATCGACCCGGCAGTAATCAGCGTATCGGCACCTTCACTCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGttt             >  2:1431661/4‑150 (MQ=255)
                                                                                                                                     tgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtg            >  2:273908/3‑150 (MQ=255)
                                                                                                                                     aCTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtg            <  2:549983/150‑1 (MQ=255)
                                                                                                                                        gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTc                      <  2:1251352/137‑1 (MQ=255)
                                                                                                                                        gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTc                      >  1:1251352/1‑137 (MQ=255)
                                                                                                                                        gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtgatg         >  1:592357/1‑150 (MQ=255)
                                                                                                                                        gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtgatg         >  1:220006/1‑150 (MQ=255)
                                                                                                                                        gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtgatg         >  1:683433/1‑150 (MQ=255)
                                                                                                                                         cGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtgatga        >  2:919582/1‑150 (MQ=255)
                                                                                                                                          ggCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGaa       <  1:13613/150‑1 (MQ=255)
                                                                                                                                            cAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGaaaa     <  2:1145983/150‑1 (MQ=255)
                                                                                                                                            cAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGaaaa     <  2:626893/150‑1 (MQ=255)
                                                                                                                                               tCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAAttt  >  1:754085/1‑150 (MQ=255)
                                                                                                                                               tCTGCCGCACATGGTAAGACTCGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTCACGCGAGAAATTCCAGCTAACGTAATTTATTGCCGCCAATTGCCATGGCGGTGACGTCATCCCGTTTGATGAAAAttt  >  2:1400738/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
TGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTT  >  NC_000913/1999088‑1999380

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: ‑