Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,496,513 C→T 6.8% intergenic (+23/‑163) 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,5130CT 6.8% 154.9 / 2.3 59intergenic (+23/‑163)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 C (27/28);  new base T (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

GACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTA  >  NC_000913/4496365‑4496656
                                                                                                                                                    |                                                                                                                                               
gACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAAc                                                                                                                                                 >  1:464546/1‑149 (MQ=255)
 tcgggTACAATCCACAACGTCTGCCGCATACCCTTTGCGGGCCGAAGGGGCGAAAGGGGTAGACGCAGTTGATTCAAAATCACCCGTAGAAATACGTGCCGGTTCGGGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACt                                                                                                                                                <  2:34210/146‑1 (MQ=255)
 aCAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACt                                                                                                                                                >  1:521465/1‑149 (MQ=255)
   aGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGGATGTAAATAGACCTCAGCTGa                                                                                                                                              >  1:230004/1‑149 (MQ=255)
    ggTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCTTTCGGCATCAAAAGTATGTAAATAGACCTCAACTGAg                                                                                                                                             <  2:654516/149‑1 (MQ=255)
    ggTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAg                                                                                                                                             <  1:332387/149‑1 (MQ=255)
     gTATAATCAACAACGATTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTGCATGCAAAATCAACCGTAGAAATACATGCAGGTTAGAGTACGGCCTTCGGCACCGAAAGTATGAAAAGAAAACTCAACTGAgg                                                                                                                                            >  2:603803/1‑149 (MQ=255)
      tataATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAggc                                                                                                                                           >  1:323285/1‑148 (MQ=255)
      tataATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGACTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGt                                                                                                                                           <  1:154590/149‑1 (MQ=255)
          aTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCTTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttt                                                                                                                                       <  2:530159/149‑1 (MQ=255)
           tCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTCtatt                                                                                                                                      <  2:644093/149‑1 (MQ=255)
           tCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCtttt                                                                                                                                      <  2:564553/149‑1 (MQ=255)
           tCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAAAAGACCTCAGTTGAGGTCtatt                                                                                                                                      <  1:459508/149‑1 (MQ=255)
              acaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttttttt                                                                                                                                   <  2:99025/149‑1 (MQ=255)
                 aCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAAGCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATATACCCCAACTGAGGGCTTTTTTTATg                                                                                                                                >  1:481490/1‑149 (MQ=255)
                       tCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGaa                                                                                                                          >  2:303566/1‑149 (MQ=255)
                       tCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACCGAGGTCTTTTTTTATGCCTGaa                                                                                                                          >  2:303652/1‑149 (MQ=255)
                            aTACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCa                                                                                                                     <  1:274850/149‑1 (MQ=255)
                                atcttacgTGACGCAGTGCCGAAACCAGTAGACGCAGTTGATTCCAACCCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGTGGTCTTTTTTTATGCCTGAAATCCAgtgt                                                                                                                 <  1:550220/142‑1 (MQ=255)
                                      aGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATCCGTGCCGGTTCGATTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTAATGCCTGAATTCCAGTTTTtatca                                                                                                           >  1:360056/1‑148 (MQ=255)
                                      aGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTtatct                                                                                                           >  1:245749/1‑149 (MQ=255)
                                           cGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCtt                                                                                                      <  2:36939/149‑1 (MQ=255)
                                                tGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGTCTGAAATCCAGTGTTTATCTATCTTTCCCg                                                                                                 <  1:408328/149‑1 (MQ=255)
                                                 ggCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGc                                                                                                >  1:121300/1‑149 (MQ=255)
                                                   cGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCta                                                                                              >  2:698398/1‑149 (MQ=255)
                                                   cGAAATCGGCAGACGCAGTTGATTCAAAATCCACCGTAGAAATACGTCCCGGTTCGAGTCCGGCCTTCGGCACCCAAAGCAGGCCAATCGACCTCACCTCAGATCTTTTTTTATGCCAGAACTCCAATGTTTACCATCCTTTCCCGccg                                                                                              >  2:355468/1‑147 (MQ=255)
                                                    gAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGCAATACGTACCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTGTGCCTGAAATCCAGTGGGTATCTACCTTTTCCCCtat                                                                                             >  2:585071/1‑149 (MQ=255)
                                                           gTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCGGCTATATTAAct                                                                                      >  1:154715/1‑149 (MQ=255)
                                                            tAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAActc                                                                                     >  2:180403/1‑149 (MQ=255)
                                                                cGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTGTGTAAATAGACCTCAACTGAGGTCTTTTTTTCTGCCTGAAATCCGGTGTTTATCTATCTTTCCCGCTATCTTACCtctctc                                                                                 >  1:134731/1‑149 (MQ=255)
                                                                cGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCCCAACTGCGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACtctctc                                                                                 <  1:165714/149‑1 (MQ=255)
                                                                  cAGTTGATTCAAAATCAACCGTAAAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGATGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCaa                                                                               <  1:233937/149‑1 (MQ=255)
                                                                   aGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTGTATCTATCTCTCCCGCTATATTAACTCTCTCAAg                                                                              <  1:396332/149‑1 (MQ=255)
                                                                        aTTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTGTATTAACTCCTTCAAGGTCaa                                                                         >  1:489500/1‑149 (MQ=255)
                                                                         ttCCAAATCAACAGTCGAACTACATCCAGGTTCGAATCCGGCCATACGCACCAAAAGTATGTGAATAGACCTCAACTGAGGTCTTATTTTATGCCTGCTTTCCAGGGTTTTTCTATCTTTCCCGCTATATTACCCCTATCAAGGTCAAc                                                                        >  2:247169/1‑149 (MQ=255)
                                                                         ttCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAAc                                                                        >  1:338682/1‑149 (MQ=255)
                                                                               aTCAACCGTAGAAATACGTGCCTGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTCAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGatat                                                                  <  2:581687/149‑1 (MQ=255)
                                                                                        aGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGGACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACa                                                         >  1:571698/1‑149 (MQ=255)
                                                                                            aTACGTGACGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTa                                                     <  1:400254/149‑1 (MQ=255)
                                                                                                       ttCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCGGTCGTGCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGt                                          <  1:150443/149‑1 (MQ=255)
                                                                                                          gAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTgg                                       >  1:607961/1‑149 (MQ=255)
                                                                                                             tCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAGGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTAc                                    >  1:673102/1‑149 (MQ=255)
                                                                                                                  acTTCTGCACCAAAAGTATGTAAATAGACCGCAACTTAGGGCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAgatg                               <  2:101252/148‑1 (MQ=255)
                                                                                                                    ctgggCATCAAAAGTATGTAAACAGACCCTAACTGAGGTCCGTTTATCCGCCTTAAATGCAGTGTTTATCTACCTTGCCCGCTATATTGTCTCTGTCAAGGTCAACCGATATCAACGTACATCTACCATCATATGTTGGTACAggtgat                             <  2:114341/146‑1 (MQ=255)
                                                                                                                      cGGCACCAAAATGATGTAAAAAGCCCTCAACGTAAGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTCCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATgg                           <  1:686189/149‑1 (MQ=255)
                                                                                                                      cGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATgg                           <  2:338682/149‑1 (MQ=255)
                                                                                                                      cGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATgg                           <  1:166860/149‑1 (MQ=255)
                                                                                                                      aggCACCCATAGTATGTAAATACACCTCAACTTAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCCTTCCCGCTATTCAACCTCTCTCAAGGTCACCCGATATCAACGTTCATCTACCAACATATGTTGGTACAGATGATgg                           <  2:616160/148‑1 (MQ=255)
                                                                                                                       ggCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCAGCTATATTAACTCTCTCAAGGTCAACCCATATCAACGTACATCTACCAACATCTGTTGGTACAGACGACGGt                          >  2:120660/1‑149 (MQ=255)
                                                                                                                                 gTATGTAAATAGGCCTCAACTGAGGGCTTTTTTTATGCCTGAAAGCCAGTGTTTATCTATCTTTCCCGCTATATTAACGCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGtt                <  2:25938/149‑1 (MQ=255)
                                                                                                                                    tGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGa             <  2:127626/149‑1 (MQ=255)
                                                                                                                                     aaaaaaaGACCTCAGTTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGAt            >  2:499797/7‑149 (MQ=255)
                                                                                                                                     aaaaaaaGACCTCAGTTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGAt            >  2:499836/7‑149 (MQ=255)
                                                                                                                                       aaaTAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATaa          <  2:141315/149‑1 (MQ=255)
                                                                                                                                         aTAGACCTCAACTGAGGTCTTTTTTTATTCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATGAACTCTCTGAAGGTCAACCGATGTCAACGTACATCTAACAACATATGTTGGTACAGATACTGGTGTGTCCGGTTCGCTaagg        >  2:314221/1‑147 (MQ=255)
                                                                                                                                            gACCTCAACGGAGGTCTTTTTTTATGCCTGAAATCCAGTTTTTATCTATCTTGCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCtt     <  1:39469/149‑1 (MQ=255)
                                                                                                                                              ccTCAACTGAGGTCTTTTTTTTTTTTTTTAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTAAAGGTCAACCGATCTCAACGTACATCTACCAACATATGTTGATTCTTATGATTGTATTTCCGGTTCGATAATGCTTGt   >  2:385152/1‑149 (MQ=255)
                                                                                                                                              ccTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATTAGTTGGTCAAGATGATGGGATTTCCCGTTCGATAATACTTGt   >  2:121215/1‑149 (MQ=255)
                                                                                                                                               cTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAAACCTTGTa  >  1:295163/1‑149 (MQ=255)
                                                                                                                                               cTCAACCGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTa  <  1:627332/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGCTTGTA  >  NC_000913/4496365‑4496656

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