Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,998,332 T→A 6.1% pseudogene (497/1137 nt) ygeQ ← pseudogene, glycosyl hydrolase family 15, part of T3SS PAI ETT2 remnant

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,998,3320TA 6.1% 979.8 / 7.4 331pseudogene (497/1137 nt)ygeQpseudogene, glycosyl hydrolase family 15, part of T3SS PAI ETT2 remnant
Reads supporting (aligned to +/- strand):  ref base T (127/168);  new base A (28/3);  total (160/171)
Fisher's exact test for biased strand distribution p-value = 2.21e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 334 total aligned reads displayed.

AAGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATAGTGGGTTTCT  >  NC_000913/2998183‑2998481
                                                                                                                                                     |                                                                                                                                                     
  gCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGc                                                                                                                                                     <  2:7132091/150‑1 (MQ=255)
  gCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGc                                                                                                                                                     <  1:2757448/150‑1 (MQ=255)
   catcatTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTTGGTTGCTCATTACGTACAGCa                                                                                                                                                    >  1:6164933/1‑150 (MQ=255)
     tcatTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAAATACGTCCTGCATc                                                                                                                                                  >  2:51423/1‑150 (MQ=255)
        ttttGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATTACGTCCTGCATCCGc                                                                                                                                               >  1:3585884/1‑150 (MQ=255)
             tttATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAAt                                                                                                                                          <  1:6715243/150‑1 (MQ=255)
             tttATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAAt                                                                                                                                          <  1:4704096/150‑1 (MQ=255)
             tttATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAAt                                                                                                                                          <  1:2461514/150‑1 (MQ=255)
             tttATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAAt                                                                                                                                          <  2:5866186/150‑1 (MQ=255)
                aTGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTGTTTGGTTGCTAATTACGTCCTGCAAACGCTTAATCTg                                                                                                                                       >  2:4869196/1‑150 (MQ=255)
                  ggTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCCCTTAATCTGGt                                                                                                                                     >  2:5200942/1‑150 (MQ=255)
                   gTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTc                                                                                                                                    <  2:1087143/150‑1 (MQ=255)
                   gTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTGTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTc                                                                                                                                    >  2:6984545/1‑150 (MQ=255)
                     tCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGTCGTCCTGCAACCGCTTAATATGGTCgg                                                                                                                                  >  1:981693/1‑150 (MQ=255)
                     tCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCgg                                                                                                                                  <  1:2055323/150‑1 (MQ=255)
                      ccACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCggg                                                                                                                                 <  1:5347657/150‑1 (MQ=255)
                      ccACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTGTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCggg                                                                                                                                 >  2:5857340/1‑150 (MQ=255)
                       cacaGCTGCGGCTTGCCTTATTCCTGCACATCTGCCATCACGGTTGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGc                                                                                                                                <  1:2968483/150‑1 (MQ=255)
                        acaGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGTCGTCCTGCATCCGCTTAATCTGGTCGGGCg                                                                                                                               >  2:2428743/1‑150 (MQ=255)
                        acaGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCg                                                                                                                               >  2:1716586/1‑150 (MQ=255)
                          aGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTa                                                                                                                             <  2:7233808/150‑1 (MQ=255)
                          aGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTa                                                                                                                             <  2:116989/150‑1 (MQ=255)
                               cGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACat                                                                                                                        <  2:3905845/150‑1 (MQ=255)
                                ggCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTGTTGGGTTGCTAATGACGTCCTGCATACGCTTAATATGGTCGGGCGTAGACata                                                                                                                       >  2:6060213/1‑150 (MQ=255)
                                   ttGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTGTTGGGTTGCTAATGACGTCCAGCAACCGCTTAATCTGGTCGGGCGTCGACCTATaa                                                                                                                    >  1:1394611/1‑150 (MQ=255)
                                      ccTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATcc                                                                                                                 <  1:3413991/150‑1 (MQ=255)
                                      ccTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATcc                                                                                                                 <  1:3259248/150‑1 (MQ=255)
                                       cttcttCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCGAATGACGTCCGGCATCCCCTTAATCTGGTCGGGCGTAGACATATAATccc                                                                                                                >  1:7074301/1‑150 (MQ=255)
                                        ttcttcCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTGTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATcaca                                                                                                               >  2:5039581/1‑150 (MQ=255)
                                         tcttcCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTGTTGGGGTGCGAATGACGTCCTGCATCCGCGTAATCGGGTCGGGCGTAGACATATAATCCcac                                                                                                              >  1:1108560/1‑150 (MQ=255)
                                            tcCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAgg                                                                                                           <  1:5299605/150‑1 (MQ=255)
                                             ccTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTTGGTTGCTAATGAAGTCCTGCATACGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAggg                                                                                                          >  1:5670709/1‑150 (MQ=255)
                                             ccTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAggg                                                                                                          <  2:5651546/150‑1 (MQ=255)
                                             ccTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAggg                                                                                                          <  1:769687/150‑1 (MQ=255)
                                                gCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGa                                                                                                       <  1:6033163/150‑1 (MQ=255)
                                                   caTCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCa                                                                                                    <  1:1253424/150‑1 (MQ=255)
                                                    aTCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAg                                                                                                   <  2:940406/150‑1 (MQ=255)
                                                     acTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGa                                                                                                  <  1:7170806/149‑1 (MQ=255)
                                                      cTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGaa                                                                                                 <  1:2921020/150‑1 (MQ=255)
                                                       tGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAAc                                                                                                <  2:4308572/150‑1 (MQ=255)
                                                       tGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAAACCCACAGGGTGAGCAGAAc                                                                                                >  2:437439/1‑150 (MQ=255)
                                                         ccATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATACGCTTAATCTGGTCGGGCGTAGACCAATAATCCCACAGGGGGGGCAGAACCt                                                                                              >  2:5547768/1‑150 (MQ=255)
                                                         ccATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGAAGTCCTGCATCCGCGTAATAAGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCt                                                                                              >  1:4410579/1‑150 (MQ=255)
                                                          cATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATACCACAGGGTGAGCAGAACCtt                                                                                             >  2:326337/1‑150 (MQ=255)
                                                          cATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTGTTTGGTTGCTAATGACGTCCTGCCTACGCTTAAACTGGTCGGGCGGAGACAAATAATCCCACAGGGGGAGCAGAAccgt                                                                                             >  2:1313568/1‑148 (MQ=255)
                                                          aaTCACGGGGGAGTTGCTGTCAAATAGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTTGGTTGCTAATTTAGAACATCAAAACGATAATCTTGTCCGGCGTATAAAAATAAACCCACAGGGTGAGCAGAAccat                                                                                             >  2:2084419/2‑148 (MQ=255)
                                                            tCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTTGGTTGCTAATTACGTCCAGCATCCGCTTAATCTGGTCGGGCGGAGACATATAATCCCACAGGGTGAGCAGAACCattt                                                                                           >  1:4573336/1‑150 (MQ=255)
                                                            tCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTACTGCCTACCCTTAATCTGGTCGGGCGTAGACATATAATACCACAGGGTGAGCAGAACCtttt                                                                                           >  2:866898/1‑150 (MQ=255)
                                                             cACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCttttt                                                                                          <  1:2406549/150‑1 (MQ=255)
                                                             cACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCttttt                                                                                          <  1:1688873/150‑1 (MQ=255)
                                                              aCGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCtttttt                                                                                         <  1:1988446/150‑1 (MQ=255)
                                                                gggggAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGc                                                                                       <  1:1981178/150‑1 (MQ=255)
                                                                  gggAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCgg                                                                                     <  2:7171889/150‑1 (MQ=255)
                                                                  gggAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCgg                                                                                     <  1:6317726/150‑1 (MQ=255)
                                                                   ggAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGc                                                                                    <  1:5136500/150‑1 (MQ=255)
                                                                     aGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTa                                                                                  <  2:699996/150‑1 (MQ=255)
                                                                       ttGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATATGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCGTTTTTGCGGCTaca                                                                                >  2:3504791/1‑150 (MQ=255)
                                                                       ttGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGAAGTCCTGCATACGCTTAATCTGGTCGGGCGTAGACATATAATACCACAGGGTGAGCAGAACCTTTTTTGCGGCTaca                                                                                >  2:2903526/1‑150 (MQ=255)
                                                                        tGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATTACGTCCTGCATCCGCTTAATCTGGTAGGGCGTAGTCATATAATCCCACAGGGTGAGCAGAACCGTTTTTGCGGCTacaa                                                                               >  2:2587361/1‑149 (MQ=255)
                                                                        tGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTacac                                                                               <  1:889935/150‑1 (MQ=255)
                                                                          cTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTa                                                                             <  2:5720811/150‑1 (MQ=255)
                                                                          cTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATACGCTTAATCTGGTCGGGCGTAGACATATAAACCCACAGGGTGATCAGAACCTTTTTTTAGGCTacacaa                                                                             >  2:2159467/1‑148 (MQ=255)
                                                                               aaaTCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTTGGTTGCTAATGACGTCCAGCATCCGCTTAATCTGGTCGGGCGTAGACATATAAAACCACAGGGTGTGCAGAACCTTTTTTTGGGCGACAAAAtttac                                                                        >  1:3029317/1‑148 (MQ=255)
                                                                               aaaTCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATTACGTCCAGCAACCGCTTAATAAGGTCGGGCGGAGACCTATAAAACCACAGGGTGAGCAGAAACATTTTTTCGGCTACACAATTTcc                                                                        >  1:4131048/1‑150 (MQ=255)
                                                                                aaTCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTccc                                                                       >  2:5341547/1‑150 (MQ=255)
                                                                                 aTCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTTGGTTGCGAATGTCGTCCAGCATACCCAAAATATGGTCGGGCGTAGACAAATAATAACACAGGGTGAGCCGAACCATTTTTTCGGCGACACAAATTACct                                                                      >  1:6776410/1‑150 (MQ=255)
                                                                                 aTCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAAATACGTACTGCATACGCTTAATATTGTCGGGCGTAGTAATATAAAACCACAGGGTGAGCAGAAACATTTTTTCGGCTACACAAATTCCct                                                                      >  1:4843994/1‑150 (MQ=255)
                                                                                  tCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCAACCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCATTTTTGCGGCTACACTATTTCCctg                                                                     >  2:6016945/1‑150 (MQ=255)
                                                                                        ttGCACAGCATTCATTTGCTCGGGAATTCCTTACAGTCGTTTTGGTTTGCTTATGACGTCCTGCATACGCTTTATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCctgatctga                                                               <  2:3318681/149‑1 (MQ=255)
                                                                                        atGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCctgatctga                                                               <  2:3003769/150‑1 (MQ=255)
                                                                                        atGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCctgatctga                                                               >  2:531434/1‑150 (MQ=255)
                                                                                        atGCACAGAATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGGTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGTGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCctgatctga                                                               <  2:4957587/150‑1 (MQ=255)
                                                                                            acaGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACAATTTCCCTTATATGATACt                                                           >  2:4758503/1‑150 (MQ=255)
                                                                                             caGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTa                                                          <  2:3555175/150‑1 (MQ=255)
                                                                                                cattcattTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATg                                                       <  1:5960460/150‑1 (MQ=255)
                                                                                                          cACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGcc                                             <  2:1560216/150‑1 (MQ=255)
                                                                                                            cGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGcccc                                           <  2:3439826/150‑1 (MQ=255)
                                                                                                                  ttCCATCCAGTCGTTTTGGGTTGCGAATGACGTCCTGCAACCGCATAATCTGGTCGGGCGTAGAAATATAATCCCACACGGTGAGCAGAACCTTTTTTGCGGCTACACAATATCCCTGATCTGATACAAAAACCAAAAAACCCCAaaaaa                                     >  2:6390356/1‑149 (MQ=255)
                                                                                                                     catccaGTCGTTTTGGGTTGCTAATTACGTCCTGCATCCGCTTAATCTGGTCGGGCGTATACATATAAACCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTAt                                  >  2:550142/1‑150 (MQ=255)
                                                                                                                             cGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAAc                          <  1:625912/150‑1 (MQ=255)
                                                                                                                             cGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAAc                          <  2:2810828/150‑1 (MQ=255)
                                                                                                                                tttGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAg                       >  1:4556001/1‑150 (MQ=255)
                                                                                                                                         cTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGcatca              >  1:3863015/1‑150 (MQ=255)
                                                                                                                                         cTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGcatca              >  2:3288857/1‑150 (MQ=255)
                                                                                                                                         cTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGcatca              >  2:1437796/1‑150 (MQ=255)
                                                                                                                                         cTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGcatca              >  1:1810931/1‑150 (MQ=255)
                                                                                                                                         cTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGcatca              >  2:3964054/1‑150 (MQ=255)
                                                                                                                                          tAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGcatcat             >  1:1317608/1‑150 (MQ=255)
                                                                                                                                          tAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGcatcat             >  2:2329992/1‑150 (MQ=255)
                                                                                                                                           aaTGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATa            <  1:3314621/150‑1 (MQ=255)
                                                                                                                                           aaTGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATa            <  1:445/150‑1 (MQ=255)
                                                                                                                                           aaTGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATa            <  2:4131048/150‑1 (MQ=255)
                                                                                                                                             tGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATAGt          <  2:5009791/150‑1 (MQ=255)
                                                                                                                                             tGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATAGt          >  2:4724542/1‑150 (MQ=255)
                                                                                                                                              gACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATAGTg         <  2:6265801/150‑1 (MQ=255)
                                                                                                                                               aCGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATAGTgg        >  2:4535320/1‑150 (MQ=255)
                                                                                                                                               aCGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATAGTgg        <  1:6562114/150‑1 (MQ=255)
                                                                                                                                                  tCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAAGCTTTTTTGCGGCTACACTATTTGCCTGATCTGTTACTAATGCAATATAGAACCATAATCTATCCCTCATCCAGACAGCATCATAGTGGGtt     >  2:6477870/1‑150 (MQ=255)
                                                                                                                                                  tCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATAGTGGGtt     >  2:6760158/1‑150 (MQ=255)
                                                                                                                                                    cTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATAGTGGGTTTc   <  2:1804058/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGCATCATTTTGTTTATGGTTCCACAGCTGCGGCTTGCCTTCTTCCTGCACATCTGCCATCACGGGGGAGTTGCTGTCAAATCGGATATGCACAGCATTCATTTGCACGGGAATTCCATCCAGTCGTTTTGGGTTGCTAATGACGTCCTGCATCCGCTTAATCTGGTCGGGCGTAGACATATAATCCCACAGGGTGAGCAGAACCTTTTTTGCGGCTACACTATTTCCCTGATCTGATACTAATGCCATATAGCCCCATAAACTATCCCTCAACCAGACAGCATCATAGTGGGTTTCT  >  NC_000913/2998183‑2998481

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: