Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,519,314 T→G 6.9% intergenic (+90/+24) 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,3140TG 6.9% 462.5 / 20.8 161intergenic (+90/+24)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 T (75/74);  new base G (4/7);  total (79/82)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

Only 100 of 163 total aligned reads displayed.

TCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAA  >  NC_000913/4519165‑4519463
                                                                                                                                                     |                                                                                                                                                     
tCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCt                                                                                                                                                       <  2:786558/150‑1 (MQ=255)
 cGGTGGAGTTGCATGACAAGGTCATCTGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGGCATTACCGCCCATGTTTATTTCATTCATTGTGGGTATTGTTCTTATTATCGCCGCTAATCATTAAAATCCTGCCCCATATCTa                                                                                                                                                      >  2:128878/1‑150 (MQ=37)
  ggTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAACCACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCACATATCTaa                                                                                                                                                     >  1:915363/1‑149 (MQ=255)
   gTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACa                                                                                                                                                    >  2:534424/1‑150 (MQ=255)
   gTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATAGCCGCTAATCAATAAAATCCTGCCCCATATCTACa                                                                                                                                                    >  1:81064/1‑150 (MQ=255)
        gTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATgggg                                                                                                                                               <  2:1441748/150‑1 (MQ=255)
         ttGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGc                                                                                                                                              <  2:625492/150‑1 (MQ=255)
         ttGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGc                                                                                                                                              <  1:1075214/150‑1 (MQ=255)
           gCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAACTCCTGCCCCATGTAGATCTGGGGCAg                                                                                                                                            <  1:1714078/150‑1 (MQ=37)
           gCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATGTAGATATGGGGCAg                                                                                                                                            <  2:941907/150‑1 (MQ=255)
           gCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAg                                                                                                                                            >  1:531583/1‑150 (MQ=255)
           gCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAg                                                                                                                                            >  2:50452/1‑150 (MQ=255)
            cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATGTAGATATGGGGCAgg                                                                                                                                           <  1:143292/150‑2 (MQ=255)
            cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATGTAGATATGGGGCAgg                                                                                                                                           <  1:1684134/150‑2 (MQ=255)
            cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgt                                                                                                                                           >  2:714705/1‑150 (MQ=255)
              tGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATGTAGATATGGGGCAggat                                                                                                                                         <  1:1273337/150‑4 (MQ=37)
              tGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAgttg                                                                                                                                         <  2:1482334/150‑1 (MQ=255)
                  aaGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCa                                                                                                                                     >  1:1134810/1‑150 (MQ=255)
                   aGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCAACATGGGGCAGTTGTTCAt                                                                                                                                    <  1:1496700/150‑1 (MQ=255)
                      tCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCt                                                                                                                                 <  1:347836/150‑1 (MQ=255)
                       cATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGCTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCtt                                                                                                                                >  2:1529211/1‑150 (MQ=255)
                        aTCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCttt                                                                                                                               >  1:1526839/1‑150 (MQ=255)
                         tCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCtttt                                                                                                                              <  2:1199897/150‑1 (MQ=255)
                              cagtATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAgtgt                                                                                                                         >  2:222150/4‑150 (MQ=255)
                              cATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGTGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAgtgt                                                                                                                         <  1:1420995/150‑1 (MQ=255)
                              aaTTATCTGAACATAAAACACTATCAGTAAGTTTGAGTCATTACCGACTATGTTTATTTTTTAAATTTTTGGTATTGTTCTTATTATATCAGTTAATGAATTAAATACTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAgtgt                                                                                                                         <  2:346197/149‑1 (MQ=25)
                               aTTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGtgtg                                                                                                                        <  2:1460442/150‑1 (MQ=255)
                               aTTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGATAATCAATAAAATCCTGCCCCATATCTACATGTGGAATTTGTTCATTCTTTTAgtttt                                                                                                                        >  1:356664/1‑147 (MQ=38)
                                     tGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAAtt                                                                                                                  >  1:1033865/1‑150 (MQ=255)
                                     tGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCATGCCCCATATCTACATGGGGCAGTTGTGCATTCTTTTAGTGTGGTAAtt                                                                                                                  >  1:740834/1‑150 (MQ=255)
                                      gAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTc                                                                                                                 >  2:1690897/1‑150 (MQ=255)
                                      gAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTc                                                                                                                 <  2:1033865/150‑1 (MQ=255)
                                       aaCATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTca                                                                                                                <  1:542412/150‑1 (MQ=255)
                                       aaCATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTca                                                                                                                <  1:1153899/150‑1 (MQ=255)
                                       aaCATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTca                                                                                                                <  1:1015420/150‑1 (MQ=255)
                                         cATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTcaca                                                                                                              >  1:462166/1‑150 (MQ=255)
                                         cATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTcaca                                                                                                              >  2:956407/1‑150 (MQ=255)
                                          aTAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGTGCAGTTGTTCATTCTTTTAGTGTGGTAATTCacac                                                                                                             >  1:369562/1‑150 (MQ=255)
                                          aTAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTAGTGTGGTAATTCACACg                                                                                                            <  2:1094492/150‑1 (MQ=255)
                                          aTAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAATTGTTCATTCTTTTAGTGTGGTAATTCacac                                                                                                             <  2:157266/150‑1 (MQ=255)
                                             aaaCACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGcc                                                                                                          >  1:5386/1‑150 (MQ=255)
                                                    aTCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTTGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCaaaa                                                                                                   <  2:1553755/150‑1 (MQ=255)
                                                          aaGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTca                                                                                                                >  1:1627477/1‑131 (MQ=255)
                                                             ttGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGACg                                                                                          >  2:436520/1‑150 (MQ=255)
                                                              tGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGt                                                                                         <  2:820092/150‑1 (MQ=255)
                                                               ggAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGtt                                                                                        >  2:1634352/1‑150 (MQ=255)
                                                                gAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTTTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTc                                                                                       >  2:685524/1‑150 (MQ=255)
                                                                  gTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCt                                                                                     >  1:983872/1‑150 (MQ=255)
                                                                     aTTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtca                                                                                  <  2:462166/150‑1 (MQ=255)
                                                                       tACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtcatc                                                                                <  2:274050/150‑1 (MQ=255)
                                                                         ccGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaa                                                                              >  2:796050/1‑150 (MQ=255)
                                                                         acGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaa                                                                              <  1:1244176/149‑1 (MQ=255)
                                                                          cGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaac                                                                             <  2:420571/150‑1 (MQ=255)
                                                                           gACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaaca                                                                            >  2:581247/1‑150 (MQ=255)
                                                                               aTGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCa                                                                        <  1:50452/150‑1 (MQ=255)
                                                                               aTGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCa                                                                        <  2:875808/150‑1 (MQ=255)
                                                                                tGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAg                                                                       >  1:680653/1‑150 (MQ=255)
                                                                                     aTTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCg                                                                  >  1:1250505/1‑150 (MQ=255)
                                                                                       ttCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTg                                                                <  2:736492/150‑1 (MQ=255)
                                                                                        tCATACATTGTGTGTATTGTTCTTATTATCGCCGCTAATAAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGa                                                               <  2:294037/150‑1 (MQ=255)
                                                                                        tCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGa                                                               <  1:71099/150‑1 (MQ=255)
                                                                                         cATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGAc                                                              >  2:1031972/1‑150 (MQ=255)
                                                                                          aTACATTGTGTGTATTGTTATTATTATCTTCTCTAATAAATATAATTCTGCCAAATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCTCACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACa                                                             <  1:99603/150‑1 (MQ=255)
                                                                                           tACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACat                                                            <  2:634083/150‑1 (MQ=255)
                                                                                              aTTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACCCGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATAtc                                                         >  2:1520410/1‑150 (MQ=255)
                                                                                                tgtgGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATAtcct                                                       >  2:135400/1‑150 (MQ=255)
                                                                                                 gtgGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCAAGCAAAAACTCTGACGTTCCTTCATCAACAATCAGGTCCGTGACATAtcctc                                                      >  2:1389164/1‑150 (MQ=255)
                                                                                                  tgGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCGTTCATCAACAATCAGGTCCGTGACCTAtcctcc                                                     >  1:760016/1‑150 (MQ=255)
                                                                                                  tgGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATAtcctcc                                                     <  2:657131/150‑1 (MQ=255)
                                                                                                   gggTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTccc                                                    >  2:956244/1‑150 (MQ=255)
                                                                                                   gggTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTccc                                                    <  1:882493/150‑1 (MQ=255)
                                                                                                    ggTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCca                                                   <  2:983872/150‑1 (MQ=255)
                                                                                                    ggTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCca                                                   <  1:568555/150‑1 (MQ=255)
                                                                                                        ttGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCcagcag                                               <  1:944310/150‑1 (MQ=255)
                                                                                                          gTTCTTATAATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAggg                                             <  2:216559/150‑1 (MQ=255)
                                                                                                           ttCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGc                                            <  1:1305916/150‑1 (MQ=255)
                                                                                                            tCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCa                                           >  2:1168191/1‑150 (MQ=255)
                                                                                                                attatCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCa                                       <  1:1549013/150‑1 (MQ=255)
                                                                                                                   atCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACg                                    >  1:1484932/1‑150 (MQ=255)
                                                                                                                   atCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACg                                    <  2:5386/150‑1 (MQ=255)
                                                                                                                      gccgcTAATCAATAAAATCCTGCACCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACCCGCCAGCAAAAACTCTGCCGTTCCTTCATCATCAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACCTTg                                 >  1:30437/1‑150 (MQ=255)
                                                                                                                          cTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTc                             <  1:47962/150‑1 (MQ=255)
                                                                                                                            aaTCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCAt                           <  1:956407/150‑1 (MQ=255)
                                                                                                                               cAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTACTTCATAAACAATCAGGTCCGTGACAAATCCTCCCAGCAGGGCACCCAACGTTGAGTCATAGc                        >  1:859210/1‑150 (MQ=255)
                                                                                                                               aaaTAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGc                        <  2:1521308/149‑1 (MQ=255)
                                                                                                                                    aaaTCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCtt                   >  1:98116/1‑150 (MQ=255)
                                                                                                                                     aaTCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCATCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCttt                  <  1:669818/150‑1 (MQ=255)
                                                                                                                                       tgcTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGTTAATTCACACGCCAGCAATAACTCTGCCGTTCCTTCATCAACAATCATGTCCGTGACATATCCTCCCAGAATGGCACCCAACGTTGCGTCATAGCCTCTTTcc                <  1:1326216/148‑1 (MQ=255)
                                                                                                                                       aacTGCCCCATGTAGATATGGGGCAGTTGTTCATTATTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATACTACCAGGAGGGCACCCAACGTTGCGTCATAGCCTCTTTcc                >  2:285233/3‑150 (MQ=255)
                                                                                                                                        acTGCCCCATGTAGATATGGGGCAGTTGTACATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCCTCAACAATCAGGTCCGTGACATATCCTCCCAGCAGCGCCCCCAACGTTGCTGCCTAGACTCTTTccc               >  1:1538439/2‑150 (MQ=255)
                                                                                                                                        acTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTccc               >  2:1684441/2‑150 (MQ=255)
                                                                                                                                         cTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcccc              <  1:222150/150‑1 (MQ=255)
                                                                                                                                          tGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCcca             >  2:1406828/1‑150 (MQ=255)
                                                                                                                                           gCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccac            >  1:1619122/1‑150 (MQ=255)
                                                                                                                                             cccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCCGCAAAAACACTGCCGTTCCTTCATCACCTATCAGGTCCGTGCCATATCCTCCCAGCAGGGCCCCGAACGTTGCGTCATCGCCTCGGTCCccacca          >  1:698793/1‑150 (MQ=255)
                                                                                                                                             cccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATAtcctc                                                      <  2:1292856/106‑1 (MQ=255)
                                                                                                                                               cATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGc        <  1:585496/150‑1 (MQ=255)
                                                                                                                                                 tatCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCaa      >  2:497529/1‑150 (MQ=255)
                                                                                                                                                   tCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAgg    >  1:1472717/1‑150 (MQ=255)
                                                                                                                                                     tACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCATCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGaa  >  1:727101/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAAGGAA  >  NC_000913/4519165‑4519463

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: