Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 314,258 G→A 5.9% 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 5.9% 409.6 / 8.1 128intergenic (‑453/‑99)ykgP/eaeHpseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified
Reads supporting (aligned to +/- strand):  ref base G (65/55);  new base A (2/6);  total (67/61)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 131 total aligned reads displayed.

ATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGT  >  NC_000913/314111‑314406
                                                                                                                                                   |                                                                                                                                                    
aTCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCgg                                                                                                                                                     >  4:941400/1‑149 (MQ=255)
aTCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCGGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCgg                                                                                                                                                     >  1:134697/1‑149 (MQ=255)
 tCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGa                                                                                                                                                    >  4:320141/1‑149 (MQ=255)
 tCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGa                                                                                                                                                    >  1:54625/1‑149 (MQ=255)
    aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGTGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                 >  2:901553/1‑149 (MQ=255)
    aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                 <  3:110927/149‑1 (MQ=255)
    aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                 <  4:72327/149‑1 (MQ=255)
    aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                 <  1:159602/149‑1 (MQ=255)
     tttAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAg                                                                                                                                                >  3:786593/1‑149 (MQ=255)
       tAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGa                                                                                                                                              <  3:373972/149‑1 (MQ=255)
        aaTGACTTGCACAAAAAGCTAAATTTTCCCCCATATATAACAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGTAGAAAAGGaa                                                                                                                                             <  4:704067/149‑1 (MQ=255)
        aaTGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAATAATATAATCACGCGCCCAACCACCTGATGAGTGGCTATAGGGACTGGATTTATTAGGTGGCGGTGCACTCTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaa                                                                                                                                             <  3:672923/149‑1 (MQ=255)
           gACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGAAAta                                                                                                                                          <  1:343096/149‑2 (MQ=255)
            aCTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCg                                                                                                                                         >  2:70206/1‑149 (MQ=255)
             cTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                        <  3:320141/149‑1 (MQ=255)
              ttGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCggg                                                                                                                                       >  1:671852/1‑149 (MQ=255)
              ttGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                        <  4:735178/148‑1 (MQ=255)
               tGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGa                                                                                                                                      <  4:970200/149‑1 (MQ=255)
                    aaaaaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCAAACCACCTGATGAGTGGCTATAGGCACTGGATATACTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcc                                                                                                                                 >  4:871734/1‑149 (MQ=255)
                      aaaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcccc                                                                                                                               <  2:666087/149‑1 (MQ=255)
                       aaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCg                                                                                                                              <  1:60928/149‑1 (MQ=255)
                          cTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGtt                                                                                                                           >  1:341507/1‑149 (MQ=255)
                              aTTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGATATCGGGAAATCCCCGGTTTTTc                                                                                                                       <  2:451757/149‑1 (MQ=255)
                                    ccccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACaa                                                                                                                 <  4:553403/149‑1 (MQ=255)
                                      ccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGc                                                                                                               >  2:1052388/1‑149 (MQ=255)
                                        ataaataaAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAg                                                                                                             <  2:777926/149‑1 (MQ=255)
                                        ataaataaAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAg                                                                                                             >  1:952909/1‑149 (MQ=255)
                                             taaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGcc                                                                                                        >  2:523600/1‑149 (MQ=255)
                                             taaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGcc                                                                                                        >  3:923976/1‑149 (MQ=255)
                                                aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatt                                                                                                     >  1:435388/1‑149 (MQ=255)
                                                 aaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatta                                                                                                    >  3:860911/1‑149 (MQ=255)
                                                  atatAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCattat                                                                                                   >  1:344232/1‑149 (MQ=255)
                                                  atatAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCCGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCattat                                                                                                   <  1:523600/149‑1 (MQ=255)
                                                   tataATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCattatt                                                                                                  >  2:766260/1‑149 (MQ=255)
                                                      aaTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtgt                                                                                               <  1:637096/149‑1 (MQ=255)
                                                          ccGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                 <  2:273249/95‑1 (MQ=255)
                                                          ccGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCTGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCt                                                                                           <  2:290624/149‑1 (MQ=255)
                                                          ccGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCt                                                                                           >  2:735279/1‑149 (MQ=255)
                                                          ccGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCt                                                                                           >  2:797349/1‑149 (MQ=255)
                                                            gcgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGc                                                                                         <  4:656741/149‑1 (MQ=255)
                                                                    ccaccTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTc                                                                                 <  2:1069454/149‑1 (MQ=255)
                                                                      accTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGt                                                                               <  1:215472/149‑1 (MQ=255)
                                                                       ccTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGtt                                                                              <  3:1040296/149‑1 (MQ=255)
                                                                         tgatgaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTaa                                                                            >  3:70770/1‑149 (MQ=255)
                                                                         tgatgaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTaa                                                                            >  4:5355/1‑149 (MQ=255)
                                                                         tgatgaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTaa                                                                            >  1:351428/1‑149 (MQ=255)
                                                                             gaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCggg                                                                                                                                       >  3:895208/1‑86 (MQ=255)
                                                                             gaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCggg                                                                                                                                       <  4:895208/86‑1 (MQ=255)
                                                                             gaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCg                                                                        >  2:956262/1‑149 (MQ=255)
                                                                              aGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGt                                                                       <  4:464785/149‑1 (MQ=255)
                                                                               gTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGtt                                                                      >  1:49401/1‑149 (MQ=255)
                                                                                tGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtc                                                                     >  3:792505/1‑149 (MQ=255)
                                                                                    tataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatca                                                                 <  4:639976/149‑1 (MQ=255)
                                                                                    tataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatca                                                                 <  3:767948/149‑1 (MQ=255)
                                                                                    tataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatca                                                                 <  3:838533/149‑1 (MQ=255)
                                                                                     ataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAg                                                                <  3:23804/149‑1 (MQ=255)
                                                                                        ggCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGaa                                                             >  4:975571/1‑149 (MQ=255)
                                                                                         gCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGAc                                                                                                           <  4:30898/102‑1 (MQ=255)
                                                                                           aCTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAAtt                                                          >  4:159882/1‑149 (MQ=255)
                                                                                           aCTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAAtt                                                          >  2:422355/1‑149 (MQ=255)
                                                                                              ggATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtc                                                       >  3:989705/1‑149 (MQ=255)
                                                                                              ggATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtc                                                       >  4:324691/1‑149 (MQ=255)
                                                                                               gATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtct                                                      >  2:365136/1‑149 (MQ=255)
                                                                                                  atatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCaa                                                   <  3:5355/149‑1 (MQ=255)
                                                                                                  atatTAGGTGGCGGTGCACCTTCTTACATAAAGGTATTTCCTTTTCTGCAGAAAAGGAAAACGGCAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCaa                                                   <  2:789694/149‑1 (MQ=255)
                                                                                                   tatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                  <  1:956262/149‑1 (MQ=255)
                                                                                                      tAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTc                                               <  4:923976/149‑1 (MQ=255)
                                                                                                      tAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTc                                               <  2:134697/149‑1 (MQ=255)
                                                                                                               gTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaa                                      >  3:611052/1‑149 (MQ=255)
                                                                                                                tGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaaa                                     >  4:1003388/1‑149 (MQ=255)
                                                                                                                 gCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaaaa                                    <  2:341507/149‑1 (MQ=255)
                                                                                                                   aCTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACa                                  <  3:159882/149‑1 (MQ=255)
                                                                                                                   aCTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACa                                  <  1:735279/149‑1 (MQ=255)
                                                                                                                   aCTGTCTTAGATAAAGGTATTTACTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCATTGTCACATTATAAAACa                                  >  4:890461/1‑149 (MQ=255)
                                                                                                                    cTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAg                                 <  4:965297/149‑1 (MQ=255)
                                                                                                                     ttttttACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCAAATTATAAAACAgg                                >  4:920727/1‑149 (MQ=255)
                                                                                                                       tCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                                <  1:1034008/147‑1 (MQ=255)
                                                                                                                       tCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                                >  2:1034008/1‑147 (MQ=255)
                                                                                                                       tCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTc                              <  2:478356/149‑1 (MQ=255)
                                                                                                                            cATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTTACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATaaa                         <  2:944366/149‑1 (MQ=255)
                                                                                                                             aTAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGTAAGTCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCATTGTCACATTATAAAACAGTTCATAaac                        >  2:389366/1‑149 (MQ=255)
                                                                                                                             aTAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTAttt                                                                                                 >  1:1053327/1‑76 (MQ=255)
                                                                                                                             aTAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAaac                        >  3:589584/1‑149 (MQ=255)
                                                                                                                                 aGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAAcc                    >  1:201230/1‑149 (MQ=255)
                                                                                                                                 aGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAAcc                    >  2:462256/1‑149 (MQ=255)
                                                                                                                                   gTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACa                                  >  4:507377/1‑133 (MQ=255)
                                                                                                                                   gTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACa                                  <  3:507377/133‑1 (MQ=255)
                                                                                                                                    gaTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACg                 >  1:1002583/2‑149 (MQ=255)
                                                                                                                                      tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGt                               >  1:909652/1‑133 (MQ=255)
                                                                                                                                      tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGt                               <  2:909652/133‑1 (MQ=255)
                                                                                                                                      tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAt               >  2:45853/1‑149 (MQ=255)
                                                                                                                                         ccTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTc            >  1:93474/1‑149 (MQ=255)
                                                                                                                                          cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTAACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCg           <  1:901553/149‑1 (MQ=255)
                                                                                                                                           ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtctc                                                     >  1:276638/1‑106 (MQ=255)
                                                                                                                                           ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGt          >  4:568010/1‑149 (MQ=255)
                                                                                                                                           ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGt          >  3:45901/1‑149 (MQ=255)
                                                                                                                                               cTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTAtt      <  1:326559/149‑1 (MQ=255)
                                                                                                                                                tGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTc     <  1:783661/149‑1 (MQ=255)
                                                                                                                                                  cGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAg   <  1:766260/149‑1 (MQ=255)
                                                                                                                                                   ggAAAAAGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGt  <  2:318361/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGT  >  NC_000913/314111‑314406

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: