Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,519,316 C→T 7.7% 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 7.7% 267.1 / 8.7 94intergenic (+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 (38/47);  new base T (5/2);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAAGA  >  NC_000913/4519167‑4519465
                                                                                                                                                     |                                                                                                                                                     
ggTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTAc                                                                                                                                                       >  1:1236593/1‑150 (MQ=255)
 gTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACa                                                                                                                                                      >  1:1265219/1‑150 (MQ=255)
    gAGTTGCATGAAAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATgg                                                                                                                                                   <  2:940126/150‑1 (MQ=255)
      gTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATgggg                                                                                                                                                 >  1:1072870/1‑150 (MQ=255)
       ttGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGc                                                                                                                                                <  1:547369/150‑1 (MQ=255)
       ttGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGc                                                                                                                                                <  2:277724/150‑1 (MQ=255)
          caggACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATATCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATAAATAAAATCCTGCCCCATGTAGATATGGGGCAgg                                                                                                                                             <  2:746888/147‑2 (MQ=37)
          cATGACTATGTCTTCTGGCATTATCTCAAAATAAAACTCTATCAGACAGTTTGAGTAATTAACGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATTAAATCCTGCCCCATGTAGAGATGGGGCAgg                                                                                                                                             <  2:641795/150‑2 (MQ=255)
          cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATGTAGATATGGGGCAgg                                                                                                                                             <  1:877917/150‑2 (MQ=255)
          cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgg                                                                                                                                             <  2:949686/150‑2 (MQ=255)
          cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGGATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTTCCCCATATCTACATGGGGCAgt                                                                                                                                             <  2:939118/150‑1 (MQ=255)
          aaTGAGAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTTTTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgt                                                                                                                                             <  1:108588/149‑1 (MQ=255)
             gACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTGCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgttgt                                                                                                                                          >  2:923849/1‑150 (MQ=255)
               cAAGGTCATGGGGCATTATCTCAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTc                                                                                                                                        <  2:644880/150‑1 (MQ=255)
               cAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTc                                                                                                                                        <  2:946547/150‑1 (MQ=255)
               cAAGGACATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGCAGTCATTACCGCCCATGAATATTTCATACATAGTGGGTATTGTTCTTATTATCGCAGCTAATCAATAAAATCCTGACCCATATCTACATGGGGAAgttgtta                                                                                                                                        >  1:847422/1‑149 (MQ=37)
                  ggTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCAtt                                                                                                                                     >  2:1073756/1‑150 (MQ=255)
                    gcATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCt                                                                                                                                   <  2:1249504/149‑1 (MQ=255)
                       tCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCtttt                                                                                                                                >  2:713769/1‑150 (MQ=255)
                       acGGGCATTATCTGAACATAAACCCCTATCAGTAACTTTGAGTCATTACCGACCATTTTTATTTCAAACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCtttt                                                                                                                                <  2:1218283/149‑1 (MQ=39)
                         gggCATTATCTGAACATATAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTTTGGGTTTTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAg                                                                                                                              <  2:798545/150‑1 (MQ=255)
                           gCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAgtg                                                                                                                            >  1:614972/1‑150 (MQ=255)
                           gCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAgtg                                                                                                                            <  2:1007066/150‑1 (MQ=255)
                            cATTATCTGAACATAAAACACTATCTTTAAGTTGGAGTCATTACCGTCCATTTTTATTTCATACATTGTGTGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCGCCATTTCTACATTGGGCAGTTGTTCATTCTTTTAgtgt                                                                                                                           <  1:1016097/150‑1 (MQ=38)
                              ttATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTgg                                                                                                                         <  2:699963/150‑1 (MQ=255)
                                agcTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTa                                                                                                                       <  2:851840/148‑1 (MQ=255)
                                aTCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCATTTGTTCATTATTTTATTGTGGTa                                                                                                                       >  2:238484/1‑150 (MQ=255)
                                   tGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAAtt                                                                                                                    <  2:373138/150‑1 (MQ=255)
                                    taaCATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTc                                                                                                                   <  1:608769/149‑1 (MQ=255)
                                    gAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCAGATCTACATGGGGCAGTTGTTCATTCTTTTATTGTGGTAATTc                                                                                                                   >  2:130330/1‑150 (MQ=255)
                                     aaCATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTca                                                                                                                  >  2:1041795/1‑150 (MQ=255)
                                       cATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTcaca                                                                                                                <  2:98037/150‑1 (MQ=255)
                                            aaCACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCa                                                                                                           <  1:278735/150‑1 (MQ=255)
                                               acTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCa                                                                                                        >  1:493919/1‑150 (MQ=255)
                                                 attgaGTAATGTGGAGTCATTACCGACCATGTTTAATTCATACATTGTGGTTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCaaa                                                                                                      <  1:196436/146‑1 (MQ=255)
                                                        aaGTTGGATTCATTACCGACCATTTTTATTTCATACATTGTTGGTATTTTTATTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTg                                                                                               <  1:640881/150‑1 (MQ=255)
                                                           ttGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTTTGGTAATTCACACGCCAGCAAAAACTCTGCCg                                                                                            >  2:564946/1‑150 (MQ=255)
                                                            tGGAGTCATTACCGCCCATGTTTATTTCATAAATACTGGGTATTGTTCTTATTATCGCCGCTTATCAATAAAATCCTGCCCCATATCTAGATGGCGCAGCTGGACATTCTTTTAGTGTGGTAAATCACACGCCAGCACAAACTCTGCCGt                                                                                           >  2:471938/1‑150 (MQ=255)
                                                                gTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCt                                                                                       >  1:292380/1‑150 (MQ=255)
                                                                gTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCt                                                                                       >  2:439429/1‑150 (MQ=255)
                                                                    ttACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtcat                                                                                   <  1:341502/150‑1 (MQ=255)
                                                                    ttACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtcat                                                                                   <  2:717019/150‑1 (MQ=255)
                                                                      aCCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtcatca                                                                                 >  2:378282/1‑150 (MQ=255)
                                                                       ccGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaa                                                                                <  1:124511/150‑1 (MQ=255)
                                                                         gACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAACATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaaca                                                                              <  2:808417/150‑1 (MQ=255)
                                                                         gACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaaca                                                                              <  1:731778/150‑1 (MQ=255)
                                                                            cATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAAATCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATc                                                                           <  2:686293/150‑1 (MQ=255)
                                                                             aTGTTTATTTCATACATTGTGTTTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCAATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCa                                                                          <  1:168592/150‑1 (MQ=255)
                                                                             aTGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCa                                                                          <  2:1298114/150‑1 (MQ=255)
                                                                             aTGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGAAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCAGTTCCTTCATCAACAATCa                                                                          >  1:755990/1‑150 (MQ=255)
                                                                               tttaatatCATACATTGTGGGTATGGTCCTTATTATCGCCGCTACTCAATAAACTCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCATTTACTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCCTCAACAATCAgg                                                                        <  2:30933/143‑1 (MQ=255)
                                                                               gTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAgg                                                                        >  1:903543/1‑150 (MQ=255)
                                                                               gTTAATATCCTACATTGTGGGTCTTGATCTTATTACCGCCGCTAATCAATAACATCCTTCCCCATATCTACATGGGGCAGTTGTGCATTCTTTTAGTGTGGTAAATCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAgg                                                                        <  2:150662/150‑1 (MQ=255)
                                                                                 ttATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTc                                                                      <  1:313544/150‑1 (MQ=255)
                                                                                   aTTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCg                                                                    <  2:1122232/150‑1 (MQ=255)
                                                                                   aTTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCg                                                                    <  2:561192/150‑1 (MQ=255)
                                                                                     ttCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTg                                                                  >  1:489794/1‑150 (MQ=255)
                                                                                       cATACAAAGTGGGTATACTTCTTATTATCGCCGCAAATCAAAAAAATCCTGCCCCATATCTACATGGGGCAGTTGTACATTCTTTTAGTGTGGTATTTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAAACAGGTCCGTGAc                                                                >  2:22424/1‑150 (MQ=255)
                                                                                        aTACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGaaa                                                               >  1:926957/1‑148 (MQ=255)
                                                                                        aTACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACa                                                               <  2:1236593/150‑1 (MQ=255)
                                                                                             ttGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATAtcc                                                          <  2:35127/150‑1 (MQ=255)
                                                                                             ttGTGGGTATTGTTCTTATTATCGCCGCTAATCAAAAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATAtcc                                                          <  1:238484/150‑1 (MQ=255)
                                                                                                tgGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCAAATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACCCGCCAGCAAAAACTCTGCCGTTACTTCATCAACAATCAGGTCCGTGACATAtcctcc                                                       >  2:833514/1‑150 (MQ=255)
                                                                                                 gggTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTccc                                                      <  1:57912/150‑1 (MQ=255)
                                                                                                 gggTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTGCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTccc                                                      <  2:108478/150‑1 (MQ=255)
                                                                                                   gTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCcag                                                    >  2:321760/1‑150 (MQ=255)
                                                                                                   gTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCcag                                                    >  2:739223/1‑150 (MQ=255)
                                                                                                       tGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAgg                                                <  1:263355/150‑1 (MQ=255)
                                                                                                        gTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAggg                                               >  2:831939/1‑150 (MQ=255)
                                                                                                        gTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAggg                                               <  2:330524/150‑1 (MQ=255)
                                                                                                          tCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAACAAATCTGCCGTTCCTTCATCACCAATCAGGTCCGTGACAAATCCTCCCAGCAGGGCa                                             >  2:949320/1‑150 (MQ=255)
                                                                                                                 atCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACg                                      >  2:697976/1‑150 (MQ=255)
                                                                                                                   cgccgcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGGGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCAGCACCAATCAGGGCCGTGCGATATCCTCCCAGCAGGGCCCCCAACGtt                                    >  2:208862/1‑150 (MQ=255)
                                                                                                                    gccgcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTg                                   <  1:1078095/150‑1 (MQ=255)
                                                                                                                     ccgcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGc                                  <  2:551841/150‑1 (MQ=255)
                                                                                                                          aaTCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATAAACAATAAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCaa                             >  1:1053230/1‑149 (MQ=255)
                                                                                                                           aTCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATa                            <  2:355659/150‑1 (MQ=255)
                                                                                                                            tCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAg                           >  2:198444/1‑150 (MQ=255)
                                                                                                                              aaTAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACGCACCGTTGCGTCATAGcc                         >  1:931111/1‑150 (MQ=255)
                                                                                                                               aTAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAGCATTCAGGTCCGTGACATATCCTCCCAGCAGGGCCCCCTAGCTTGCGTCCTGGccg                        >  2:208498/1‑149 (MQ=255)
                                                                                                                               aTAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCct                        >  2:1128306/1‑150 (MQ=255)
                                                                                                                                   acaacTGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGGTCCTTCATCAACAATCAGGTCCGGGACATATCCTCCCAGCAGGGCACCCAACGTGGCGTCATTGACGCttt                    >  1:1097664/5‑150 (MQ=255)
                                                                                                                                    caacTGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTc                   >  1:904885/4‑150 (MQ=255)
                                                                                                                                    aTCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTc                   <  1:1261154/150‑1 (MQ=255)
                                                                                                                                     tCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcc                  <  1:1073756/150‑1 (MQ=255)
                                                                                                                                     aacTGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcc                  >  1:138439/3‑150 (MQ=255)
                                                                                                                                     aacTGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcc                  >  2:1211734/3‑150 (MQ=255)
                                                                                                                                      ccTGCCCCATAGCTACATGGGGCAGTTGGTCATTCTGTTAGTGTGGTAATTCACACGCCAGCAAAACCTCTGCCGTTCCTGCATCAGCAAGCCGGGCCGGGCCATATCCGCCCCGCAGGGCCCCCACCGTGGCGGCCTAGCGGCGTTccc                 >  2:523345/1‑150 (MQ=255)
                                                                                                                                      acTGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTccc                 >  2:585453/2‑150 (MQ=255)
                                                                                                                                         gCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCCTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccac              >  2:1146816/1‑150 (MQ=255)
                                                                                                                                          ccccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccacc             <  2:1072870/150‑1 (MQ=255)
                                                                                                                                          ccccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccacc             >  2:297685/1‑150 (MQ=255)
                                                                                                                                           cccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccacca            >  2:477221/1‑150 (MQ=255)
                                                                                                                                            ccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTACCCGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAg           >  1:75410/1‑150 (MQ=255)
                                                                                                                                                  cTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGa     >  1:1291001/1‑150 (MQ=255)
                                                                                                                                                     cATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTTCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAAGa  <  2:788082/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAAGA  >  NC_000913/4519167‑4519465

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: