Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,519,316 C→T 5.6% intergenic (+92/+22) insB1 → / ← yjhU IS1 transposase B;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,519,3160CT 5.6% 386.8 / 5.8 129intergenic (+92/+22)insB1/yjhUIS1 transposase B;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element
Reads supporting (aligned to +/- strand):  ref base C (45/73);  new base T (4/3);  total (53/76)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

TGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAAG  >  NC_000913/4519169‑4519464
                                                                                                                                                   |                                                                                                                                                    
tGGAGTTGCATGACAAGGTCATCGGGCTTTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATTTTTAATTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTTCAt                                                                                                                                                    <  1:1113506/150‑1 (MQ=255)
   aGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATggg                                                                                                                                                 <  2:519821/150‑1 (MQ=255)
     ttGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGc                                                                                                                                               <  1:695944/150‑1 (MQ=255)
     ttGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCAAATCAATAAAATCCTGCCCCATGTAGATATGGGGc                                                                                                                                               <  2:931045/150‑1 (MQ=37)
      tGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCa                                                                                                                                              >  1:1367331/1‑150 (MQ=255)
        cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATGTAGATATGGGGCAgg                                                                                                                                            <  2:1512350/150‑2 (MQ=255)
        cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgt                                                                                                                                            >  1:1675775/1‑150 (MQ=255)
        cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCAATATCTACATGGGGCAgt                                                                                                                                            >  1:33606/1‑150 (MQ=255)
          tGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgttg                                                                                                                                          >  1:1570020/1‑150 (MQ=255)
          tGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATAAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATGTAGATATGGGGCAggat                                                                                                                                          <  2:528478/150‑4 (MQ=37)
           gACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgttgt                                                                                                                                         <  2:545372/150‑1 (MQ=255)
            aCAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgttgtt                                                                                                                                        >  1:1258014/1‑150 (MQ=255)
            aCAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgttgtt                                                                                                                                        >  1:1391897/1‑150 (MQ=255)
             cAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgttgtta                                                                                                                                       >  1:1562957/1‑149 (MQ=255)
                 gTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTc                                                                                                                                   <  1:1563869/150‑1 (MQ=255)
                   cATCGGGCATTATCTGAACATAAAACACTATCAGAAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCtt                                                                                                                                 <  1:802582/150‑1 (MQ=255)
                    aTCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCttt                                                                                                                                <  2:1447082/150‑1 (MQ=255)
                    aTCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCttt                                                                                                                                >  2:716369/1‑150 (MQ=255)
                     tCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCtttt                                                                                                                               <  2:1154284/150‑1 (MQ=255)
                     tCGGGCATTATCTGAACATAAAACACTATAAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATTCATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCtttt                                                                                                                               <  1:1445455/150‑1 (MQ=255)
                      cGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTa                                                                                                                              >  2:804875/1‑150 (MQ=255)
                       gggCATTATCTGAACATAAAACCCTATCAGTAAGCTGGAGGCCTTCCCGACCCTGTTTCATTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAg                                                                                                                             <  1:254896/150‑1 (MQ=255)
                       gggCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAg                                                                                                                             <  2:1665333/150‑1 (MQ=255)
                         gCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAgtg                                                                                                                           >  1:1273181/1‑150 (MQ=255)
                         gCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGGCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAgtg                                                                                                                           <  1:1385638/150‑1 (MQ=255)
                            ttATCTGAACATAAAACACTATCAGTAAGTTGGATTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTgg                                                                                                                        <  1:1163397/150‑1 (MQ=255)
                            ttATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTgg                                                                                                                        <  1:922056/150‑1 (MQ=255)
                             tATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGt                                                                                                                       <  2:1410337/150‑1 (MQ=255)
                             tATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGt                                                                                                                       <  1:1623845/150‑1 (MQ=255)
                                  gAACATAACACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAAAAAACTCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTc                                                                                                                  <  2:901041/150‑1 (MQ=255)
                                  gAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAAtta                                                                                                                  >  1:1744004/1‑149 (MQ=255)
                                    aCATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCTGTTGTTCATTCTTTTAGTGTGGTAATTcac                                                                                                                >  1:1179220/1‑150 (MQ=255)
                                     cATAATACTCTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTcaca                                                                                                               <  1:716682/150‑1 (MQ=255)
                                      aTAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCacac                                                                                                              >  2:176626/1‑150 (MQ=255)
                                       tAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACg                                                                                                             >  1:1693361/1‑150 (MQ=255)
                                          aaCACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATACTGCCAAATATCTAAATGGGGCAGTTGTTCATTCTTTTATTGTGGTAATTCACACGCCa                                                                                                          >  1:516801/1‑150 (MQ=255)
                                              cTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCaa                                                                                                      >  2:1353851/1‑150 (MQ=255)
                                                aTCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCAAATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACAAGCCAGCaaaa                                                                                                    >  2:202828/1‑150 (MQ=255)
                                                    gTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAActc                                                                                                >  2:1235575/1‑150 (MQ=255)
                                                    gTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAActc                                                                                                >  1:1044323/1‑150 (MQ=255)
                                                        tttGTAGTCATTACCGACAATGTTTATTTGTTACATTGTGGGTATTGTTCTTATTATCGCCGGTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGcc                                                                                            <  2:583782/149‑1 (MQ=255)
                                                        gTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGcc                                                                                            <  1:414329/150‑1 (MQ=255)
                                                            gAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTc                                                                                        <  1:949726/150‑1 (MQ=255)
                                                            gAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTc                                                                                        <  1:1071430/150‑1 (MQ=255)
                                                               tCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCGCATATCTACATGGGGCAGTTGTTCATTCTTATAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCtt                                                                                     >  2:1664025/1‑150 (MQ=255)
                                                                cATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTCTCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtc                                                                                    <  2:1153175/150‑1 (MQ=255)
                                                                 aTTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtca                                                                                   >  1:521845/1‑150 (MQ=255)
                                                                   tACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtcatc                                                                                 <  2:1675775/150‑1 (MQ=255)
                                                                   tACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtcatc                                                                                 <  2:1031129/150‑1 (MQ=255)
                                                                     ccGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaa                                                                               <  2:752490/150‑1 (MQ=255)
                                                                       gACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCAAATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCATCAAAAACTCTGCCGTTCCTTCATcaaca                                                                             <  1:395707/150‑1 (MQ=255)
                                                                         tcATGTTTTTTTCATAAATTGTGGGTTTTGTTCTTATTATAGCGGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTATTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAAt                                                                           <  2:1306227/149‑1 (MQ=255)
                                                                          cATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATc                                                                          <  2:1391897/150‑1 (MQ=255)
                                                                           aTGTTTATTTCATTCATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCa                                                                         <  2:1066734/150‑1 (MQ=255)
                                                                           aTGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCa                                                                         <  1:1502050/150‑1 (MQ=255)
                                                                           aTGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCACCATATCTACATGTGGCAGTTGTTCATTCTTTTAGTGTGTTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCa                                                                         <  1:976432/150‑1 (MQ=255)
                                                                           aTGTTTATTTAATACATTGTGGGTATTGTTCTTATTATCTCAGCTAATCTATAAAATCCTGCCCCATATCTACATTGGGCAGTTGTTCATTCTTTTAGTGTTGTAATTCACACGCCAGCAAAAACTCTGCCGTACCTTCATCAACAATCa                                                                         <  1:291020/150‑1 (MQ=255)
                                                                                  tttCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGATcctt                                                                  >  2:155068/1‑148 (MQ=255)
                                                                                   ttCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTg                                                                 <  1:855442/150‑1 (MQ=255)
                                                                                     cATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGAc                                                               <  2:1394623/150‑1 (MQ=255)
                                                                                             gtgGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATAtcctc                                                       <  1:176626/150‑1 (MQ=255)
                                                                                                ggTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCca                                                    <  1:876777/150‑1 (MQ=255)
                                                                                                  tATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATATACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCATGTCCGTGACATATCATCCcagc                                                  >  1:169063/1‑150 (MQ=255)
                                                                                                      gTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAggg                                              >  2:1239929/1‑150 (MQ=255)
                                                                                                      gTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAggg                                              >  2:1056760/1‑150 (MQ=255)
                                                                                                       ttCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCAAATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCAGCAGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCACAGCAGGGc                                             >  1:1372399/1‑150 (MQ=255)
                                                                                                        tCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGCGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCa                                            <  2:1592084/150‑1 (MQ=255)
                                                                                                         cTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTTGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCATGTCCGTGACATATCCTCCCAGCAGGGCAc                                           <  1:1312055/150‑1 (MQ=255)
                                                                                                         cTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGGAATTCACACGCACGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGccc                                           >  1:1151625/1‑148 (MQ=255)
                                                                                                          ttattaTCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCAcc                                          >  1:708899/1‑150 (MQ=255)
                                                                                                           tattatCGCCGCTAATCAATAAAATCCTGCCCCATATCTTGATGGGGCAGTTGTTCATTCTTGTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCGTTCATCAACACTAAGGTCCGTGACATATCCACCCAGCAGGGCGccc                                         >  2:567764/1‑150 (MQ=255)
                                                                                                           tattatCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGGCCGTGACATATCCTCCCAGCAGGGCAccc                                         <  2:277994/150‑1 (MQ=255)
                                                                                                                tCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCACCGt                                    >  1:1738048/1‑150 (MQ=255)
                                                                                                                 cgccgcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGtt                                   >  1:1418469/1‑150 (MQ=255)
                                                                                                                  gccgcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTg                                  <  2:1314006/150‑1 (MQ=255)
                                                                                                                   ccgcTAATCAATAAAATCCTTCCCCATATCTACATGGGGCAGTTGGTCTTTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGc                                 <  1:882874/150‑1 (MQ=255)
                                                                                                                   ccgcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGc                                 <  2:1273181/150‑1 (MQ=255)
                                                                                                                    cgcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCg                                <  1:312256/150‑1 (MQ=255)
                                                                                                                     gcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGt                               >  1:1467353/1‑150 (MQ=255)
                                                                                                                     gcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGt                               <  1:1303589/150‑1 (MQ=255)
                                                                                                                      cTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTc                              <  2:1419497/150‑1 (MQ=255)
                                                                                                                        aaTCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCAt                            >  2:568027/1‑150 (MQ=255)
                                                                                                                             aTAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCct                       >  1:1821996/1‑150 (MQ=255)
                                                                                                                              tAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCctc                      >  1:93969/1‑150 (MQ=255)
                                                                                                                               aaaaTCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCctct                     >  2:323389/1‑150 (MQ=255)
                                                                                                                                   tCCTTCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcc                 <  1:967693/150‑1 (MQ=255)
                                                                                                                                   tCCTGCCCCATATCTACATGGGGCAGTTGTTCACTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcc                 <  1:1352342/150‑1 (MQ=255)
                                                                                                                                    ccTGCCCCATATCTCCATGGGGCAGTTGTTCATTCTTTTAGTGTGTTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTccc                <  2:1116790/150‑1 (MQ=255)
                                                                                                                                    acTGCCCCATGTAGATATGTGGCAGTTGTTCATTATTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATAAACAATCAGGTACGACACGTATCCTCCCAGCAGGGCACCCAACGTTGCGTCAAAGCCGCTTTccc                >  1:466474/2‑150 (MQ=255)
                                                                                                                                     cTGCCGCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACTATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcccc               <  2:498514/150‑1 (MQ=255)
                                                                                                                                     cTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcccc               <  2:1173049/150‑1 (MQ=255)
                                                                                                                                     cTGCCCCAAATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcccc               <  1:276784/150‑1 (MQ=255)
                                                                                                                                         cccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccacca           <  1:1036322/150‑1 (MQ=255)
                                                                                                                                          ccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAg          <  1:1400384/150‑1 (MQ=255)
                                                                                                                                             tatCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAAGTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCaa       <  2:1626990/150‑1 (MQ=255)
                                                                                                                                             tatCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATAAGGTCCGTGACATATCCTCACAGCAGGGCACACAACGTTGCGTCATTGCCTCTTTCCCCACCTGCaa       >  1:893415/1‑150 (MQ=255)
                                                                                                                                               tCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAgg     <  2:521845/150‑1 (MQ=255)
                                                                                                                                                cTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACAAGCCAGAAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGGGACATATCCTCCCAGCAGTGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAgta    <  2:1519804/150‑3 (MQ=255)
                                                                                                                                                  aCATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAAg  >  2:1449593/1‑150 (MQ=255)
                                                                                                                                                  aCATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAAg  <  2:1127779/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAAG  >  NC_000913/4519169‑4519464

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: