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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,999,2310TA 8.6% 541.3 / 32.9 162T84S (ACT→TCT) dcyDD‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (73/75);  new base A (8/6);  total (81/81)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

Only 100 of 166 total aligned reads displayed.

GTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATT  >  NC_000913/1999083‑1999379
                                                                                                                                                    |                                                                                                                                                    
 tCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCTAGCAACAAGCGATTGCCGTTGGCTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTc                                                                                                                                                     >  3:1490667/1‑149 (MQ=255)
         tCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGcac                                                                                                                                             >  3:1223829/1‑149 (MQ=255)
          cGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGcaca                                                                                                                                            <  4:2328/149‑1 (MQ=255)
             acaTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGGCGCACATgg                                                                                                                                         <  1:39475/149‑1 (MQ=255)
             acaTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATgg                                                                                                                                         >  4:87110/1‑149 (MQ=255)
                  tCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGa                                                                                                                                    <  3:1500250/149‑1 (MQ=255)
                   cAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGAc                                                                                                                                   <  4:148509/149‑1 (MQ=255)
                       cTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGa                                                                                                                               <  1:564699/149‑1 (MQ=255)
                         gggTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATc                                                                                                                             >  2:1453279/1‑149 (MQ=255)
                           gTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGc                                                                                                                           <  3:17480/149‑1 (MQ=255)
                            tATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGcc                                                                                                                          >  2:213154/1‑149 (MQ=255)
                             aTTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGccc                                                                                                                         <  4:639478/149‑1 (MQ=255)
                              aagAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGcccc                                                                                                                        <  4:738628/147‑1 (MQ=255)
                                 aaCAGATCCGGCTACAAGCGATTGCCGTTTGATAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACTCAGTGCAGACTGATTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTTAATCGCCCCGGc                                                                                                                     <  4:356031/149‑1 (MQ=255)
                                 aaCAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGc                                                                                                                     >  1:600591/1‑149 (MQ=255)
                                 aaCAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGc                                                                                                                     >  2:1486668/1‑149 (MQ=255)
                                  aCAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcag                                                                                                                                             <  4:464150/124‑2 (MQ=255)
                                      aTCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTaa                                                                                                                <  1:734059/149‑1 (MQ=255)
                                         cAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGTAGTCTGCCGCACATGGtt                                                                                                                                       <  1:497843/123‑1 (MQ=255)
                                         cAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGtt                                                                                                                                       >  2:497843/1‑123 (MQ=255)
                                          aGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAg                                                                                                            >  2:1672957/1‑149 (MQ=255)
                                          aGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAg                                                                                                            >  4:1049919/1‑149 (MQ=255)
                                          aGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAg                                                                                                            >  3:338198/1‑149 (MQ=255)
                                           gCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGc                                                                                                           <  2:1655024/149‑1 (MQ=255)
                                            caacaaGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCg                                                                                                          <  2:368404/149‑1 (MQ=255)
                                                aaGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATc                                                                                                      >  1:1238992/1‑149 (MQ=255)
                                                    gATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCa                                                                                                  >  1:1095429/1‑149 (MQ=255)
                                                      ttGTCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCAcc                                                                                                >  2:211482/1‑149 (MQ=255)
                                                       tGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCt                                                                                               >  2:576863/1‑149 (MQ=255)
                                                           gTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTgg                                                                                                                                >  4:1113451/1‑112 (MQ=255)
                                                           gTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTgg                                                                                                                                <  3:1113451/112‑1 (MQ=255)
                                                           gTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCAc                                                                                           <  3:268213/149‑1 (MQ=255)
                                                             tGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGc                                                                                         <  3:87110/149‑1 (MQ=255)
                                                                ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                                                                      >  4:635383/1‑149 (MQ=255)
                                                                ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                                                                      >  4:487248/1‑149 (MQ=255)
                                                                ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                                                                      <  1:1305092/149‑1 (MQ=255)
                                                                ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                                                                      >  4:1225919/1‑149 (MQ=255)
                                                                  aaaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagt                                                                                                                                            >  1:387106/1‑91 (MQ=255)
                                                                  aaaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagt                                                                                                                                            <  2:387106/93‑3 (MQ=255)
                                                                   aaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCa                                                                                   <  2:1436971/149‑1 (MQ=255)
                                                                   aaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCa                                                                                   >  3:1588979/1‑149 (MQ=255)
                                                                    aTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCAt                                                                                  >  1:223188/1‑149 (MQ=255)
                                                                    aTAGCTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCAt                                                                                  >  4:1561473/1‑149 (MQ=255)
                                                                     tAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATc                                                                                 >  1:1014787/1‑149 (MQ=255)
                                                                       gTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATc                                                                                                              <  1:1507389/118‑1 (MQ=255)
                                                                         tttCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCAc                                                                                           <  2:950592/135‑1 (MQ=255)
                                                                         tttCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCAc                                                                                           >  1:950592/1‑135 (MQ=255)
                                                                          ttCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCcg                                                                            <  2:831169/149‑1 (MQ=255)
                                                                           tCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCcgc                                                                           <  2:361873/149‑1 (MQ=255)
                                                                                     tgCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTcc                                                                 >  3:620335/1‑149 (MQ=255)
                                                                                      gCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGc                                                                                                           <  4:1682100/106‑1 (MQ=255)
                                                                                      gCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGc                                                                                                           >  3:1682100/1‑106 (MQ=255)
                                                                                       ccAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg                                                               <  1:1586794/149‑1 (MQ=255)
                                                                                         aaTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTa                                                                                                                                      >  2:908524/1‑76 (MQ=255)
                                                                                         aaTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATc                                                                                                                             >  4:1312432/1‑85 (MQ=255)
                                                                                         aaTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGAACATGGTTAGACTGGATc                                                                                                                             <  3:1312432/85‑1 (MQ=255)
                                                                                          aTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCtt                                                            <  2:105366/149‑1 (MQ=255)
                                                                                           tAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTa                                                           >  4:1622759/1‑149 (MQ=255)
                                                                                             ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACg                                                         >  3:509799/1‑149 (MQ=255)
                                                                                              gATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa                                                                >  3:1685027/1‑141 (MQ=255)
                                                                                              gATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGt                                                        >  2:1075120/1‑149 (MQ=255)
                                                                                                ttttCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTaa                                                      <  3:580639/149‑1 (MQ=255)
                                                                                                     cagcagCGCCACGCAATGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAt                                                 >  4:1616344/1‑149 (MQ=255)
                                                                                                      agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAtt                                                >  1:1122271/1‑149 (MQ=255)
                                                                                                      agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAtt                                                >  3:209620/1‑149 (MQ=255)
                                                                                                      agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAtt                                                <  3:635383/149‑1 (MQ=255)
                                                                                                       gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACAt                                                                                                                                           <  2:1529824/57‑1 (MQ=255)
                                                                                                        cagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgc                                              >  1:1525141/1‑149 (MQ=255)
                                                                                                            gcCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGTCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgccgcc                                          <  4:1143509/149‑1 (MQ=255)
                                                                                                             ccACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGccc                                         <  1:1224991/149‑1 (MQ=255)
                                                                                                                 acAGTGCAGACCGAGTTTGGGAGCGACTGCGGCAGTGTGCCGAAAATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTTCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTg                                     >  4:108092/2‑149 (MQ=255)
                                                                                                                   aGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGcc                                   >  4:1076084/1‑149 (MQ=255)
                                                                                                                     tGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCAt                                 <  2:639459/149‑1 (MQ=255)
                                                                                                                     tGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCAt                                 >  1:1154692/1‑149 (MQ=255)
                                                                                                                     tGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCAt                                 <  4:572593/149‑1 (MQ=255)
                                                                                                                        aGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATggg                              <  1:1060185/149‑1 (MQ=255)
                                                                                                                        aGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATggg                              <  1:272851/149‑1 (MQ=255)
                                                                                                                         gACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATgggg                             <  1:685330/149‑1 (MQ=255)
                                                                                                                         gACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATgggg                             <  1:903984/149‑1 (MQ=255)
                                                                                                                           ccGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGt                           >  1:117449/1‑149 (MQ=255)
                                                                                                                                   ggCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCAt                   <  1:843059/149‑1 (MQ=255)
                                                                                                                                    gCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa                                                                >  4:1417709/1‑103 (MQ=255)
                                                                                                                                     cAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATcc                 <  1:1135367/149‑1 (MQ=255)
                                                                                                                                     cAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATcc                 <  4:811675/149‑1 (MQ=255)
                                                                                                                                      aGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATccc                >  2:1598656/1‑149 (MQ=255)
                                                                                                                                      aGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATccc                >  4:168949/1‑149 (MQ=255)
                                                                                                                                       ccatgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCg               >  3:928176/6‑149 (MQ=255)
                                                                                                                                        catgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGAc                        >  3:79566/5‑139 (MQ=255)
                                                                                                                                        catgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGAc                        >  3:1084801/5‑139 (MQ=255)
                                                                                                                                        catgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGt              >  1:799479/5‑149 (MQ=255)
                                                                                                                                         gACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGtt             <  1:1452183/149‑1 (MQ=255)
                                                                                                                                           gtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtg           >  3:50453/2‑149 (MQ=255)
                                                                                                                                           cTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa                                                                <  1:314283/96‑1 (MQ=255)
                                                                                                                                            tGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtga          <  2:1095429/149‑1 (MQ=255)
                                                                                                                                            tGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtga          <  3:1049919/149‑1 (MQ=255)
                                                                                                                                              cGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtgatg        <  4:928176/149‑1 (MQ=255)
                                                                                                                                                 cAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGaaa     >  2:1317201/1‑149 (MQ=255)
                                                                                                                                                 cAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGaaa     <  3:1225919/149‑1 (MQ=255)
                                                                                                                                                  aGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGaaaa    <  2:1265986/149‑1 (MQ=255)
                                                                                                                                                    tCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAAtt  >  1:874145/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATT  >  NC_000913/1999083‑1999379

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: