Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 907,069 A→C 5.2% L405R (CTT→CGT)  ybjT ← putative NAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913907,0690AC 5.2% 271.2 / 7.6 77L405R (CTT→CGT) ybjTputative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base C (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 = 4.11e-01

GGGATCATCAATAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCGTGTCTGCCACAAAATATTGCCAA  >  NC_000913/906922‑907216
                                                                                                                                                   |                                                                                                                                                   
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:610733/1‑93 (MQ=255)
                                                                                   tttATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCAc                                                                                                                         >  2:584803/1‑93 (MQ=255)
                                                                                   tttATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCAc                                                                                                                         <  1:610733/93‑1 (MQ=255)
                                                                                   tttATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCAc                                                                                                                         <  1:584803/93‑1 (MQ=255)
                                                                                     tATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTc                                                               <  2:102369/149‑1 (MQ=255)
                                                                                     tATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTc                                                               <  2:198268/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                                                       <  3:345173/149‑1 (MQ=255)
                                                                                             aGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCtt                                                       >  1:160476/1‑149 (MQ=255)
                                                                                               gCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTAt                                                     <  2:455055/149‑1 (MQ=255)
                                                                                                 tAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGa                                                   <  4:201502/149‑1 (MQ=255)
                                                                                                 tAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGa                                                   >  4:125717/1‑149 (MQ=255)
                                                                                                 tAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGa                                                   >  2:774493/1‑149 (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                                                                                              >  2:436994/1‑84 (MQ=255)
                                                                                                                       cGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGt                                                                                              <  1:436994/84‑1 (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                        <  1:564622/142‑1 (MQ=255)
                                                                                                                                   gCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCgt                        >  2:564622/1‑142 (MQ=255)
                                                                                                                                      aaaTAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCAt                                  >  2:762158/1‑129 (MQ=255)
                                                                                                                                      aaaTAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCAt                                  <  1:762158/129‑1 (MQ=255)
                                                                                                                                        aTAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCGTGTCTGCCACaaa            <  3:238821/149‑1 (MQ=255)
                                                                                                                                             aacGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCg                                                                          >  1:512741/1‑82 (MQ=255)
                                                                                                                                             aacGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCg                                                                          <  4:125649/82‑1 (MQ=255)
                                                                                                                                             aacGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCg                                                                          <  2:512741/82‑1 (MQ=255)
                                                                                                                                             aacGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCg                                                                          >  3:125649/1‑82 (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)
                                                                                                                                                  aaGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCGTGTCTGCCACAAAATATTGCCaa  >  1:327130/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGGATCATCAATAACCAGTAAAACAGCCCCGGCATACCGTGCGGATGCCAGAAAGCGCGGACATCGATAGTACGATAGTCGCCTTTATCTTCCAGGCTAAAACACAGTCGTCCCAGCCCCGGCGCTTTCATGCCAAATAACAACGTAAGTTGTTTTTCCGGTTCAACGACAATCACTTTCCAGCTATCCACCGCATCGCCAGTCTGTAAATATTCGCGCTCCGGGCGGCCTTTCGCCAGCTTATGACCGATCGCGCGGTCCATCAACGCCCGTGTCTGCCACAAAATATTGCCAA  >  NC_000913/906922‑907216

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: