Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,760,053 C→T 6.2% L157L (CTG→CTA sufS ← cysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,760,0530CT 6.2% 262.6 / 4.7 81L157L (CTG→CTAsufScysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (40/36);  new base T (2/3);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

TCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACG  >  NC_000913/1759908‑1760197
                                                                                                                                                 |                                                                                                                                                
tCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGCGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCgtaggtag                                                                                                                                                >  8:315704/1‑148 (MQ=255)
tCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAg                                                                                                                                                <  1:199536/148‑1 (MQ=255)
tCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAg                                                                                                                                                <  6:405195/148‑1 (MQ=255)
          aGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGt                                                                                                                                      <  5:153022/148‑1 (MQ=255)
               cATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGTAGCGTCTCCAGTTGCaa                                                                                                                                 <  7:315704/148‑1 (MQ=255)
                  ccaccaGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACg                                                                                                                               >  4:170271/1‑147 (MQ=255)
                   caccaGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTa                                                                                                                             >  6:229700/1‑148 (MQ=255)
                    accaGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTAc                                                                                                                            >  8:217349/1‑148 (MQ=255)
                     ccaGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTAcc                                                                                                                           >  4:382731/1‑148 (MQ=255)
                      caGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCa                                                                                                                          <  8:199892/148‑1 (MQ=255)
                       aGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGt                                                                                                                              >  8:455184/1‑143 (MQ=255)
                       aGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGt                                                                                                                              <  7:455184/143‑1 (MQ=255)
                        gCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATc                                                                                                                        <  7:246921/148‑1 (MQ=255)
                             tttGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGAt                                                                                                                   >  4:23193/1‑148 (MQ=255)
                              ttGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGa                                                                                                                    <  8:398362/146‑1 (MQ=255)
                                  gcCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGc                                                                                                              >  6:15295/1‑148 (MQ=255)
                                     aTGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCgg                                                                                                            >  2:185764/1‑147 (MQ=255)
                                             gcgcAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGAt                                                                                                         >  4:251074/1‑142 (MQ=255)
                                             gcgcAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGAt                                                                                                         <  3:251074/142‑1 (MQ=255)
                                                     gTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTg                                                                                           <  1:233525/148‑1 (MQ=255)
                                                       gATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTgcg                                                                                         >  2:430704/1‑148 (MQ=255)
                                                           aTTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCa                                                                                      <  5:15295/147‑1 (MQ=255)
                                                             ttCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAAcg                                                                                   >  1:239962/1‑148 (MQ=255)
                                                                    aGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGAACCATCGGGATTGAGCGGGATCACACGCAGATCTGCGCCAACGCGTgcg                                                                             >  5:113480/1‑147 (MQ=255)
                                                                      tGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCa                                                                           >  4:6938/1‑147 (MQ=255)
                                                                       ggATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGTAGCGTCTCCAGTTGCAACGTACCATCg                                                                                                                       >  1:404603/1‑102 (MQ=255)
                                                                       ggATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGTAGCGTCTCCAGTTGCAACGTACCATCg                                                                                                                       <  2:404603/102‑1 (MQ=255)
                                                                         aTTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCACAAAGc                                                                       <  3:170271/148‑1 (MQ=255)
                                                                          ttttCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCa                                                                      <  2:239962/148‑1 (MQ=255)
                                                                          ttttCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCa                                                                      >  6:253428/1‑148 (MQ=255)
                                                                                 tgCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTg                                                                  <  5:253428/145‑1 (MQ=255)
                                                                                 tgCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGcc                                                                <  7:217349/147‑1 (MQ=255)
                                                                                   ccAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGTATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAgg                                                             <  3:6938/148‑1 (MQ=255)
                                                                                   ccAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAg                                                              <  3:382731/147‑1 (MQ=255)
                                                                                    cAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTACTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCAGCGCCAACGCGTGCGCAAAGCATCTGCCAgg                                                             >  8:130923/1‑147 (MQ=255)
                                                                                      aGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGaa                                                          >  1:384882/1‑148 (MQ=255)
                                                                                      aGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGaa                                                          >  1:333134/1‑148 (MQ=255)
                                                                                        gctgTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGTAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACa                                                        <  1:106548/145‑1 (MQ=255)
                                                                                           gTTGGAGACATGAGTAATTGCCAGCAGGCGAGATTTCTCATCAAACATCGTAGGCAGCGTCTCCAGTTGCAACGTACCAACGCGATTGAGCGGGATCAAACGCAGCTCTGCTCCAACGCGTGCGCAAAGCATCTGCCACGGAACAATg                                                     >  2:3190/1‑148 (MQ=255)
                                                                                            ttGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCCCACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGCACAATGt                                                    >  5:192656/1‑148 (MQ=255)
                                                                                             tGGAGACATGAGTAATTGCCAGAAGGCGAGTTTTCTCATCAAACAGCGTAGGTAGCGTCTCCAGTTGCAACGTACCATAGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCCTCCTCAAATCATCTGCCAGGGAACaa                                                       >  2:106548/1‑144 (MQ=255)
                                                                                               gagaCATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAg                                                 <  1:430704/148‑1 (MQ=255)
                                                                                                   cATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCgtg                                             <  1:185764/148‑1 (MQ=255)
                                                                                                     tGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCa                                                               >  6:72924/1‑128 (MQ=255)
                                                                                                     tGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCa                                                               <  5:72924/128‑1 (MQ=255)
                                                                                                     tGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCgtggt                                           >  2:16848/1‑148 (MQ=255)
                                                                                                      gAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCggg                                                                                                                     >  4:218604/1‑73 (MQ=255)
                                                                                                      gAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCggg                                                                                                                     <  3:218604/73‑1 (MQ=255)
                                                                                                      gAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCgtggtg                                          >  7:197120/1‑148 (MQ=255)
                                                                                                         tAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTc                                       >  6:139513/1‑148 (MQ=255)
                                                                                                              gCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATc                                   <  8:197120/147‑1 (MQ=255)
                                                                                                                 agcagGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATctgac                               >  7:59950/1‑148 (MQ=255)
                                                                                                                  gcagGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCAt                                                                     <  5:146456/109‑1 (MQ=255)
                                                                                                                  gcagGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCAt                                                                     >  6:146456/1‑109 (MQ=255)
                                                                                                                  gcagGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTCGCgtggt                                           >  2:262424/1‑135 (MQ=255)
                                                                                                                  gcagGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATctgact                              >  5:250707/1‑148 (MQ=255)
                                                                                                                  gcagGCGAGGTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCgtggt                                           <  1:262424/135‑1 (MQ=255)
                                                                                                                   cagGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATctgact                              >  8:17569/1‑147 (MQ=255)
                                                                                                                   cagGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGAGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATctga                                >  5:167256/1‑145 (MQ=255)
                                                                                                                     ggCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCTAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACtgat                           <  3:77450/148‑1 (MQ=255)
                                                                                                                     ggCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCTAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACtgat                           >  4:77450/1‑148 (MQ=255)
                                                                                                                       cGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCgtg                                             <  5:261486/128‑1 (MQ=255)
                                                                                                                       cGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCgtg                                             >  6:261486/1‑128 (MQ=255)
                                                                                                                        gAGTTTTCTCATCAATCAGGGTATGCAGGGTCTCCAGTTGCAACGTACCATCGGGTTTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCTTTGTGCTCCATCTGACtgatgat                        <  3:398264/148‑1 (MQ=255)
                                                                                                                         aGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACtgatgat                        >  2:246470/1‑147 (MQ=255)
                                                                                                                           ttttCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTgatgatgat                     >  7:410682/1‑148 (MQ=255)
                                                                                                                           ttttCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTgatgatgat                     <  2:333134/148‑1 (MQ=255)
                                                                                                                              tctcATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGt                   <  5:229700/147‑1 (MQ=255)
                                                                                                                              tctcATCAAACAGCGTAGGCAGCGTCTCCAGGTGCAACGTACCATCGGGATTGCGCGGGCTCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGt                   <  3:23193/147‑1 (MQ=255)
                                                                                                                                 catcaAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTATc               >  2:257119/1‑148 (MQ=255)
                                                                                                                                   tcaAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTATcgc             >  4:72164/1‑148 (MQ=255)
                                                                                                                                    caAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTa                 >  3:112442/1‑143 (MQ=255)
                                                                                                                                    caAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTa                 <  4:112442/143‑1 (MQ=255)
                                                                                                                                    caAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTAt                <  6:271064/144‑1 (MQ=255)
                                                                                                                                    caAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTAt                >  5:271064/1‑144 (MQ=255)
                                                                                                                                           cGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTATCGCccgcccgc     <  7:17569/148‑1 (MQ=255)
                                                                                                                                            gTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAgg                                                             >  1:255003/1‑91 (MQ=255)
                                                                                                                                            gTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAgg                                                             <  2:255003/91‑1 (MQ=255)
                                                                                                                                               ggCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACg  <  1:257119/147‑1 (MQ=255)
                                                                                                                                               ggCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACg  <  3:72164/147‑1 (MQ=255)
                                                                                                                                                gCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTATcgcccgcc        <  2:62072/140‑1 (MQ=255)
                                                                                                                                                gCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTATcgcccgcc        >  1:62072/1‑140 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCCAGGGAACAATGTTAGCGTGGTGCTCCATCTGACTGATGATGATGTTATCGCCCGCCCGCACG  >  NC_000913/1759908‑1760197

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