Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,629,437 T→C 6.5% intergenic (+19/+16) ydfZ → / ← ydfI selenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,629,4370TC 6.5% 316.2 / 8.0 107intergenic (+19/+16)ydfZ/ydfIselenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (55/45);  new base C (2/5);  total (57/50)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 109 total aligned reads displayed.

TGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCA  >  NC_000913/1629289‑1629586
                                                                                                                                                    |                                                                                                                                                     
tGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATg                                                                                                                                                      >  2:938830/1‑150 (MQ=255)
 ggCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACGTGATTCGTCTCGGCATGATCTATGGGTGTGAATGCCGCCGAggg                                                                                                                                                     >  1:57761/1‑150 (MQ=255)
 ggCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATgg                                                                                                                                                     >  2:101588/1‑150 (MQ=255)
     aGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcggc                                                                                                                                                 <  1:1141798/150‑1 (MQ=255)
     aGATCCTGTCGATTGATACTGAAGGGCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcggc                                                                                                                                                 <  2:1142020/150‑1 (MQ=255)
        tCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAtt                                                                                                                                              >  1:228313/1‑150 (MQ=255)
        tCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAtt                                                                                                                                              <  1:1460051/150‑1 (MQ=255)
         ccTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGGCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttc                                                                                                                                             <  2:102733/150‑2 (MQ=255)
         ccTGTCGATTGACACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTg                                                                                                                                             <  1:949129/150‑1 (MQ=255)
           tGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACATGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttcac                                                                                                                                           <  2:852201/150‑4 (MQ=255)
           tGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAAACGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACATGATTCGTCTCGGCATGAACTAAGAGAGTGAAATCCGCCGATGGCGTCATTGCt                                                                                                                                           >  1:1379573/1‑150 (MQ=255)
            gTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttcaca                                                                                                                                          <  1:1463667/150‑5 (MQ=255)
               gATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAAGTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttttt                                                                                                                                       <  2:134748/150‑1 (MQ=255)
               gATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTAGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCGTGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttttt                                                                                                                                       >  1:732689/1‑150 (MQ=255)
                aTTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCTGAAAAACAGTAGTAGTTGAAGGTTGTGAAGAGAACCTGGCACCACTGGACATGATTCGTCTAGTCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTTCtttttt                                                                                                                                      <  1:710872/150‑1 (MQ=255)
                 ttGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTa                                                                                                                                     >  1:602142/1‑150 (MQ=255)
                   gATACTGAAGGTCTGCCCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGTCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCGTTTTTACt                                                                                                                                   >  2:720765/1‑150 (MQ=255)
                     tACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTc                                                                                                                                 >  1:371702/1‑150 (MQ=255)
                     tACTGAAGGTCTGACCGCTGAGCAAATACGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAGTGCTTTTTTACTTc                                                                                                                                 >  1:377898/1‑150 (MQ=255)
                      aCTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCa                                                                                                                                <  1:645198/150‑1 (MQ=255)
                       cTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCAc                                                                                                                               >  2:1373211/1‑150 (MQ=255)
                       cTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGGTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCAc                                                                                                                               <  1:1199453/150‑1 (MQ=255)
                         gAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACgg                                                                                                                             >  2:1268276/1‑150 (MQ=255)
                         gAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACgg                                                                                                                             >  1:913958/1‑150 (MQ=255)
                         gAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGGTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACgg                                                                                                                             <  1:105043/150‑1 (MQ=255)
                               cTGCCCGCAGCGGAAATGCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGTCCTGATTAGTATCCGCATGAACTAAGCGTGTTAATGCAGCCGATGGCGGCATTGCTTGTTTAATTCCCGGAATAtt                                                                                                                       >  2:439235/1‑150 (MQ=255)
                                  aCCGCTGAGCACATCCTGCGCGGAAACACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAAATAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTg                                                                                                                    <  2:868567/150‑1 (MQ=255)
                                   ccGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGc                                                                                                                   >  1:359197/1‑150 (MQ=255)
                                   ccGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGc                                                                                                                   <  2:1221425/150‑1 (MQ=255)
                                   ccGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGc                                                                                                                   <  1:520579/150‑1 (MQ=255)
                                    cGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGcc                                                                                                                  <  1:405425/150‑1 (MQ=255)
                                         aGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGt                                                                                                             >  2:142920/1‑150 (MQ=255)
                                           cAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCg                                                                                                           <  2:1481904/150‑1 (MQ=255)
                                             aaTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCt                                                                                                         >  2:79932/1‑150 (MQ=255)
                                               tCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCttt                                                                                                       >  2:1294568/1‑150 (MQ=255)
                                                ccgccgCGGAAAAACCGTAGTTGTTGACGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTc                                                                                                      <  1:572652/150‑1 (MQ=255)
                                                ccgccgCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTc                                                                                                      <  1:565509/150‑1 (MQ=255)
                                                 cgccgcGGAAAAACCGTAGTTGTTGAATGTTGTGAAGAGAAACTGTCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                     <  1:1032664/150‑1 (MQ=255)
                                                 cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                     <  1:1330671/150‑1 (MQ=255)
                                                     gcgGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGcc                                                                                                 >  2:878315/1‑150 (MQ=255)
                                                        gAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATg                                                                                              >  1:1389261/1‑150 (MQ=255)
                                                            aaCCGTAGTTGTTGCAGGTTGTGAAGCGAACCTGGCACCACTGGACCTGAGTCGGCTCGGCATGTACTAAGCGTGTGTAAGCCGCCGATGGCGGCATTGCTTTTTTACTTCACTGAATATTTTGCCACGGTCGCATTCGCGCCATGCGCt                                                                                          <  1:1092140/150‑1 (MQ=255)
                                                            aaCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGCTGGCGGCAATGCTTTTTTACTACACGGAATATTTTGCCACTGTCGCTTTCGCGCCATGCGCt                                                                                          >  1:911332/1‑150 (MQ=255)
                                                               cGTAGTTGTTGCAGGTTGGGAAGAGTAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTTAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaat                                                                                       <  1:1176608/150‑1 (MQ=255)
                                                                gTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaata                                                                                      <  2:486524/150‑1 (MQ=255)
                                                                gTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCAAGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaata                                                                                      >  2:1163869/1‑150 (MQ=255)
                                                                 tAGTTGTTGAAGGTTGTGAAGAGAAAATGGCACCACTGGACCTGATTCGTCTGGTCATGTACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACTGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaataa                                                                                     <  2:1172003/150‑1 (MQ=255)
                                                                     tgttgAAGGTTGTGAAGAGAAACTGGCCCCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGAc                                                                                 >  2:498328/1‑150 (MQ=255)
                                                                        tgAAGGTTGTGAAGAGAAACTGGCACCACTGGCCCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAg                                                                              >  1:168144/1‑150 (MQ=255)
                                                                        tgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAg                                                                              <  1:142920/150‑1 (MQ=255)
                                                                        tgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAg                                                                              <  1:1522365/150‑1 (MQ=255)
                                                                        tgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAg                                                                              >  1:1615850/1‑150 (MQ=255)
                                                                           aGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTAc                                                                           >  2:1345617/1‑150 (MQ=255)
                                                                                   aaGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGt                                                                   >  2:922030/1‑150 (MQ=255)
                                                                                   aaGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGt                                                                   <  2:913958/150‑1 (MQ=255)
                                                                                    agagAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTc                                                                  <  1:88363/150‑1 (MQ=255)
                                                                                    agagAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTc                                                                  >  2:1156741/1‑150 (MQ=255)
                                                                                     gagaAACTGGATCCACTGGACCTGATTCGTCTCGGCATCAACTAAGCGTGTGATTGCCGCCGATGGCGGCATTGCTTTTTTACTTCTCGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCa                                                                 <  1:1159256/150‑1 (MQ=255)
                                                                                      agaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCTTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCAc                                                                <  1:931089/150‑1 (MQ=255)
                                                                                       gaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCCCGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCAct                                                               >  2:72085/1‑150 (MQ=255)
                                                                                       gaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCAct                                                               >  1:1581277/1‑150 (MQ=255)
                                                                                             ggCACCTCTGGACATGATTCGTCTCGGCGTGAACTAAGCGTGGTTATGCATACGATTGCGGGATTGCTTTTTTACTTGACGGAATATTTTGCAACGGTCGCTTTCGCTACATGCGCTAATAAAGACAAGTACTTTTCCGTCACTCTTGCa                                                         <  2:732689/150‑1 (MQ=255)
                                                                                              gCACCACTGGACCTGATTCTGCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAg                                                        <  2:347750/150‑1 (MQ=255)
                                                                                              gCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAg                                                        <  1:938830/150‑1 (MQ=255)
                                                                                               caccacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGt                                                       >  1:631755/1‑150 (MQ=255)
                                                                                                 ccacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaa                                                     >  1:716335/1‑150 (MQ=255)
                                                                                                 ccacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaa                                                     <  2:410157/150‑1 (MQ=255)
                                                                                                      ggACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaaa                                                >  1:901868/1‑150 (MQ=255)
                                                                                                      ggACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaaa                                                >  1:423761/1‑150 (MQ=255)
                                                                                                       gACCTGATTCGTCTCGGCATGAACTAAGCTTGTGAATGCCGCCGATGGAGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaaac                                               <  2:1504615/150‑1 (MQ=255)
                                                                                                         ccTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTa                                             >  1:1190643/1‑150 (MQ=255)
                                                                                                            gATTCGTCTCGTCAAGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCGTTTTTACTTCACGGAATATTTTGCCACGGGCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTg                                          <  1:927314/150‑1 (MQ=255)
                                                                                                            gATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTg                                          <  2:602142/150‑1 (MQ=255)
                                                                                                               tCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCt                                       <  1:922030/150‑1 (MQ=255)
                                                                                                                cGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCCCTCTTGCCGTAAACAAACTATTGTCTg                                      >  2:1159882/1‑150 (MQ=255)
                                                                                                                cGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTGGCCACGGTCGCTTTCGCGCCAGTGGCTAATAAAGCCAAGTACGTTTCCGTCACTCTGGCTGTAAACCAACTATTGGCTg                                      >  2:812679/1‑150 (MQ=255)
                                                                                                                   ctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGAATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCa                                   >  2:112452/1‑150 (MQ=255)
                                                                                                                      ggCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAat                                >  2:755554/1‑150 (MQ=255)
                                                                                                                       gCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatc                               >  1:594561/1‑150 (MQ=255)
                                                                                                                       gCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatc                               >  1:504167/1‑150 (MQ=255)
                                                                                                                        cATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatca                              <  2:168144/150‑1 (MQ=255)
                                                                                                                        cATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatca                              <  2:1379573/150‑1 (MQ=255)
                                                                                                                        cATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatca                              >  1:1508395/1‑150 (MQ=255)
                                                                                                                             aCTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCAcc                         >  2:477546/1‑150 (MQ=255)
                                                                                                                             aCTAAGCGTGTGAATGCCGCCGATGGCGGCATTCCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCTCCATCCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCAcc                         >  2:1048203/1‑150 (MQ=255)
                                                                                                                               tAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCaa                       <  2:1436336/150‑1 (MQ=255)
                                                                                                                                      gtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTTATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc                >  2:1097196/1‑150 (MQ=255)
                                                                                                                                      agcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc                >  2:297335/4‑150 (MQ=255)
                                                                                                                                        gAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCt              >  2:1381787/1‑150 (MQ=255)
                                                                                                                                        caaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAACCTATTGTCTGGCAAATCATCACCAAAGATCGCCt              >  2:1187729/2‑150 (MQ=255)
                                                                                                                                            gccgccGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAAt          >  1:262953/1‑150 (MQ=255)
                                                                                                                                             ccgccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATc         >  1:529385/1‑150 (MQ=255)
                                                                                                                                             ccgccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATc         <  1:1345617/150‑1 (MQ=255)
                                                                                                                                              cgccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCg        <  1:968672/150‑1 (MQ=255)
                                                                                                                                              cgccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCg        <  1:190779/150‑1 (MQ=255)
                                                                                                                                               gccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGc       >  1:1151340/1‑150 (MQ=255)
                                                                                                                                                ccgATGGCGGCATTGCGTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGcc      >  1:1397620/1‑150 (MQ=255)
                                                                                                                                                  gATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAg    <  1:101588/150‑1 (MQ=255)
                                                                                                                                                   aTGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGc   <  2:872514/150‑1 (MQ=255)
                                                                                                                                                    tGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTAGCGCCATTCGCTAATAAAGACAAGTAAGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAAAAATCACAAAAACTAGCCTTAATAGCCAGCa  >  1:1656439/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCA  >  NC_000913/1629289‑1629586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: