Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,759,548 A→G 5.9% Y326H (TAT→CAT)  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,759,5480AG 5.9% 208.4 / 3.4 68Y326H (TAT→CAT) sufScysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base G (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACC  >  NC_000913/1759405‑1759695
                                                                                                                                               |                                                                                                                                                   
gTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATagagag                                                                                                                                                >  3:332154/1‑149 (MQ=255)
   gTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTg                                                                                                                                              <  6:390753/148‑1 (MQ=255)
    tAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTga                                                                                                                                             <  3:135651/148‑1 (MQ=255)
     aaTAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGt                                                                                                                                               >  3:123481/1‑145 (MQ=255)
       tAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATAAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATc                                                                                                                                         <  3:164405/149‑1 (MQ=255)
         ggCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCg                                                                                                                                       <  4:381095/149‑1 (MQ=255)
          gCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCgg                                                                                                                                      >  4:340638/1‑149 (MQ=255)
            cATCAATGGCATTTCGCAGTGATGTCCGGTACGCACAGCAATTCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGATATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTa                                                                                                                                    <  1:227205/149‑1 (MQ=255)
            cATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTa                                                                                                                                    <  6:105755/149‑1 (MQ=255)
             aTCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTa                                                                                                                                    >  2:70766/1‑148 (MQ=255)
                  tGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGAtt                                                                                                                              >  7:192726/1‑149 (MQ=255)
                     cATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCa                                                                                                                           <  3:384686/149‑1 (MQ=255)
                      aTTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCACGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAg                                                                                                                          >  2:104607/1‑149 (MQ=255)
                         gcgcAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCtg                                                                                                                       >  6:9595/1‑149 (MQ=255)
                            cAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGa                                                                                                                    <  5:385986/149‑1 (MQ=255)
                             aGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGa                                                                                                                    >  8:164414/1‑148 (MQ=255)
                             aGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGa                                                                                                                    >  5:238363/1‑148 (MQ=255)
                              gTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATa                                                                                                                  <  4:332154/149‑1 (MQ=255)
                              gTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATa                                                                                                                  <  3:16695/149‑1 (MQ=255)
                              gTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATa                                                                                                                  <  3:194394/149‑1 (MQ=255)
                               tGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTTGGCCATAGATAGTGAGCTCCGGTACAGATTACAGCTGTGATAg                                                                                                                 >  6:213656/1‑149 (MQ=255)
                               tGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGa                                                                                                                    >  2:183205/1‑146 (MQ=255)
                               tGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGa                                                                                                                    <  1:183205/146‑1 (MQ=255)
                               tGATGTCCGGTACGCACAGCAATGCAGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGAAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGAGTCCAGCTGTGATAg                                                                                                                 >  2:293476/1‑149 (MQ=255)
                                  tGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAgatg                                                                                                              >  5:146714/1‑146 (MQ=255)
                                      cGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATaa                                                                                                          <  8:192726/149‑1 (MQ=255)
                                       ggTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATGGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAAt                                                                                                         >  6:10156/1‑149 (MQ=255)
                                       ggTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATGGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAAt                                                                                                         >  6:4559/1‑149 (MQ=255)
                                       ggTACGCAAAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCTTAGGCGTGGGGTTTACCGAGATTATAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAAt                                                                                                         <  6:223499/149‑1 (MQ=255)
                                        gTACGCACAGCAATGCCGTAATTATCGAGATAACTGCCAACATCATAGGCGTTGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATATTg                                                                                                        <  6:358695/149‑1 (MQ=255)
                                                    aTGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTg                                                                                            >  1:371104/1‑149 (MQ=255)
                                                     tGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGt                                                                                           >  1:135628/1‑149 (MQ=255)
                                                        cGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGttga                                                                                        >  1:248550/1‑147 (MQ=255)
                                                            aTTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATGGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACt                                                                                    <  5:4559/149‑1 (MQ=255)
                                                            aTTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATGGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACt                                                                                    <  5:10156/149‑1 (MQ=255)
                                                              tATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCg                                                                                  >  8:15356/1‑149 (MQ=255)
                                                                tCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCa                                                                                        <  1:300776/141‑1 (MQ=255)
                                                                tCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCa                                                                                        >  2:300776/1‑141 (MQ=255)
                                                                  gataaaaCTTCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAATCCTGTTTTTTGGGCCATAGAGAGTGATATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTa                                                                              <  4:361537/146‑1 (MQ=255)
                                                                  gagaAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATaa                                                                                                          >  2:277301/1‑121 (MQ=255)
                                                                  gagaAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATaa                                                                                                          <  1:277301/121‑1 (MQ=255)
                                                                     aaaaCTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGt                                                                           >  2:281840/1‑149 (MQ=255)
                                                                        aCTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGttat                                                                        <  2:218806/149‑1 (MQ=255)
                                                                         cTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGttatt                                                                       <  1:104607/149‑1 (MQ=255)
                                                                               aCATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATACTGTTCATACTCGGATATGTTATTAAGccc                                                                 >  4:107726/1‑149 (MQ=255)
                                                                                       ggCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGcc                                                          >  8:136641/1‑148 (MQ=255)
                                                                                       ggCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGcc                                                          >  5:368774/1‑148 (MQ=255)
                                                                                          gtggtgTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGaaa                                                      <  7:369353/149‑1 (MQ=255)
                                                                                          gtggtgTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGaaa                                                      <  1:70766/149‑1 (MQ=255)
                                                                                           tggtgTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGaaa                                                      >  2:383449/1‑148 (MQ=255)
                                                                                             gtgtTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACAt                                                   <  5:9595/149‑1 (MQ=255)
                                                                                              tgtTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCAATGTTATTAAGCCCCAGCGCCGAAACATa                                                  <  2:135628/149‑1 (MQ=255)
                                                                                                tttACCGAGATTATAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCAGAAACATACt                                                >  1:195040/1‑149 (MQ=255)
                                                                                                tttACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCACAGCGCCGAAACATACt                                                >  2:175657/1‑149 (MQ=255)
                                                                                                     cGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAgc                                           >  5:288992/1‑149 (MQ=255)
                                                                                                          ttAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCCGCGCCGAAACATACTCCAgcgccgc                                      >  2:218473/1‑149 (MQ=255)
                                                                                                               aGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCaa                                 <  4:123481/149‑1 (MQ=255)
                                                                                                               aGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCaa                                 <  2:371104/149‑1 (MQ=255)
                                                                                                                 cAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGACCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGa                               >  8:78701/1‑149 (MQ=255)
                                                                                                                  aaTAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGAc                              >  5:67218/1‑149 (MQ=255)
                                                                                                                  aaTAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGAc                              >  2:200791/1‑149 (MQ=255)
                                                                                                                      aCGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAat                          <  6:238363/149‑1 (MQ=255)
                                                                                                                        gCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatga                        <  1:281840/149‑1 (MQ=255)
                                                                                                                         ccAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATAAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgat                       <  3:340638/149‑1 (MQ=255)
                                                                                                                                     tgtgGGCAATAGAGAGTGAGATCCGGTACAGCTTCCAGCTGTGATAGCGCATAATGCATCAGATTCAGTTCATACTCGGCTATGTTATTAAGCCCGAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTAt            <  5:213656/148‑1 (MQ=255)
                                                                                                                                       tgGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTgg         <  6:67218/149‑1 (MQ=255)
                                                                                                                                       tgGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTgg         <  1:383449/149‑1 (MQ=255)
                                                                                                                                         ggCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGgt       >  8:272426/1‑149 (MQ=255)
                                                                                                                                              tAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTAcc  >  1:261026/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACC  >  NC_000913/1759405‑1759695

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: