Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,519,316 C→T 5.5% 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.5% 471.9 / 9.4 148intergenic (+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 (67/71);  new base T (6/2);  total (73/75)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 148 total aligned reads displayed.

GTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAAG  >  NC_000913/4519168‑4519464
                                                                                                                                                    |                                                                                                                                                    
gTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACa                                                                                                                                                     >  2:282658/1‑150 (MQ=255)
  ggAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATg                                                                                                                                                   <  1:1625750/150‑1 (MQ=255)
    aGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATggg                                                                                                                                                 >  2:1225538/1‑150 (MQ=255)
       tGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGCGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCa                                                                                                                                              >  2:77690/1‑150 (MQ=255)
       tGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCa                                                                                                                                              >  2:1067786/1‑150 (MQ=255)
        gCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATGTAGATATGGGGCAg                                                                                                                                             <  2:477990/150‑1 (MQ=255)
        gCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAg                                                                                                                                             >  2:1295242/1‑150 (MQ=255)
        gCATGACAAGGTCATCGGGCATAATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAg                                                                                                                                             <  1:1453045/150‑1 (MQ=255)
         cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgt                                                                                                                                            >  1:1865974/1‑150 (MQ=255)
          aTGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAATTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCAAAATCTACATGGGTCAgtt                                                                                                                                           >  2:1008689/1‑150 (MQ=255)
           tGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgttg                                                                                                                                          >  2:1491405/1‑150 (MQ=255)
           tGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgttg                                                                                                                                          >  1:1171261/1‑150 (MQ=255)
              cAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgttgtgc                                                                                                                                       >  1:476638/1‑148 (MQ=255)
               aaGGTCATCGTGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCa                                                                                                                                      <  1:1371622/150‑1 (MQ=255)
               aaGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTACATACATTGTGGCTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCa                                                                                                                                      <  1:79761/150‑1 (MQ=255)
                aGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCAt                                                                                                                                     <  1:854567/150‑1 (MQ=255)
                  gTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTc                                                                                                                                   <  2:1999622/150‑1 (MQ=255)
                   tCATCGTGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCt                                                                                                                                  <  1:164268/150‑1 (MQ=255)
                    cATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCtt                                                                                                                                 >  2:1687562/1‑150 (MQ=255)
                      tCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCtttt                                                                                                                               <  1:1684497/150‑1 (MQ=255)
                          gCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAgtg                                                                                                                           >  2:525369/1‑150 (MQ=255)
                            aTTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGtgtg                                                                                                                         <  2:1571153/150‑1 (MQ=255)
                            aTTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGtgtg                                                                                                                         <  1:1046933/150‑1 (MQ=255)
                             ttATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTgg                                                                                                                        >  1:986570/1‑150 (MQ=255)
                                 cTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAAt                                                                                                                    <  1:1060461/150‑1 (MQ=255)
                                  tGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAAtt                                                                                                                   <  2:1869717/150‑1 (MQ=255)
                                   gAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTc                                                                                                                  >  1:1269922/1‑150 (MQ=255)
                                     aCATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTcac                                                                                                                <  2:1369356/150‑1 (MQ=255)
                                     aCATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTcac                                                                                                                >  1:1973708/1‑150 (MQ=255)
                                      cATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTcaca                                                                                                               <  2:908012/150‑1 (MQ=255)
                                       aTAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCacac                                                                                                              >  1:1972489/1‑150 (MQ=255)
                                        tAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCacacc                                                                                                             >  2:866180/1‑149 (MQ=255)
                                        tAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACg                                                                                                             >  2:1357070/1‑150 (MQ=255)
                                            acacTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAg                                                                                                         <  2:1553946/150‑1 (MQ=255)
                                             cacTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGc                                                                                                        >  1:728974/1‑150 (MQ=255)
                                             cacTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGc                                                                                                        <  1:634447/150‑1 (MQ=255)
                                               cTATAAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCaa                                                                                                      >  2:1185745/1‑150 (MQ=255)
                                                    aGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAAct                                                                                                 >  2:44209/1‑150 (MQ=255)
                                                       aaGTTGGAGTCATTACCGACCATGTTTATTTCCTACATTGTTGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTg                                                                                              <  1:1789943/150‑1 (MQ=255)
                                                       aaGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTg                                                                                              >  2:761116/1‑150 (MQ=255)
                                                       aaGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTg                                                                                              <  2:887000/150‑1 (MQ=255)
                                                       aaGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTAATTCTTTTAGTGTGTTAATTCACACGCCAGCAAAAACTCTg                                                                                              >  2:1738353/1‑150 (MQ=255)
                                                        aGTTGGAGTCATTACCGACCATGTATATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGc                                                                                             <  1:954621/150‑1 (MQ=255)
                                                               gTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACTCCAGCAAAAACTCTGCAGTTCCt                                                                                      >  1:1267583/1‑150 (MQ=255)
                                                                 cATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtc                                                                                    <  2:1935364/150‑1 (MQ=255)
                                                                   ttACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtcat                                                                                  >  1:1033823/1‑150 (MQ=255)
                                                                      ccGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaa                                                                               <  1:676974/150‑1 (MQ=255)
                                                                            aTGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCa                                                                         <  1:531855/150‑1 (MQ=255)
                                                                            aTGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCa                                                                         <  2:8449/150‑1 (MQ=255)
                                                                             tGTTTATTTCATAGATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCACATATCTACATGGGGAAGTTGTTCATTCTTTTAGTGTGGTTATTCACACGCCAGCAAAAACTCTGCCGTGCCTTCATCAACAATCAg                                                                        >  1:270227/1‑150 (MQ=255)
                                                                              ttttATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCACCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAgg                                                                       <  1:36834/149‑1 (MQ=255)
                                                                              gttaaattCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAACTCCTGCCCCAAATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCggg                                                                       <  2:271046/144‑1 (MQ=255)
                                                                              gTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAgg                                                                       <  1:570120/150‑1 (MQ=255)
                                                                              gTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAgg                                                                       >  1:1395744/1‑150 (MQ=255)
                                                                                  aTTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTACg                                                                   >  2:1786766/1‑150 (MQ=255)
                                                                                    ttCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCATTAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTg                                                                 <  1:60955/150‑1 (MQ=255)
                                                                                      cATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCACCATATCTACATGGGGAAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGAc                                                               >  1:298687/1‑150 (MQ=255)
                                                                                         aCATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACata                                                            <  1:704012/150‑1 (MQ=255)
                                                                                              gtgGGTTTTGTTCTTATTATCGCCGCTAATAAATAAAATCCTGCCCCATATCTAGATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCAAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATAtcctc                                                       <  2:1078008/150‑1 (MQ=255)
                                                                                               tgGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATAtcctcc                                                      <  1:407302/150‑1 (MQ=255)
                                                                                                 ggTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCca                                                    <  2:1033823/150‑1 (MQ=255)
                                                                                                 ggTATTGTTCATATTATCGCCGCTAATCAATATAATCCTGCCCCATATCTACATTGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCca                                                    <  1:438425/150‑1 (MQ=255)
                                                                                                      tGTTCTTATTATCGCCGCTAATCAATAAAATACTGACCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAgg                                               <  2:109988/150‑1 (MQ=255)
                                                                                                       gTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAggg                                              <  2:67874/150‑1 (MQ=255)
                                                                                                       gTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAggg                                              <  1:221463/150‑1 (MQ=255)
                                                                                                       gTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAggg                                              >  2:987182/1‑150 (MQ=255)
                                                                                                          cTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCGGGGCAc                                           >  2:1283637/1‑150 (MQ=255)
                                                                                                          cTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCAc                                           >  2:893371/1‑150 (MQ=255)
                                                                                                          cTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCAc                                           >  2:1408845/1‑150 (MQ=255)
                                                                                                                 tCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGt                                    >  1:1674548/1‑150 (MQ=255)
                                                                                                                  cgccgcTTATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAGCAATCAGGTCCGAGACATCTCCTCCCAGCAGGGCCCCCCACGtt                                   >  2:1230080/1‑150 (MQ=255)
                                                                                                                  cgccgcTAATCAATAAAATCCTGCCCCATATCTACAAGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGtt                                   >  2:899519/1‑150 (MQ=255)
                                                                                                                  cgccgcTAATCAATAAAATCCTGCCCCATATCTACAAGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGtt                                   >  2:907311/1‑150 (MQ=255)
                                                                                                                   gccgcTAATCAATAAAATCCGGCCCCATTTCTAAATGGGGCAGTTGTGCATTCTTTTATGTTGGTAATTCACACGCCATCAAAAACTCTGCCGGTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTg                                  <  2:1832704/150‑1 (MQ=255)
                                                                                                                     cgcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCg                                <  1:1225538/150‑1 (MQ=255)
                                                                                                                     cacTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCg                                <  1:595634/148‑1 (MQ=255)
                                                                                                                        taataaATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGAAATTCACACGCCAGCAAAAACTCAGCCGATCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCCCCCAACGTTGCGGCa                             >  1:1094487/1‑150 (MQ=255)
                                                                                                                         aaTCAATAAAATCCTGCCCCATATCAACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATAAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCAt                            >  1:544366/1‑150 (MQ=255)
                                                                                                                          aTCAATAAAATCCAGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATa                           <  1:828046/150‑1 (MQ=255)
                                                                                                                            cAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGc                         <  1:761116/150‑1 (MQ=255)
                                                                                                                                 aacaacTGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCAGTTCCTTAAGCAAAAATCAGGTACGTGAAATAACCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCctcat                    >  2:385240/6‑148 (MQ=255)
                                                                                                                                 aaaTCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGTCATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCtt                    >  2:987742/1‑150 (MQ=255)
                                                                                                                                  acaacTGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCttt                   >  1:1076875/5‑150 (MQ=255)
                                                                                                                                    tCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcc                 <  1:467266/150‑1 (MQ=255)
                                                                                                                                    tCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcc                 <  1:1615111/150‑1 (MQ=255)
                                                                                                                                    tCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcc                 <  2:1171261/150‑1 (MQ=255)
                                                                                                                                    tCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTCGTAAATCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTGCCAGCAGGGCACCGAACGTTGCGTCATAGCCTCTTTcc                 <  1:282658/150‑1 (MQ=255)
                                                                                                                                     acTGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTccc                >  1:237217/2‑150 (MQ=255)
                                                                                                                                     acTGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTccc                >  2:1659593/2‑150 (MQ=255)
                                                                                                                                       tGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTccccc              >  1:1520735/1‑149 (MQ=255)
                                                                                                                                       tGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCcca              <  2:728974/150‑1 (MQ=255)
                                                                                                                                       tGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCcca              >  2:766980/1‑150 (MQ=255)
                                                                                                                                        gCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccac             >  1:868122/1‑150 (MQ=255)
                                                                                                                                          cccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccaccc           >  2:1985920/1‑149 (MQ=255)
                                                                                                                                          cccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccacca           >  2:1785754/1‑150 (MQ=255)
                                                                                                                                           ccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCAAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccaccac          >  2:783083/1‑149 (MQ=255)
                                                                                                                                            cATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGc         >  1:1862159/1‑150 (MQ=255)
                                                                                                                                              tatCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCaa       >  2:851870/1‑150 (MQ=255)
                                                                                                                                               atCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAg      <  1:478518/150‑1 (MQ=255)
                                                                                                                                                   aCATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAAg  <  2:1741481/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAAG  >  NC_000913/4519168‑4519464

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