Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,519,313 C→A 8.6% intergenic (+89/+25) 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,3130CA 8.6% 237.5 / 12.4 94intergenic (+89/+25)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 (39/47);  new base A (5/3);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-01

GGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGA  >  NC_000913/4519167‑4519462
                                                                                                                                                  |                                                                                                                                                     
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                                                                                                                                             <  2:277724/150‑1 (MQ=255)
       ttGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGc                                                                                                                                             <  1:547369/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                                                                                                                         <  2:1007066/150‑1 (MQ=255)
                           gCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAgtg                                                                                                                         >  1:614972/1‑150 (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                                                                                    >  2:439429/1‑150 (MQ=255)
                                                                gTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCt                                                                                    >  1:292380/1‑150 (MQ=255)
                                                                    ttACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtcat                                                                                <  2:717019/150‑1 (MQ=255)
                                                                    ttACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtcat                                                                                <  1:341502/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:561192/150‑1 (MQ=255)
                                                                                   aTTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCg                                                                 <  2:1122232/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:739223/1‑150 (MQ=255)
                                                                                                   gTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCcag                                                 >  2:321760/1‑150 (MQ=255)
                                                                                                       tGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAgg                                             <  1:263355/150‑1 (MQ=255)
                                                                                                        gTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAggg                                            <  2:330524/150‑1 (MQ=255)
                                                                                                        gTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAggg                                            >  2:831939/1‑150 (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               >  2:1211734/3‑150 (MQ=255)
                                                                                                                                     aacTGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcc               >  1:138439/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:297685/1‑150 (MQ=255)
                                                                                                                                          ccccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccacc          <  2:1072870/150‑1 (MQ=255)
                                                                                                                                           cccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccacca         >  2:477221/1‑150 (MQ=255)
                                                                                                                                            ccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTACCCGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAg        >  1:75410/1‑150 (MQ=255)
                                                                                                                                                  cTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGa  >  1:1291001/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGA  >  NC_000913/4519167‑4519462

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: