Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,496,511 A→G 10.1% intergenic (+22/‑164) leuX → / → intB tRNA‑Leu/pseudogene, integrase homology;IS, phage, Tn; Phage‑related functions and prophages; KpLE2 phage‑like element; P4‑like integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,496,5110AG10.1% 335.5 / 24.3 130intergenic (+22/‑164)leuX/intBtRNA‑Leu/pseudogene, integrase homology;IS, phage, Tn; Phage‑related functions and prophages; KpLE2 phage‑like element; P4‑like integrase
Reads supporting (aligned to +/- strand):  ref base A (53/63);  new base G (5/8);  total (58/72)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01

Only 100 of 135 total aligned reads displayed.

GCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAAC  >  NC_000913/4496362‑4496660
                                                                                                                                                     |                                                                                                                                                     
gCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCaa                                                                                                                                                       <  1:1013808/150‑1 (MQ=255)
gCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCaa                                                                                                                                                       >  1:693822/1‑150 (MQ=255)
 cGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCaaa                                                                                                                                                      >  1:775203/1‑149 (MQ=255)
 cGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAAc                                                                                                                                                      >  2:560430/1‑150 (MQ=255)
   aCAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTg                                                                                                                                                    <  2:447580/150‑1 (MQ=255)
   aCAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTg                                                                                                                                                    <  1:1031544/150‑1 (MQ=255)
   aCAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTg                                                                                                                                                    >  1:1299278/1‑150 (MQ=255)
      ggTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAgg                                                                                                                                                 >  2:628549/1‑150 (MQ=255)
      ggTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAgg                                                                                                                                                 <  2:20939/150‑1 (MQ=255)
        tataATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTc                                                                                                                                               <  2:1082375/150‑1 (MQ=255)
        tataATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTc                                                                                                                                               <  1:1142275/150‑1 (MQ=255)
            aTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTCtatt                                                                                                                                           <  1:1052994/150‑1 (MQ=255)
            aTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTCtatt                                                                                                                                           <  1:1375664/150‑1 (MQ=255)
            aTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCtttt                                                                                                                                           <  2:685103/150‑1 (MQ=255)
            aTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCtttt                                                                                                                                           >  1:952/1‑150 (MQ=255)
            aTCCACAACGTTTTCCGACTACCTCTTCAGTGCCGAAGTGGCGAAATCTGTAGAAGCAGTTGATTCAAAATCAAACGTAGAAATACGTTCCGGTTCGAGTCCTGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTCtatt                                                                                                                                           <  1:1072770/150‑1 (MQ=255)
             tCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttttt                                                                                                                                          >  1:1186711/1‑150 (MQ=255)
             tCCACAACGTGTTGCGCATACCTCTTCAGTGCCGTAGTGGCGAAATCGGTTGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTGCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTCtattt                                                                                                                                          <  2:1097382/150‑1 (MQ=255)
               cacaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTCtatttac                                                                                                                                        <  1:1542705/150‑3 (MQ=255)
                acaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTCtatttaca                                                                                                                                       <  2:905824/150‑4 (MQ=255)
                    cGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGcc                                                                                                                                   >  2:612165/1‑150 (MQ=255)
                      ttttCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAGCCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTg                                                                                                                                 <  1:1666706/150‑1 (MQ=255)
                      ttttCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTg                                                                                                                                 <  1:984312/150‑1 (MQ=255)
                          ccGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCAGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAAt                                                                                                                             >  2:71179/1‑150 (MQ=255)
                            gCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGATGTCTTTTTTTATGCCTGAAATcc                                                                                                                           <  2:1042209/150‑1 (MQ=255)
                            gCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATcc                                                                                                                           <  2:882316/150‑1 (MQ=255)
                              aTACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCTTAGAAATACGTGCCGGTTCGAGTGCGGCCTTCGGCACCAAAACTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAg                                                                                                                         <  1:706948/150‑1 (MQ=255)
                               tACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAgt                                                                                                                        >  1:34444/1‑150 (MQ=255)
                                  ctctTCAGTGCCGAAGTGTCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGTGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGtt                                                                                                                     <  2:1308335/150‑1 (MQ=255)
                                      tCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTtatc                                                                                                                 <  2:1476244/150‑1 (MQ=255)
                                       cAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTGAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTtatct                                                                                                                <  1:207849/150‑1 (MQ=255)
                                        aGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTtatcta                                                                                                               >  2:787446/1‑150 (MQ=255)
                                            ccGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACATCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCtt                                                                                                           >  1:89425/1‑150 (MQ=255)
                                                aGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTccc                                                                                                       <  2:952/150‑1 (MQ=255)
                                                aGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTGTTATTGGCTGAAATCCAGTGTTTATCTATCTTTccc                                                                                                       >  1:347334/1‑150 (MQ=255)
                                                  tGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGc                                                                                                     >  2:299086/1‑150 (MQ=255)
                                                   ggCGAAATCGGTAGACGCAGTTGATTCAAAATCACCCGTAGAAATACGTGCCGGTTCGCGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGGCTTTTTTTATGCCGGAAATCCAGGGTTTATCTATCGTTCCCggt                                                                                                    >  1:320590/1‑148 (MQ=255)
                                                   ggCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCt                                                                                                    >  1:1613866/1‑150 (MQ=255)
                                                      gataTCGGTAGACGCAGTTGATTCAATATCAACCGTAGAAATACGTGGCGGTTCGTGTCTTGCCTTCGGCACGAAAAGTATGTTAATAGACCTCAACTGAGGTCGTTTTTTATGCCTGATATCCAGTGTTTATCTATCTTTCCCGCtata                                                                                                 <  1:79998/147‑1 (MQ=255)
                                                          tCGGTAGACGCAGTTGATTCATAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTaa                                                                                             >  2:1254441/1‑150 (MQ=255)
                                                               aGACGCAGTTGATTCAAAATCAACAGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTAGTTTTTTTTGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAActctc                                                                                        >  1:1065932/1‑150 (MQ=255)
                                                                gACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAActctct                                                                                       >  2:546762/1‑150 (MQ=255)
                                                                gACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAActctct                                                                                       >  2:395570/1‑150 (MQ=255)
                                                                  cGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCa                                                                                     <  1:977376/150‑1 (MQ=255)
                                                                   gCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCaa                                                                                    >  1:709507/1‑150 (MQ=255)
                                                                     aGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAgg                                                                                  >  1:157661/1‑150 (MQ=255)
                                                                      gTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGt                                                                                 >  1:1203402/1‑150 (MQ=255)
                                                                      gTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGt                                                                                 <  2:1151765/150‑1 (MQ=255)
                                                                        tGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCa                                                                               <  2:1399383/150‑1 (MQ=255)
                                                                           ttCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAAcc                                                                            >  2:378521/1‑150 (MQ=255)
                                                                           ttCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCATGAAATCCAGTGTTTATCTATCTTTCCCCCTATATTAACTCTCTCAAGGTAAAcc                                                                            >  2:732332/1‑150 (MQ=255)
                                                                            tCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCg                                                                           >  2:771163/1‑150 (MQ=255)
                                                                             cAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCATGAAATCCAGTGTTTATCTATCTTTCCAGCTATATTAACTCTATCAAGGTCAACCGa                                                                          >  2:1014396/1‑150 (MQ=255)
                                                                              aaaaTCTACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGat                                                                         <  1:150324/150‑1 (MQ=255)
                                                                               aaaTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGata                                                                        >  2:1332796/1‑150 (MQ=255)
                                                                                aaTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGatat                                                                       <  2:371801/150‑1 (MQ=255)
                                                                                 aTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATc                                                                      <  2:1299278/150‑1 (MQ=255)
                                                                                    aaCCGTAGAAATACGTGCCGGTTCGAGTCCGGCCGTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAAc                                                                   <  2:101289/150‑1 (MQ=255)
                                                                                       cGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTa                                                                <  2:1186711/150‑1 (MQ=255)
                                                                                        gTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTAc                                                               >  2:1449382/1‑150 (MQ=255)
                                                                                          aGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACAt                                                             >  1:662339/1‑150 (MQ=255)
                                                                                          aGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACATCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTAcct                                                             >  2:1310326/1‑148 (MQ=255)
                                                                                           gAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATc                                                            >  2:325254/1‑150 (MQ=255)
                                                                                            aaaTACGTGCCGGTTCGAGTCCTGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCt                                                           <  1:1172398/150‑1 (MQ=255)
                                                                                            aaaTACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCt                                                           <  2:1332551/150‑1 (MQ=255)
                                                                                                   tGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGATGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACa                                                    <  1:1211689/150‑1 (MQ=255)
                                                                                                       ggTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATg                                                >  1:940792/1‑150 (MQ=255)
                                                                                                         ttCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGtt                                              <  1:560430/150‑1 (MQ=255)
                                                                                                         ttCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGATTATCTATCTTTCCCGCTATATTAAGTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGtt                                              <  1:286654/150‑1 (MQ=255)
                                                                                                            gAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGt                                           <  1:612165/150‑1 (MQ=255)
                                                                                                              gTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTAc                                         >  1:1361162/1‑150 (MQ=255)
                                                                                                                ccGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAg                                       >  1:1434996/1‑150 (MQ=255)
                                                                                                                 cGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAga                                      >  1:28269/1‑150 (MQ=255)
                                                                                                                   gCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAgatg                                    >  2:474408/1‑150 (MQ=255)
                                                                                                                    ccTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATAATAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAgatga                                   <  2:17321/150‑1 (MQ=255)
                                                                                                                        cGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGt                               <  1:156247/150‑1 (MQ=255)
                                                                                                                         ggCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCTCTATATTATCTCTCTCAAGGTCATCCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTa                              <  1:99727/150‑1 (MQ=255)
                                                                                                                           cACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTAtt                            <  2:865485/150‑1 (MQ=255)
                                                                                                                             ccAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTc                          >  2:305338/1‑150 (MQ=255)
                                                                                                                               aaaaGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTTTCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCg                        <  1:560319/150‑1 (MQ=255)
                                                                                                                               aaaaGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCg                        <  1:1169711/150‑1 (MQ=255)
                                                                                                                               aaaaGTATGTAAATAGACCTCAACTGAGGTATTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCACTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCg                        >  2:956988/1‑150 (MQ=255)
                                                                                                                                aatgTATGTAAATATACCTCAACTTAGGTCTTTTTTTATGCCTGAATTCGAGTGTTTCTCTATGTTTCGCGCTATATTATCTCTCGGTATGTCTAGAGATATCAACGTTCATCTACCAACATATGTTGGTACAGATGATGGTATTTCCgg                       <  1:325163/147‑1 (MQ=255)
                                                                                                                                 aaGTATGTAAATAGAACTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACACATGTTGGTACAGATGATGGTATTTCCGGt                      <  2:775203/150‑1 (MQ=255)
                                                                                                                                  ggcataaaaaaAGACCTGAGTTGAGGTCTATTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTATCTCTCTAAAGGTCAACCGATATAAACGTACATCTACCACCATATGTTGGTAAAGATGATGGTATTTCAGGtt                     >  1:102504/8‑150 (MQ=255)
                                                                                                                                   gTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTc                    <  1:305338/150‑1 (MQ=255)
                                                                                                                                   gTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTc                    <  2:1613866/150‑1 (MQ=255)
                                                                                                                                         aaaaAGACCTCAGTTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAAt              >  1:457959/1‑150 (MQ=255)
                                                                                                                                         aaaaAGACCTCAGTTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAAt              >  2:1340929/1‑150 (MQ=255)
                                                                                                                                         aaaaAGACCTCAGTTGAGGTCTTTTTTTAAGCCTGAAATCCAGTGTTTATCTATATTTCCCGCTATATTAACTCTCTCAAGCTCAACCGATATCAACGTACATCTACCAACATATGTTTGTACAGATCATGGTATTTCCGGTACGATAAt              >  2:1636593/1‑150 (MQ=255)
                                                                                                                                          aaTAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATg             <  2:34444/150‑1 (MQ=255)
                                                                                                                                          aaTAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCAGGTGCGATAATg             >  1:1490270/1‑150 (MQ=255)
                                                                                                                                            tAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCt           <  2:709507/150‑1 (MQ=255)
                                                                                                                                              gACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTg         >  1:804313/1‑150 (MQ=255)
                                                                                                                                                  tCAACTGATGTCTTTTTTTATTCCTGAAATCCAGTTTTTATCTATCTTTCCAGCTATATTAACTCTCTCAAGGTCAACAGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTAcc     <  1:329661/150‑1 (MQ=255)
                                                                                                                                                  tCAACTGATGTATTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTTTTTCCGGTTCGATAATGCTTGTAcc     <  1:876911/150‑1 (MQ=255)
                                                                                                                                                  tCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTAcc     <  1:1332796/150‑1 (MQ=255)
                                                                                                                                                  gcAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAGCCGATATCATCGTACATCTACCAACATATGTTGGTACAGATGATCGTATTTCCGGTTCGATAATACTTGTagc     >  1:221202/2‑148 (MQ=255)
                                                                                                                                                     aCTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAAc  <  2:1485215/150‑1 (MQ=255)
                                                                                                                                                     aCTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAAc  >  2:21964/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCAAC  >  NC_000913/4496362‑4496660

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: