Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 691,390 G→C 5.5% A132G (GCC→GGC)  ybeX ← putative ion transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913691,3900GC 5.5% 314.7 / 7.6 91A132G (GCC→GGC) ybeXputative ion transport
Reads supporting (aligned to +/- strand):  ref base G (38/48);  new base C (2/3);  total (40/51)
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.37e-01

GGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTT  >  NC_000913/691245‑691538
                                                                                                                                                 |                                                                                                                                                    
ggTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCtt                                                                                                                                                   <  1:230059/149‑1 (MQ=255)
  tCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCa                                                                                                                                                 >  3:88031/1‑149 (MQ=255)
  tCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCa                                                                                                                                                 >  4:108942/1‑149 (MQ=255)
    aCCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGc                                                                                                                                               >  3:228658/1‑149 (MQ=255)
    aCCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGc                                                                                                                                               >  4:464331/1‑149 (MQ=255)
      cAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCAt                                                                                                                                             <  2:284229/149‑1 (MQ=255)
      cAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCAt                                                                                                                                             >  4:216857/1‑149 (MQ=255)
      cAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCAt                                                                                                                                             >  1:364485/1‑149 (MQ=255)
       aGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtg                                                                                                                                            <  3:108942/149‑2 (MQ=255)
       aGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtg                                                                                                                                            <  3:54329/149‑2 (MQ=255)
            ggAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgc                                                                                                                                       <  3:375905/149‑1 (MQ=255)
                 ccccACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATa                                                                                                                                  >  3:47862/1‑149 (MQ=255)
                   ccACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATaaa                                                                                                                                >  2:477443/1‑149 (MQ=255)
                   ccACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATaaa                                                                                                                                <  1:411184/149‑1 (MQ=255)
                    cACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGCCGGCTTCAGCATCGCTGCGCATAAAc                                                                                                                               <  4:252861/149‑1 (MQ=255)
                    cACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAAc                                                                                                                               >  3:403054/1‑149 (MQ=255)
                        gAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgca                                                                                                                           <  4:93116/149‑1 (MQ=255)
                        gAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgca                                                                                                                           <  4:433680/149‑1 (MQ=255)
                        gAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGgca                                                                                                                           <  2:221318/149‑1 (MQ=255)
                                 aaTAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCtt                                                                                                                  >  4:375193/1‑149 (MQ=255)
                                 aaTAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCtt                                                                                                                  >  1:49312/1‑149 (MQ=255)
                                  aTAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCttt                                                                                                                 >  4:63053/1‑149 (MQ=255)
                                  aTAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCttt                                                                                                                 <  2:187207/149‑1 (MQ=255)
                                    aaCGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCg                                                                                                               >  4:259624/1‑149 (MQ=255)
                                     aCGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGc                                                                                                              <  3:216857/149‑1 (MQ=255)
                                      cGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGcc                                                                                                             <  3:63053/149‑1 (MQ=255)
                                       gATCGCCATGTTGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCa                                                                                                            >  3:416704/1‑149 (MQ=255)
                                          cGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCa                                                                                                         <  4:355482/149‑1 (MQ=255)
                                           gCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAg                                                                                                        >  3:80592/1‑149 (MQ=255)
                                             cATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGaa                                                                                                      <  2:44033/149‑1 (MQ=255)
                                                gtgGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATcc                                                                                                   <  1:477443/149‑1 (MQ=255)
                                                 tgGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATccc                                                                                                  <  4:88031/149‑1 (MQ=255)
                                                  ggTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCt                                                                                                 <  3:178670/149‑1 (MQ=255)
                                                            agagCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGAt                                                                                       <  4:429855/149‑1 (MQ=255)
                                                                cGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCttt                                                                                   <  3:464331/149‑1 (MQ=255)
                                                                cGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCttt                                                                                   >  1:169647/1‑149 (MQ=255)
                                                                cGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCttt                                                                                   <  4:47862/149‑1 (MQ=255)
                                                                 gAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGc                                                                                                                                                     <  1:96682/82‑1 (MQ=255)
                                                                  aaaCTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGt                                                                                 <  3:259624/149‑1 (MQ=255)
                                                                  aaaCTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGt                                                                                 <  3:375193/149‑1 (MQ=255)
                                                                  aaaCTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGt                                                                                 >  4:397191/1‑149 (MQ=255)
                                                                          tCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTa                                                                         <  2:364485/149‑1 (MQ=255)
                                                                            aGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAAt                                                                       >  2:463200/1‑149 (MQ=255)
                                                                                tCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCAcc                                                                   >  4:290358/1‑149 (MQ=255)
                                                                                tCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCAcc                                                                   <  2:348487/149‑1 (MQ=255)
                                                                                   ggTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAAc                                                                                                                               >  4:261294/1‑86 (MQ=255)
                                                                                   ggTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAAc                                                                                                                               <  3:261294/86‑1 (MQ=255)
                                                                                           cgcTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTGAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGa                                                        <  4:416704/149‑1 (MQ=255)
                                                                                             cTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCACGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGt                                                      <  4:415962/149‑1 (MQ=255)
                                                                                              ttACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAAc                                                            <  4:227286/142‑1 (MQ=255)
                                                                                              ttACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAAc                                                            >  3:227286/1‑142 (MQ=255)
                                                                                              ttACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTg                                                     <  3:221751/149‑1 (MQ=255)
                                                                                                 cTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGc                                                  <  2:356616/149‑1 (MQ=255)
                                                                                                 cTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTCAGTGGGc                                                  <  2:257524/149‑1 (MQ=255)
                                                                                                  tttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCa                                                                                                         >  4:37038/1‑93 (MQ=255)
                                                                                                  tttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCa                                                                                                         <  3:37038/93‑1 (MQ=255)
                                                                                                  tttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGt                                                          >  1:263860/1‑140 (MQ=255)
                                                                                                  tttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGt                                                          <  2:263860/140‑1 (MQ=255)
                                                                                                       ggAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTc                                            <  1:141690/149‑1 (MQ=255)
                                                                                                       ggAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTc                                            >  4:243558/1‑149 (MQ=255)
                                                                                                         aaCGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTc                                                                                                                    >  1:31705/1‑75 (MQ=255)
                                                                                                         aaCGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTc                                                                                                                    <  2:31705/75‑1 (MQ=255)
                                                                                                         aaCGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCga                                          <  2:133018/149‑1 (MQ=255)
                                                                                                         aaCGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCga                                          <  1:108653/149‑1 (MQ=255)
                                                                                                           cGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatg                                        <  4:403054/149‑1 (MQ=255)
                                                                                                              cAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCTCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGatgatg                                     <  4:54483/149‑1 (MQ=255)
                                                                                                               aaCCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCtt                                                                              <  1:446873/107‑1 (MQ=255)
                                                                                                               aaCCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCtt                                                                              >  2:446873/1‑107 (MQ=255)
                                                                                                                    ccTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCa                               >  4:473151/1‑149 (MQ=255)
                                                                                                                       gACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatga                                       >  1:121614/1‑138 (MQ=255)
                                                                                                                       gACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatga                                       <  2:121614/138‑1 (MQ=255)
                                                                                                                            aaCACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTc                       <  3:473151/149‑1 (MQ=255)
                                                                                                                             acacTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCg                      >  3:44088/1‑149 (MQ=255)
                                                                                                                               acTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTc                    <  4:178901/149‑1 (MQ=255)
                                                                                                                                cTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTcc                   >  2:404328/1‑149 (MQ=255)
                                                                                                                                cTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTcc                   >  4:20984/1‑149 (MQ=255)
                                                                                                                                     tCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGt              >  2:258254/1‑149 (MQ=255)
                                                                                                                                     tCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGt              >  3:106747/1‑149 (MQ=255)
                                                                                                                                      ccATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTc             <  1:463200/149‑1 (MQ=255)
                                                                                                                                       caagCTGAAGGCTTCAGCATCTCTCCGCAAAAACGGCAGCAAGTCTTTCGCCATCCGAATCCCTTCAATGTGATCTTTGACTTCGCTAAACACCGGGAAACGTGATTGGGCGTACTCGAAGAGCACATCAAGACATTCGTCCAGCGTCt            >  4:296114/4‑149 (MQ=255)
                                                                                                                                       cATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCt            >  1:363790/1‑149 (MQ=255)
                                                                                                                                       cATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCt            <  4:44088/149‑1 (MQ=255)
                                                                                                                                        aTGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTg           <  4:228658/149‑1 (MQ=255)
                                                                                                                                        aTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTggg                                                   <  4:40343/109‑1 (MQ=255)
                                                                                                                                        aTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTggg                                                   >  3:40343/1‑109 (MQ=255)
                                                                                                                                         tGCTGAAGGCTTCAGCATCGCTGCGCATAGACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGt                                                          <  1:397943/101‑1 (MQ=255)
                                                                                                                                         tGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTCTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGt                                                          >  2:397943/1‑101 (MQ=255)
                                                                                                                                              aaGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCg     >  1:469909/1‑149 (MQ=255)
                                                                                                                                              aaGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCg     >  1:303776/1‑149 (MQ=255)
                                                                                                                                               aGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGt    >  2:222712/1‑149 (MQ=255)
                                                                                                                                               aGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGt    <  2:303776/149‑1 (MQ=255)
                                                                                                                                                 gCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGttt  >  1:195168/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTT  >  NC_000913/691245‑691538

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