Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,496,511 A→G 5.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,5110AG 5.1% 311.1 / 5.6 97intergenic (+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 (45/47);  new base G (3/2);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-01

CGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTAC  >  NC_000913/4496363‑4496656
                                                                                                                                                    |                                                                                                                                                 
cGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCaa                                                                                                                                                   <  3:41672/149‑1 (MQ=255)
    aGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGa                                                                                                                                               <  1:673630/149‑1 (MQ=255)
      gTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATGCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAgg                                                                                                                                             >  2:703769/1‑149 (MQ=255)
       tataATCCACAACGTTTTCCGCATACCTCTTCAATGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGt                                                                                                                                            >  1:228758/1‑149 (MQ=255)
         taATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCt                                                                                                                                          >  4:680549/1‑149 (MQ=255)
          aaTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTCta                                                                                                                                         <  1:470412/149‑2 (MQ=255)
            tCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCtttt                                                                                                                                       <  4:595712/149‑1 (MQ=255)
             ccACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttttt                                                                                                                                      <  1:336268/149‑1 (MQ=255)
             ccACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAAAAGACCTCAGTTGAGGTCtattt                                                                                                                                      <  3:739634/149‑1 (MQ=255)
               acaACGTTTTCCGCATACCTCTTCAGTGCCGATGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttttttt                                                                                                                                    >  4:650431/1‑149 (MQ=255)
               acaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttttttt                                                                                                                                    >  1:590359/1‑149 (MQ=255)
               acaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttttttt                                                                                                                                    <  3:33640/149‑1 (MQ=255)
               acaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttttttt                                                                                                                                    >  3:707172/1‑149 (MQ=255)
                caACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGAGTCAAAATCAACCGTAGAAATACGTGCCGTTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACGGAGGTATTTTTGTa                                                                                                                                   >  4:632021/1‑149 (MQ=255)
                 aaCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTAt                                                                                                                                  <  2:280133/149‑1 (MQ=255)
                 aaCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTAt                                                                                                                                  <  3:207369/149‑1 (MQ=255)
                 aaCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCCCAACTGAGGTCTTTTTTTAt                                                                                                                                  <  2:691433/149‑1 (MQ=255)
                           gCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATc                                                                                                                        >  1:500064/1‑149 (MQ=255)
                             aTACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCa                                                                                                                      >  2:46873/1‑149 (MQ=255)
                               aCCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAgt                                                                                                                    >  1:335205/1‑149 (MQ=255)
                                ccTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAgtg                                                                                                                   >  4:119531/1‑149 (MQ=255)
                                 ctctTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAgtgt                                                                                                                  >  1:660659/1‑149 (MQ=255)
                                 ctctTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAgtgt                                                                                                                  >  1:157379/1‑149 (MQ=255)
                                  tctTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACATCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGtt                                                                                                                 >  4:648351/1‑149 (MQ=255)
                                   ctTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGttt                                                                                                                <  2:257954/149‑1 (MQ=255)
                                     tCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGtttct                                                                                                              <  1:177070/149‑3 (MQ=255)
                                     tCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTtat                                                                                                              <  3:407533/149‑1 (MQ=255)
                                       aGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTtatct                                                                                                            <  1:571912/149‑1 (MQ=255)
                                         tGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTtatctat                                                                                                          <  1:265219/149‑1 (MQ=255)
                                               aGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTcc                                                                                                    <  1:1840/149‑1 (MQ=255)
                                                   gCGAAATCGGTAGACGCAGTTGTTTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGCCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCt                                                                                                >  4:216044/1‑149 (MQ=255)
                                                        aTCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATAtt                                                                                           >  3:460095/1‑149 (MQ=255)
                                                           ggTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAAc                                                                                        >  1:169945/1‑149 (MQ=255)
                                                           ggTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAAc                                                                                        >  2:42396/1‑149 (MQ=255)
                                                               gACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAActctc                                                                                    <  1:724705/149‑1 (MQ=255)
                                                                aCGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTTAACTGAGGTCTTTTTTTGTGCCTGAAATCCGGTGTTTATCTATCTTTCCCGCTATATTAActctct                                                                                   >  2:188330/1‑149 (MQ=255)
                                                                aCGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAActctct                                                                                   >  3:335578/1‑149 (MQ=255)
                                                                  gCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCa                                                                                 >  3:41461/1‑149 (MQ=255)
                                                                    aGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAg                                                                               <  4:602199/149‑1 (MQ=255)
                                                                    aGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAg                                                                               <  1:728352/149‑1 (MQ=255)
                                                                     gTTGATTCAAAATCAACCGTAGAAATACGTTCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAgg                                                                              >  3:600119/1‑149 (MQ=255)
                                                                      ttGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGAATGTAAATAGACCTCAACGGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGt                                                                             >  4:370359/1‑149 (MQ=255)
                                                                       tGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTc                                                                            >  4:483991/1‑149 (MQ=255)
                                                                        gATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCa                                                                                                                      >  1:341449/1‑106 (MQ=255)
                                                                        gATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCa                                                                                                                      <  2:341449/106‑1 (MQ=255)
                                                                         aTTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCaa                                                                          >  2:698404/1‑149 (MQ=255)
                                                                          ttCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGt                                                                             <  4:669004/145‑1 (MQ=255)
                                                                          ttCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGt                                                                             >  3:669004/1‑145 (MQ=255)
                                                                          ttCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAAc                                                                         >  1:18650/1‑149 (MQ=255)
                                                                           tCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAAcc                                                                        <  1:459416/149‑1 (MQ=255)
                                                                            cAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCg                                                                       <  3:680549/149‑1 (MQ=255)
                                                                             aaaaTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGa                                                                      >  1:340589/1‑149 (MQ=255)
                                                                              aaaTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGat                                                                     <  4:364057/149‑1 (MQ=255)
                                                                                aTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGatat                                                                   >  1:157319/1‑149 (MQ=255)
                                                                                    aCCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAAc                                                               <  1:640208/149‑1 (MQ=255)
                                                                                    aCCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAAc                                                               >  3:25361/1‑149 (MQ=255)
                                                                                      cGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGt                                                             <  4:683748/149‑1 (MQ=255)
                                                                                       gTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTa                                                            >  3:317289/1‑149 (MQ=255)
                                                                                       gTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTa                                                            <  1:698404/149‑1 (MQ=255)
                                                                                          gAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAGTAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACAt                                                         >  1:500415/1‑149 (MQ=255)
                                                                                          gAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACAt                                                         <  2:157319/149‑1 (MQ=255)
                                                                                           aaaTACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATc                                                        >  3:649824/1‑149 (MQ=255)
                                                                                              tACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTAc                                                     >  2:620036/1‑149 (MQ=255)
                                                                                              tACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTAc                                                     <  4:401331/149‑1 (MQ=255)
                                                                                                 gTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCaa                                                  <  3:98515/149‑1 (MQ=255)
                                                                                                      ggTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACatat                                             <  2:660659/149‑1 (MQ=255)
                                                                                                       gTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATg                                            >  1:549258/1‑149 (MQ=255)
                                                                                                       gTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATg                                            >  3:215380/1‑149 (MQ=255)
                                                                                                         tCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGtt                                          <  1:249716/149‑1 (MQ=255)
                                                                                                          cGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTg                                         <  4:460095/149‑1 (MQ=255)
                                                                                                          cGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTg                                         <  4:707172/149‑1 (MQ=255)
                                                                                                             gTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTa                                      <  3:670040/149‑1 (MQ=255)
                                                                                                             gTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTa                                      <  2:157379/149‑1 (MQ=255)
                                                                                                                 ggCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAga                                  >  4:261033/1‑149 (MQ=255)
                                                                                                                 ggCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAga                                  >  1:516557/1‑149 (MQ=255)
                                                                                                                 ggCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAga                                  <  2:340589/149‑1 (MQ=255)
                                                                                                                 gcccTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAga                                  <  1:644770/147‑1 (MQ=255)
                                                                                                                     ttCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAgatgat                              <  2:110056/149‑1 (MQ=255)
                                                                                                                     ttCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAgatgat                              <  2:195929/149‑1 (MQ=255)
                                                                                                                     ttCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAgatgat                              <  3:119531/149‑1 (MQ=255)
                                                                                                                     ttCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAgatgat                              <  2:775536/149‑1 (MQ=255)
                                                                                                                        ggCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGt                           <  2:228758/149‑1 (MQ=255)
                                                                                                                          cACCAAAAGTATGTAACTAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTATCTCTCTCAAGGTCAACCGATATCAACGTACATCTCCCAACATATGTTGGTACAGATGATGGTAt                         <  1:703769/149‑1 (MQ=255)
                                                                                                                          cACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTAt                         >  3:555864/1‑149 (MQ=255)
                                                                                                                           aCCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTAtt                        >  4:35802/1‑149 (MQ=255)
                                                                                                                              aaaaGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTcc                     >  4:522813/1‑149 (MQ=255)
                                                                                                                              aaaaGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTcc                     >  4:274764/1‑149 (MQ=255)
                                                                                                                                aaGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTTTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATCGTATTTCCgg                   <  3:632021/149‑1 (MQ=255)
                                                                                                                                     tGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGa              >  2:521976/1‑149 (MQ=255)
                                                                                                                                        aaaaAGACCTCAGTTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATaa           >  4:426601/1‑149 (MQ=255)
                                                                                                                                        aaaaAGACCTCAGTTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATaa           >  1:733405/1‑149 (MQ=255)
                                                                                                                                        aaaaAGACCTCAGTTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATaa           >  2:558437/1‑149 (MQ=255)
                                                                                                                                            aGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCt       <  3:627743/149‑1 (MQ=255)
                                                                                                                                             gACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCtt      <  4:615720/149‑1 (MQ=255)
                                                                                                                                             gACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCtt      <  4:81792/149‑1 (MQ=255)
                                                                                                                                                cTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTa   <  4:600119/149‑1 (MQ=255)
                                                                                                                                                cTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTa   >  2:486226/1‑149 (MQ=255)
                                                                                                                                                 tCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTAc  <  2:503067/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTAC  >  NC_000913/4496363‑4496656

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: