Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 217,030 C→T 5.9% intergenic (‑27/+27) yaeF ← / ← proS putative lipoprotein/prolyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913217,0300CT 5.9% 332.8 / 8.7 101intergenic (‑27/+27)yaeF/proSputative lipoprotein/prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (47/48);  new base T (3/3);  total (50/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 = 5.46e-01

Only 100 of 103 total aligned reads displayed.

TCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAG  >  NC_000913/216891‑217176
                                                                                                                                           |                                                                                                                                                  
tCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGaaa                                                                                                                                             >  3:1135237/1‑147 (MQ=255)
tCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGaaa                                                                                                                                             <  4:1135237/147‑1 (MQ=255)
  ggCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCa                                                                                                                                         <  3:87590/149‑1 (MQ=255)
  ggCATCTACTGCTGTGGCTGATGGGTCGGGCGGGGTGATTTCAACAGTAAAAGCTGAAAGTAAAAGAAAGCTGGGGAGCAACAGGCAGCAGTACGCCGATGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCa                                                                                                                                         <  4:347677/149‑1 (MQ=255)
   gCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAg                                                                                                                                        >  4:1366151/1‑149 (MQ=255)
   gCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAg                                                                                                                                        <  4:221953/149‑1 (MQ=255)
           tgctgTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGa                                                                                                                                      >  4:193526/1‑143 (MQ=255)
           tgctgTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGa                                                                                                                                      <  3:193526/143‑1 (MQ=255)
             ctgTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                                                                              >  4:1370650/1‑149 (MQ=255)
             ctgTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                                                                              <  1:759054/149‑1 (MQ=255)
                tgGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGc                                                                                                                           <  3:1370650/149‑1 (MQ=255)
                tgGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGc                                                                                                                           >  3:1586006/1‑149 (MQ=255)
                      aTGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGc                                                                                                                     <  4:1467168/149‑1 (MQ=255)
                                ctggctgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                                                           >  4:837464/1‑149 (MQ=255)
                                  ggctgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGtt                                                                                                         <  3:837464/149‑1 (MQ=255)
                                     tgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCa                                                                                                      >  2:1192258/1‑149 (MQ=255)
                                        tatCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGt                                                                                                          <  1:1205147/142‑1 (MQ=255)
                                        tatCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGt                                                                                                          >  2:1205147/1‑142 (MQ=255)
                                          tCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGa                                                                                                 <  3:1366151/149‑1 (MQ=255)
                                            aaCGGTACAAGCCGAAAGTAAAGGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                                                           <  1:1459853/137‑1 (MQ=255)
                                            aaCGGTACAAGCCGAAAGTAAAGGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                                                           >  2:1459853/1‑137 (MQ=255)
                                            aaCGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                                                           >  2:696301/1‑137 (MQ=255)
                                            aaCGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                                                           <  1:696301/137‑1 (MQ=255)
                                                    aaGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacga                                                                                       >  4:733539/1‑149 (MQ=255)
                                                    aaGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTACCCATCTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacga                                                                                       <  1:555081/149‑1 (MQ=255)
                                                     aGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGc                                                                                                                     >  3:801140/1‑118 (MQ=255)
                                                     aGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGc                                                                                                                     <  4:801140/118‑1 (MQ=255)
                                                         gAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtt                                                                                                                                 <  1:1517524/102‑1 (MQ=255)
                                                         gAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTAAGGGCCCtt                                                                                                                                 >  2:1517524/1‑102 (MQ=255)
                                                             agtaaagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGAt                                                                                      <  2:1293122/138‑1 (MQ=255)
                                                             agtaaagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGAt                                                                                      >  1:1293122/4‑141 (MQ=255)
                                                              tAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAg                                                                                                                                        >  3:771911/1‑90 (MQ=255)
                                                              tAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAg                                                                                                                                        <  4:771911/90‑1 (MQ=255)
                                                               aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTg                                                                                                                             >  3:979314/1‑100 (MQ=255)
                                                               aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTg                                                                                                                             <  4:979314/100‑1 (MQ=255)
                                                                  agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCAc                                                                                 >  1:1179677/1‑141 (MQ=255)
                                                                  agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCAc                                                                                 <  2:1179677/141‑1 (MQ=255)
                                                                      agCTGGGGAGAAACAGGCGGCAGTACGCCATTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaatta                                                                     >  1:1491192/1‑149 (MQ=255)
                                                                       gCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                    >  4:1540374/1‑149 (MQ=255)
                                                                         tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCAt                                                                                                                         <  2:1087314/94‑1 (MQ=255)
                                                                         tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCAt                                                                                                                         >  1:1087314/1‑94 (MQ=255)
                                                                         tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGa                                                                                                 >  1:380492/1‑118 (MQ=255)
                                                                         tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGa                                                                                                 <  2:380492/118‑1 (MQ=255)
                                                                                aaaCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                                                           <  4:1586006/149‑1 (MQ=255)
                                                                                    aGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCg                                                       >  1:118756/1‑149 (MQ=255)
                                                                                       cggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg                                                    <  3:733539/149‑1 (MQ=255)
                                                                                        ggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGc                                                   <  2:118756/149‑1 (MQ=255)
                                                                                         gcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                    <  2:1822768/131‑1 (MQ=255)
                                                                                         gcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                    >  1:1822768/1‑131 (MQ=255)
                                                                                          cAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcga                                                 <  1:1192258/149‑1 (MQ=255)
                                                                                          cAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcga                                                 <  3:987924/149‑1 (MQ=255)
                                                                                                gCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAcc                                                                                                    <  4:353888/92‑1 (MQ=255)
                                                                                                gCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAcc                                                                                                    >  3:353888/1‑92 (MQ=255)
                                                                                                gCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAc                                                                                                     >  2:592746/1‑91 (MQ=255)
                                                                                                gCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAc                                                                                                     <  1:592746/91‑1 (MQ=255)
                                                                                                  cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCtt                                                                                                                  <  2:298137/76‑1 (MQ=255)
                                                                                                  cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCtt                                                                                                                  >  1:298137/1‑76 (MQ=255)
                                                                                                  cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTc                                                                            <  3:1416820/114‑1 (MQ=255)
                                                                                                  cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTc                                                                            >  4:1416820/1‑114 (MQ=255)
                                                                                                   tttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTGACTGTTTCTCGCCGTTGCGACGATAtt                                          >  1:364894/1‑147 (MQ=255)
                                                                                                   tttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTGACTGTTTCTCGCCGTTGCGACGATAtt                                          <  2:364894/147‑1 (MQ=255)
                                                                                                   tttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTta                                        >  2:519060/1‑149 (MQ=255)
                                                                                                    ttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTtat                                       <  3:1540374/149‑1 (MQ=255)
                                                                                                       gTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATAtt                                          <  1:938060/143‑1 (MQ=255)
                                                                                                       gTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATAtt                                          >  2:938060/1‑143 (MQ=255)
                                                                                                        tttATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTc                                   >  3:1340187/1‑149 (MQ=255)
                                                                                                            tCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGata                               >  4:1400938/1‑149 (MQ=255)
                                                                                                              cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGat                                                                                                <  3:214647/82‑1 (MQ=255)
                                                                                                              cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGat                                                                                                >  4:214647/1‑82 (MQ=255)
                                                                                                              cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTta                                        >  3:1431049/1‑138 (MQ=255)
                                                                                                              cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTta                                        <  4:1431049/138‑1 (MQ=255)
                                                                                                              cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtc                             <  3:1400938/149‑1 (MQ=255)
                                                                                                                          tCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgac                                                                                         >  3:567545/1‑77 (MQ=255)
                                                                                                                          tCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgac                                                                                         <  4:1660186/77‑1 (MQ=255)
                                                                                                                          tCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgac                                                                                         >  3:1660186/1‑77 (MQ=255)
                                                                                                                          tCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgac                                                                                         <  4:567545/77‑1 (MQ=255)
                                                                                                                          tCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                    <  2:273856/98‑1 (MQ=255)
                                                                                                                          tCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                    <  4:557516/98‑1 (MQ=255)
                                                                                                                          tCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                    >  1:273856/1‑98 (MQ=255)
                                                                                                                          tCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                    >  3:557516/1‑98 (MQ=255)
                                                                                                                          tCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAg                 <  4:1340187/149‑1 (MQ=255)
                                                                                                                             aTGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcgtc                          >  2:1808360/1‑137 (MQ=255)
                                                                                                                             aTGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcgtc                          <  1:1808360/137‑1 (MQ=255)
                                                                                                                             aTGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGtt              >  4:224502/1‑149 (MQ=255)
                                                                                                                                tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGt                                                                >  4:1190972/1‑96 (MQ=255)
                                                                                                                                tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGt                                                                <  3:1190972/96‑1 (MQ=255)
                                                                                                                                tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGa                  <  2:1485212/142‑1 (MQ=255)
                                                                                                                                tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGa                  >  1:1485212/1‑142 (MQ=255)
                                                                                                                                tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACg           >  2:1442849/1‑149 (MQ=255)
                                                                                                                                 aaaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACgg          >  3:1322680/1‑149 (MQ=255)
                                                                                                                                  cagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGa                                 >  4:107382/4‑125 (MQ=255)
                                                                                                                                  cagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGa                                 <  3:107382/122‑1 (MQ=255)
                                                                                                                                  cagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAg                 <  4:80066/138‑1 (MQ=255)
                                                                                                                                      aaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGcc     >  1:699546/1‑149 (MQ=255)
                                                                                                                                        gggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCg                                                       <  2:1409034/97‑1 (MQ=255)
                                                                                                                                        gggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCg                                                       >  1:1409034/1‑97 (MQ=255)
                                                                                                                                        gggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCa   <  1:1442849/149‑1 (MQ=255)
                                                                                                                                         ggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATCTTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAg  >  3:1163269/1‑149 (MQ=255)
                                                                                                                                           cccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGtt                       >  4:1695718/1‑126 (MQ=255)
                                                                                                                                           cccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGtt                       <  3:1695718/126‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
TCGGCATCTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAG  >  NC_000913/216891‑217176

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