Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,496,511 A→G 5.5% 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.5% 246.9 / 3.4 73intergenic (+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 (38/31);  new base G (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

GGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCA  >  NC_000913/4496368‑4496658
                                                                                                                                               |                                                                                                                                                   
ggTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAg                                                                                                                                                <  2:173934/149‑1 (MQ=255)
    taATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCt                                                                                                                                            <  4:511638/149‑1 (MQ=255)
      aTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACt                                                                                                                                                   >  4:1461544/1‑140 (MQ=255)
      aTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACt                                                                                                                                                   <  3:1461544/140‑1 (MQ=255)
        ccACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttttt                                                                                                                                        >  3:1084775/1‑149 (MQ=255)
          acaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttttttt                                                                                                                                      >  4:39196/1‑149 (MQ=255)
            aaCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTAt                                                                                                                                    >  2:563768/1‑149 (MQ=255)
            aaCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTAt                                                                                                                                    >  2:798756/1‑149 (MQ=255)
                      gCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATc                                                                                                                          <  3:39196/149‑1 (MQ=255)
                        aTACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTGTATGCCTGAAATCCa                                                                                                                        <  4:1084775/149‑1 (MQ=255)
                                  aGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTtatct                                                                                                              <  1:798756/149‑1 (MQ=255)
                                                          gACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCa                                                                                                                        <  4:521208/115‑1 (MQ=255)
                                                          gACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCa                                                                                                                        >  3:521208/1‑115 (MQ=255)
                                                            cGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATg                                                                                                                                   >  1:1079043/1‑102 (MQ=255)
                                                            cGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATg                                                                                                                                   <  2:1079043/102‑1 (MQ=255)
                                                             gCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTCt                                                                                                                                            >  1:930519/1‑92 (MQ=255)
                                                             gCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTCt                                                                                                                                            <  2:930519/92‑1 (MQ=255)
                                                               aGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAg                                                                                 <  1:563768/149‑1 (MQ=255)
                                                                     ttCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAAc                                                                           >  3:1407295/1‑149 (MQ=255)
                                                                      tCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGaaa                                                                                                                            <  4:1245691/99‑1 (MQ=255)
                                                                      tCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGaaa                                                                                                                            >  3:1245691/1‑99 (MQ=255)
                                                                        aaaaTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACt                                                                                                                                                   <  1:1209567/74‑1 (MQ=255)
                                                                        aaaaTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACt                                                                                                                                                   >  2:1209567/1‑74 (MQ=255)
                                                                        aaaaTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGa                                                                        >  2:1021432/1‑149 (MQ=255)
                                                                           aTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTcc                                                                                                      <  4:474541/116‑1 (MQ=255)
                                                                           aTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTcc                                                                                                      >  3:474541/1‑116 (MQ=255)
                                                                            tCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATc                                                                    >  2:1046396/1‑149 (MQ=255)
                                                                              aaCCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAAc                                                                           <  1:648764/140‑1 (MQ=255)
                                                                              aaCCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAAc                                                                           >  2:648764/1‑140 (MQ=255)
                                                                              aaCCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCaa                                                                  <  1:1046396/149‑1 (MQ=255)
                                                                                   tAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAg                                                                                                                                                >  2:275527/1‑66 (MQ=255)
                                                                                   tAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAg                                                                                                                                                <  1:275527/66‑1 (MQ=255)
                                                                                   tAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTAc                                                             >  4:1433678/1‑149 (MQ=255)
                                                                                    aGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAAc                                                                                          <  3:667768/119‑1 (MQ=255)
                                                                                    aGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAAc                                                                                          >  4:667768/1‑119 (MQ=255)
                                                                                                cGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACata                                                <  3:1433678/149‑1 (MQ=255)
                                                                                                  gTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCt                                                                                                  >  1:621255/1‑97 (MQ=255)
                                                                                                  gTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCt                                                                                                  <  2:621255/97‑1 (MQ=255)
                                                                                                   ttCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCtttt                                                                                                                                         <  1:162376/57‑1 (MQ=255)
                                                                                                   ttCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCtttt                                                                                                                                         >  2:162376/1‑57 (MQ=255)
                                                                                                       aGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCaa                                                                  <  2:235389/124‑1 (MQ=255)
                                                                                                       aGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCaa                                                                  >  1:235389/1‑124 (MQ=255)
                                                                                                           cGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTgg                                          >  4:1270067/1‑144 (MQ=255)
                                                                                                           cGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTgg                                          <  3:1270067/144‑1 (MQ=255)
                                                                                                            ggCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAga                                    <  1:1021432/149‑1 (MQ=255)
                                                                                                                ttCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAAc                                                                 >  4:17804/1‑116 (MQ=255)
                                                                                                                ttCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAAc                                                                 <  3:17804/116‑1 (MQ=255)
                                                                                                                ttCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAgatgat                                >  3:1034656/1‑149 (MQ=255)
                                                                                                                        cAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAAc                                                                                          <  3:1100793/83‑1 (MQ=255)
                                                                                                                        cAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAAc                                                                                          >  4:1100793/1‑83 (MQ=255)
                                                                                                                        cAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTAc                                                             >  3:442952/1‑112 (MQ=255)
                                                                                                                        cAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTAc                                                             <  4:442952/112‑1 (MQ=255)
                                                                                                                        cAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTc                        >  1:977988/1‑149 (MQ=255)
                                                                                                                         aaaaGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTcc                       <  4:1407295/149‑1 (MQ=255)
                                                                                                                                tGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAActct                                                                                       <  4:889248/78‑1 (MQ=255)
                                                                                                                                tGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAActct                                                                                       >  3:889248/1‑78 (MQ=255)
                                                                                                                                   aaaaAGACCTCAGTTGAGGTCTATTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTAc                                                             <  2:683296/101‑1 (MQ=39)
                                                                                                                                   aaaaAGACCTCAGTTGAGGTCTATTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTAc                                                             >  1:683296/1‑101 (MQ=39)
                                                                                                                                   aaaTAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATaa             >  2:1229929/1‑149 (MQ=255)
                                                                                                                                    aaTAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCta                                                                                                 >  2:843/1‑64 (MQ=255)
                                                                                                                                    aaTAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCta                                                                                                 <  1:843/64‑1 (MQ=255)
                                                                                                                                     aTAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTc                  >  2:507601/1‑142 (MQ=255)
                                                                                                                                     aTAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTc                  <  1:507601/142‑1 (MQ=255)
                                                                                                                                      tAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGc          >  4:500586/1‑149 (MQ=255)
                                                                                                                                          ccTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGt      >  4:1036153/1‑149 (MQ=255)
                                                                                                                                           cTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCgg                     <  1:246994/133‑1 (MQ=255)
                                                                                                                                           cTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCgg                     >  2:246994/1‑133 (MQ=255)
                                                                                                                                           cTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTa     >  3:336773/1‑149 (MQ=255)
                                                                                                                                           cTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTa     >  4:1393163/1‑149 (MQ=255)
                                                                                                                                           cTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTa     >  4:1391831/1‑149 (MQ=255)
                                                                                                                                            tCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAg                                                                                 >  1:928692/1‑72 (MQ=255)
                                                                                                                                            tCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAg                                                                                 <  2:928692/72‑1 (MQ=255)
                                                                                                                                             cAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTAcc   >  1:384317/1‑149 (MQ=255)
                                                                                                                                              aaCTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCa  >  4:478277/1‑149 (MQ=255)
                                                                                                                                              aaCTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCa  >  1:320178/1‑149 (MQ=255)
                                                                                                                                               aCTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCaa                                                    <  4:496312/98‑1 (MQ=255)
                                                                                                                                               aCTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCaa                                                    >  3:496312/1‑98 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTACCA  >  NC_000913/4496368‑4496658

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑