Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 314,258 G→A 9.3% intergenic (‑453/‑99) ykgP ← / → eaeH pseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913314,2580GA 9.3% 228.9 / 9.3 75intergenic (‑453/‑99)ykgP/eaeHpseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base A (5/2);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATT  >  NC_000913/314115‑314402
                                                                                                                                               |                                                                                                                                                
aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                             >  3:201567/1‑149 (MQ=255)
aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                             <  1:589346/149‑1 (MQ=255)
 tttAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAg                                                                                                                                            >  1:484044/1‑149 (MQ=255)
  ttAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAgg                                                                                                                                           <  4:648798/149‑1 (MQ=255)
     aTGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGaaa                                                                                                                                        <  2:484044/149‑1 (MQ=255)
       gACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATc                                                                                                                                      <  3:318288/149‑1 (MQ=255)
       gACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATc                                                                                                                                      <  3:481607/149‑1 (MQ=255)
          ttGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                             <  1:102342/139‑1 (MQ=255)
          ttGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                             >  2:102342/1‑139 (MQ=255)
           tGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGa                                                                                                                                  >  4:632292/1‑149 (MQ=255)
              acaAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAAt                                                                                                                               >  2:221890/1‑149 (MQ=255)
               caAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATc                                                                                                                              >  1:370994/1‑149 (MQ=255)
               caAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATc                                                                                                                              <  3:742997/149‑1 (MQ=255)
                aaaaaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcc                                                                                                                             >  3:338205/1‑149 (MQ=255)
                   aaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGTAAAGGAAATCGGGAAATCCCCg                                                                                                                          <  1:470181/149‑1 (MQ=255)
                    aGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGATATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAACCGGGGATTTCCCga                                                                                                                         <  4:66486/149‑2 (MQ=255)
                       tAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGttt                                                                                                                      <  1:221890/149‑1 (MQ=255)
                       tAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaaa                                                                                                                                        >  1:219105/1‑131 (MQ=255)
                       tAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaaa                                                                                                                                        <  2:219105/131‑1 (MQ=255)
                          aTTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTc                                                                                                                   <  2:543752/149‑1 (MQ=255)
                           ttttCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCt                                                                                                                  >  3:339810/1‑149 (MQ=255)
                                  ccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGc                                                                                                           >  1:497049/1‑149 (MQ=255)
                                        ataaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGc                                                                                                     <  3:677233/149‑1 (MQ=255)
                                             aaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatta                                                                                                >  2:321343/1‑149 (MQ=255)
                                               tataATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCattatt                                                                                              <  1:446299/149‑1 (MQ=255)
                                               tataATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCattatt                                                                                              >  2:291811/1‑149 (MQ=255)
                                               tataATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCattatt                                                                                              >  1:229208/1‑149 (MQ=255)
                                                ataATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTAttt                                                                                             <  4:46065/149‑1 (MQ=255)
                                                  aaTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtgt                                                                                           <  1:659356/149‑1 (MQ=255)
                                                  aaTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtgt                                                                                           >  2:571718/1‑149 (MQ=255)
                                                  aaTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCAATGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGAGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtgt                                                                                           <  3:10344/149‑1 (MQ=255)
                                                   aTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtgtg                                                                                          >  3:9258/1‑149 (MQ=255)
                                                   aTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtgtg                                                                                          >  3:345951/1‑149 (MQ=255)
                                                      ccGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCt                                                                                       >  2:563406/1‑149 (MQ=255)
                                                                  accTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaaa                                                                                                                                        >  1:413325/1‑88 (MQ=255)
                                                                  accTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaaa                                                                                                                                        <  2:413325/88‑1 (MQ=255)
                                                                          aGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGt                                                                   <  1:97902/149‑1 (MQ=255)
                                                                           gTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGtt                                                                  <  4:288409/149‑1 (MQ=255)
                                                                            tGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtc                                                                 >  2:676386/1‑149 (MQ=255)
                                                                              gCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcat                                                               <  1:321343/149‑1 (MQ=255)
                                                                                  taGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAgg                                                           >  2:285852/1‑149 (MQ=255)
                                                                                  taGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAgg                                                           >  3:403074/1‑149 (MQ=255)
                                                                                         tGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAt                                                    >  1:211492/1‑149 (MQ=255)
                                                                                         tGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAt                                                    <  4:345951/149‑1 (MQ=255)
                                                                                           gATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtct                                                  <  4:203780/149‑1 (MQ=255)
                                                                                            atatatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtctc                                                 <  1:571718/149‑1 (MQ=255)
                                                                                               tatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                              <  3:252721/149‑1 (MQ=255)
                                                                                                       ggCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACAtt                                      >  1:119455/1‑149 (MQ=255)
                                                                                                          ggTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATtata                                   >  3:592918/1‑149 (MQ=255)
                                                                                                           gTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaa                                  <  2:84871/149‑1 (MQ=255)
                                                                                                             gaaCTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaaaa                                <  4:97170/147‑1 (MQ=255)
                                                                                                             gCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaaaa                                <  3:265181/149‑1 (MQ=255)
                                                                                                                cTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAg                             <  1:557304/149‑1 (MQ=255)
                                                                                                                 tttCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                            >  3:291155/1‑149 (MQ=255)
                                                                                                                 tttCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                            >  1:491359/1‑149 (MQ=255)
                                                                                                                 tttCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                            >  1:576386/1‑149 (MQ=255)
                                                                                                                     ttACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAt                        >  1:125367/1‑149 (MQ=255)
                                                                                                                     ttACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAt                        >  2:210168/1‑149 (MQ=255)
                                                                                                                        cATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATaaa                     >  4:237140/1‑149 (MQ=255)
                                                                                                                         aTAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAaac                    >  3:75723/1‑149 (MQ=255)
                                                                                                                         aTAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAaac                    >  4:231852/1‑149 (MQ=255)
                                                                                                                           aaaGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAaacaa                  >  2:44391/1‑149 (MQ=255)
                                                                                                                               gTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCAc              <  4:201567/149‑1 (MQ=255)
                                                                                                                                 aTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGa            >  3:230580/1‑149 (MQ=255)
                                                                                                                                  tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGt                                                                   <  4:354643/93‑1 (MQ=255)
                                                                                                                                  tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGt                                                                   >  3:354643/1‑93 (MQ=255)
                                                                                                                                  tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTcac                                         <  4:125209/119‑1 (MQ=255)
                                                                                                                                  tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTcac                                         >  3:125209/1‑119 (MQ=255)
                                                                                                                                  tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAt           >  2:648387/1‑149 (MQ=255)
                                                                                                                                   ttCCTTTTCCGCAGAAAGGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCCTTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATAAGGAAATTATCTCAATGTCACATTATAAAAAAGGTCATAAACAACCACGAtt          >  2:646898/1‑149 (MQ=255)
                                                                                                                                   ttCCTTTTCCGCAGAAAAGGAAATCGGGAAATCACCGGTTTTTCTGACAAGCAGACGCCCTTATTCGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAtt          >  2:288965/1‑149 (MQ=255)
                                                                                                                                      cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCg       >  3:407877/1‑149 (MQ=255)
                                                                                                                                      cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACAACGATTTCg       >  3:24802/1‑149 (MQ=255)
                                                                                                                                       ttttCTGCGGAAAAGCAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGt      <  3:571528/149‑1 (MQ=255)
                                                                                                                                           cTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTAtt  <  4:403074/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
ATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATT  >  NC_000913/314115‑314402

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