Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 907,066 G→T 5.2% T406K (ACG→AAG)  ybjT ← putative NAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913907,0660GT 5.2% 277.3 / 4.4 77T406K (ACG→AAG) ybjTputative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (37/36);  new base T (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GGGATCATCAATAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCGTGTCTGCCACAAAATATTG  >  NC_000913/906922‑907212
                                                                                                                                                |                                                                                                                                                  
gggATCATCAATAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAg                                                                                                                                                <  1:357378/149‑1 (MQ=255)
  gATCATCAATAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAgtt                                                                                                                                              <  2:181091/149‑1 (MQ=255)
   atcatcAATAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAgttg                                                                                                                                             >  1:102369/1‑149 (MQ=255)
   atcatcAATAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAgttg                                                                                                                                             >  3:539366/1‑149 (MQ=255)
     catcaATAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAgttgtt                                                                                                                                           >  1:765076/1‑149 (MQ=255)
        caATAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGttttt                                                                                                                                        >  2:605403/1‑149 (MQ=255)
         aaTAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTc                                                                                                                                       <  4:539366/149‑1 (MQ=255)
          aTAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACTTACgttgttat                                                                                                                                         <  4:798469/146‑3 (MQ=255)
          aTAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACTTACgttgttat                                                                                                                                         >  3:798469/1‑144 (MQ=255)
          aTAACCAGTAAAACAGCCCCGGCATACCGTGCAGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTcc                                                                                                                                      >  4:685744/1‑149 (MQ=255)
            aaCCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTcc                                                                                                                                      <  2:364084/147‑1 (MQ=255)
            aaCCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTcc                                                                                                                                      >  1:364084/1‑147 (MQ=255)
               cAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAAAAACAACTTACgttgttattt                                                                                                                                       >  4:51293/1‑139 (MQ=255)
               cAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAAAAACAACTTACGTTGttattt                                                                                                                                       <  3:51293/143‑2 (MQ=255)
                  tAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAAc                                                                                                                              <  3:294535/149‑1 (MQ=255)
                    aaaCAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCTAGGCTAAAACACAGTAGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGa                                                                                                                            >  4:373710/1‑149 (MQ=255)
                     aaCAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGTCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGAc                                                                                                                           <  4:785618/149‑1 (MQ=255)
                     aaCAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGAc                                                                                                                           >  4:727476/1‑149 (MQ=255)
                          ccccGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCa                                                                                                                      <  2:765076/149‑1 (MQ=255)
                                 aTACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCa                                                                                                               <  4:149224/149‑1 (MQ=255)
                                                                cGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAgttgtt                                                                                                                                           >  1:518762/1‑90 (MQ=255)
                                                                cGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAgttgtt                                                                                                                                           <  2:518762/90‑1 (MQ=255)
                                                                   tAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTcg                                                                             <  3:685744/149‑1 (MQ=255)
                                                                     gTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTcgcg                                                                           >  3:201502/1‑149 (MQ=255)
                                                                      tACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCgcgc                                                                          >  4:238821/1‑149 (MQ=255)
                                                                      tACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCgcgc                                                                          >  2:25165/1‑149 (MQ=255)
                                                                      tACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACCTTCCAGCTATCCACCGCATCACCAGTCTGTAAATATTCgcgc                                                                          <  2:118799/149‑1 (MQ=255)
                                                                        cGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTc                                                                        <  1:323141/149‑1 (MQ=255)
                                                                             gTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGggc                                                                   <  2:131130/149‑1 (MQ=255)
                                                                                gCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGggcggc                                                                >  2:309141/1‑149 (MQ=255)
                                                                                   tttATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCAc                                                                                                                     >  2:584803/1‑93 (MQ=255)
                                                                                   tttATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCAc                                                                                                                     >  2:610733/1‑93 (MQ=255)
                                                                                   tttATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCAc                                                                                                                     <  1:610733/93‑1 (MQ=255)
                                                                                   tttATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCAc                                                                                                                     <  1:584803/93‑1 (MQ=255)
                                                                                     tATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTc                                                           <  2:198268/149‑1 (MQ=255)
                                                                                     tATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTc                                                           <  2:102369/149‑1 (MQ=255)
                                                                                      aTCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCg                                                          >  2:657116/1‑149 (MQ=255)
                                                                                         ttCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTc                                                                                                                                       <  2:371767/69‑1 (MQ=255)
                                                                                         ttCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTc                                                                                                                                       >  1:371767/1‑69 (MQ=255)
                                                                                            cAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCt                                                    >  1:742027/1‑149 (MQ=255)
                                                                                             aGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCtt                                                   >  1:160476/1‑149 (MQ=255)
                                                                                             aGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCtt                                                   <  3:345173/149‑1 (MQ=255)
                                                                                               gCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTAt                                                 <  2:455055/149‑1 (MQ=255)
                                                                                                 tAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGa                                               >  2:774493/1‑149 (MQ=255)
                                                                                                 tAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGa                                               >  4:125717/1‑149 (MQ=255)
                                                                                                 tAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGa                                               <  4:201502/149‑1 (MQ=255)
                                                                                                 tAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGa                                               <  3:125717/149‑1 (MQ=255)
                                                                                                     acacAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGa                                           <  4:298307/149‑1 (MQ=255)
                                                                                                     acacAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGa                                           <  1:774493/149‑1 (MQ=255)
                                                                                                       acaGTCGTCCCAGCCCCGGCGCTTTCATGCCAAGTAACAACGTAAGGTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATc                                         >  4:781584/1‑149 (MQ=255)
                                                                                                         aGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATcgc                                       >  1:750318/1‑149 (MQ=255)
                                                                                                         aGTCGTCCCAGCCCCGACGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATcgc                                       >  2:245128/1‑149 (MQ=255)
                                                                                                              tcCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCtt                                                   >  4:789409/1‑132 (MQ=255)
                                                                                                              tcCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCtt                                                   <  3:789409/132‑1 (MQ=255)
                                                                                                                ccAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTcc                                >  2:140059/1‑149 (MQ=255)
                                                                                                                  aGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCAt                              >  3:423123/1‑149 (MQ=255)
                                                                                                                       cGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGt                                                                                          <  1:436994/84‑1 (MQ=255)
                                                                                                                       cGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGt                                                                                          >  2:436994/1‑84 (MQ=255)
                                                                                                                             tttCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTc                                 >  1:141354/1‑135 (MQ=255)
                                                                                                                             tttCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTc                                 <  2:141354/135‑1 (MQ=255)
                                                                                                                             tttCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATAGCGCGGTCCATCAACGCCCgtg                   >  1:673222/1‑149 (MQ=255)
                                                                                                                                  tGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCGTGTCTGc              <  2:219701/149‑1 (MQ=255)
                                                                                                                                  tGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATAGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGa                                               >  4:553885/1‑116 (MQ=255)
                                                                                                                                  tGCCAAATAACAACGTAAGTTGTTTTTACGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGa                                               <  3:553885/116‑1 (MQ=255)
                                                                                                                                   gCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCgt                    >  2:564622/1‑142 (MQ=255)
                                                                                                                                   gCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCgt                    <  1:564622/142‑1 (MQ=255)
                                                                                                                                      aaaTAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCAt                              <  1:762158/129‑1 (MQ=255)
                                                                                                                                      aaaTAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCAt                              >  2:762158/1‑129 (MQ=255)
                                                                                                                                        aTAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCGTGTCTGCCACaaa        <  3:238821/149‑1 (MQ=255)
                                                                                                                                             aacGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCg                                                                      <  2:512741/82‑1 (MQ=255)
                                                                                                                                             aacGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCg                                                                      >  1:512741/1‑82 (MQ=255)
                                                                                                                                             aacGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCg                                                                      >  3:125649/1‑82 (MQ=255)
                                                                                                                                             aacGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCg                                                                      <  4:125649/82‑1 (MQ=255)
                                                                                                                                             aacGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATc                             >  2:714420/1‑123 (MQ=255)
                                                                                                                                             aacGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATc                             <  1:714420/123‑1 (MQ=255)
                                                                                                                                             aacGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCGTGTCTGCCACAAAATAtt   <  3:373710/149‑1 (MQ=255)
                                                                                                                                              acGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCGTGTCTGCCACAAAATATTg  >  3:648710/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GGGATCATCAATAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCGTGTCTGCCACAAAATATTG  >  NC_000913/906922‑907212

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