Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,519,313 C→A 8.2% intergenic (+89/+25) insB1 → / ← yjhU IS1 transposase B;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,519,3130CA 8.2% 202.6 / 3.5 87intergenic (+89/+25)insB1/yjhUIS1 transposase B;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element
Reads supporting (aligned to +/- strand):  ref base C (42/37);  new base A (6/2);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAA  >  NC_000913/4519166‑4519459
                                                                                                                                                    |                                                                                                                                                  
cGGTGGAGTTTCATGACAAGGTCAACGGGCATTATCTGAACATAAAACACAATCAGTAAGTTGGAGTGATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCt                                                                                                                                                   <  2:1027309/149‑1 (MQ=255)
  gTGGAGTTGCATGACAAGGTCATCGGGCATTAACTGAACATAAAACACTATCAGTAAGTTGGAGTAATTACCGCCCATTTTTATTTCATACATTGTGGGTTTTGTTCTTATTATCGCCGATAAGCAATAAAATCCTGCACAATATATAc                                                                                                                                                 >  2:1027735/1‑149 (MQ=255)
      aGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATgg                                                                                                                                             <  1:76811/149‑1 (MQ=255)
          gCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCa                                                                                                                                         >  1:182661/1‑149 (MQ=255)
           cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAg                                                                                                                                        >  2:179329/1‑149 (MQ=255)
           aaTGACAAGGTCATCGGGCATTATATGACCATAAAACACTAACAGTAAGTAGGAGTCATTACCGACCATGTTTATTTCATACATTTTGGGTATTGTTCTTAATATCGCCGCTAATCAATAAAATCCTGCCCCATCTAGATATGGGGCAg                                                                                                                                        <  2:261183/148‑1 (MQ=255)
                  aGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCa                                                                                                                                 <  1:157445/149‑1 (MQ=255)
                  aGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATAGCTACATCGGGCAGTTGTTCa                                                                                                                                 <  2:53091/149‑1 (MQ=255)
                    gTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCAtt                                                                                                                               >  2:761853/1‑149 (MQ=255)
                     gcATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTc                                                                                                                              <  1:797819/148‑1 (MQ=255)
                        tCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCttt                                                                                                                           <  2:377625/149‑1 (MQ=255)
                          gggCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTa                                                                                                                         <  1:308867/149‑1 (MQ=255)
                             gaTTATCTGAACATAAATAACTATCAGTAAGTTGTAGTCATTACCGACTATGTTTTTTTTTTAAATTTTGTGTATTGTTCTTATTATTTCAGATAATGAATTAAATACTGCCCCATATATACATGGGGCAGTTGTTCATTCTTTTAgtg                                                                                                                      <  2:558772/148‑1 (MQ=255)
                             cATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAgtg                                                                                                                      >  1:546986/1‑149 (MQ=255)
                               ttATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCATTAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGtgtg                                                                                                                    >  1:826023/1‑149 (MQ=255)
                                  tCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTa                                                                                                                 >  1:981229/1‑149 (MQ=255)
                                     gAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAAtt                                                                                                              >  2:25287/1‑149 (MQ=255)
                                      aaCATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTc                                                                                                             <  2:546986/149‑1 (MQ=255)
                                        cATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTcac                                                                                                           >  2:953333/1‑149 (MQ=255)
                                        cATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTATATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCAATATCTACATGGGGCATTTGTTCATTCTTTTAGTGTGGTAATTcac                                                                                                           >  2:657415/1‑149 (MQ=255)
                                         aTAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTcaca                                                                                                          >  2:660001/1‑149 (MQ=255)
                                         aTAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTCGTAATTcaca                                                                                                          >  2:865156/1‑149 (MQ=255)
                                            aaaCACTATCAGTAAGTGGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGc                                                                                                       <  1:904850/149‑1 (MQ=255)
                                             aaCACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGcc                                                                                                      >  1:820014/1‑149 (MQ=255)
                                             aaCACTATCAGTAAGTTGGAGTCATTACCGACAATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGACCCATATCTACATGGGGCAGTTTTTCATTCTTTTAGTTTGGTAATTCACAcccc                                                                                                      >  1:534072/1‑149 (MQ=255)
                                                acTATCAGTAAGTTGGAGTCATTACCGAACATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGGCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGc                                                                                                   <  2:329556/149‑1 (MQ=255)
                                                 cTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCa                                                                                                  <  1:325558/149‑1 (MQ=255)
                                                  tATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCACATATCTACATGTTGCAGTTCTTCATTCTTTTAGTGTGGTAATTCACACGccctgaa                                                                                                 >  2:429867/1‑144 (MQ=255)
                                                       gTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCTTATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAAct                                                                                            >  1:464901/1‑149 (MQ=255)
                                                        tAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTGTTTTAGTGTGGTAATTCACACGCCAGCAAAAActc                                                                                           >  2:29855/1‑149 (MQ=255)
                                                        tAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAActc                                                                                           >  2:959137/1‑149 (MQ=255)
                                                        tAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAActa                                                                                           >  1:653747/1‑148 (MQ=255)
                                                        tAAGTTGGAGTAATTACCGACCATGTTTTTTTCATAGATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAActc                                                                                           <  2:182661/149‑1 (MQ=255)
                                                          aGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTg                                                                                         >  2:143412/1‑149 (MQ=255)
                                                             tGGAGTCAGTACCGACCATGTTTAATTCATACAATCTGTGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCAGATCGACATGCTCCACTAGTTAATTCTTTTACTGTGGTAATTCACACGCCAGAAAAAACTCTCCCg                                                                                      >  2:182340/1‑149 (MQ=255)
                                                              ggAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGt                                                                                     >  1:659351/1‑149 (MQ=255)
                                                                 gTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTcc                                                                                  >  2:40762/1‑149 (MQ=255)
                                                                   cATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCtt                                                                                <  1:270805/149‑1 (MQ=255)
                                                                     ttACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTTATTAACACGCCAGCAAAAGCTCTGCCGTTCCTtca                                                                              >  1:80690/1‑149 (MQ=255)
                                                                       aCCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtcatc                                                                            <  1:111625/149‑1 (MQ=255)
                                                                        ccGACCATGTTTATTTCATACATAGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCAGCCCAATATATACATGGGGCAGTTGTTCATTCTTTAAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCGTtcatca                                                                           >  2:279850/1‑149 (MQ=255)
                                                                         cGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCGATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGACGTTCCTTCATcaa                                                                          >  1:445061/1‑149 (MQ=255)
                                                                         cGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaa                                                                          >  1:1014964/1‑149 (MQ=255)
                                                                            ccATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAAACCTGCCCCATATCTACATGGGGCAGTTGTTCATTCATTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaacaa                                                                       >  2:171174/1‑149 (MQ=255)
                                                                              atttttATTTCATAAATTGTGGGTATTGTTATTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATc                                                                     <  2:740697/146‑1 (MQ=255)
                                                                                gTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAg                                                                   <  1:909148/149‑1 (MQ=255)
                                                                                  ttATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGt                                                                 <  2:630646/149‑1 (MQ=255)
                                                                                      ttCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGt                                                             <  2:1027197/149‑1 (MQ=255)
                                                                                         aTACATTGTGGGTATTGTTCTTAATATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGAc                                                          <  2:158470/149‑1 (MQ=255)
                                                                                             aTTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACatat                                                      >  2:554351/1‑149 (MQ=255)
                                                                                                   ggTATTGTTCTTATTATCGCCGCTAATTAATAAAATCCTGCCCCTTTTCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCTCACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTccc                                                <  1:850727/149‑1 (MQ=255)
                                                                                                    gTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCca                                               <  1:330558/149‑1 (MQ=255)
                                                                                                     tATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCcag                                              >  1:338796/1‑149 (MQ=255)
                                                                                                       ttGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCcagca                                            >  1:763109/1‑149 (MQ=255)
                                                                                                        tGTTCTTATTATCGTCGCAAATCAATAAAATCCTGCCCCATATCTCCATGGGGCAGTTGATCATTCTTTTAGTGTGGTAAATCACACGCCAGCAAAAACTCTGCCGTTCCATCATCAACAATCAGGTCCGTGACATATCCTCCcagcag                                           <  2:241714/149‑1 (MQ=255)
                                                                                                        tGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCcagcag                                           <  1:906719/149‑1 (MQ=255)
                                                                                                         gTTCTTATTTTCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAgg                                          >  2:524360/1‑149 (MQ=255)
                                                                                                         gTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCGTTAATCAACAATCAGGTCCGTGACATATCCTCCCAGCAgg                                          >  1:363020/1‑149 (MQ=255)
                                                                                                              tcttatCGCCGCAAATAATATAAATCCTGCCCCTTCTCTCCATGGTGCAGTGTTTCATTCTATTAGTGTGGTATTTCACACGCAAGCCAAAACTCTGCCGGTCCTTCATCAACAATCAGGTCCGTGACAAATCCTCCCAGCAGGGCAcc                                     <  2:692216/147‑1 (MQ=255)
                                                                                                              tattatCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCATTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCAcc                                     >  1:463814/1‑149 (MQ=255)
                                                                                                               attatCGCAGCTAATCAATAAAATCCTTCCCCATATCTACATGTGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCAccc                                    <  1:4221/149‑1 (MQ=255)
                                                                                                                  atCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTTGTAATTCACACGCCAGCAAAACCTCTGCCGTTCCTTCATCAACAATCATGTCCGTGCCATATCCTCCCAGCAGGGCACCCAAc                                 >  2:829941/1‑149 (MQ=255)
                                                                                                                   acgccgcTAATCAATAAAATCCTGCCCCATAACTACATGGGGCACTTGTTCATTCTTTTAGTGAGGTAATTCTCACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCTGATCCGAGACATATCCACCCAGCAGGGCACCCATCg                                >  1:441808/2‑149 (MQ=255)
                                                                                                                       cgcTCATCAATACAATCCATCCCCATAGCAACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTACTTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACACTCAGGTCCGTGACAAATCCTCCCAGCAGGGGACCCAACGTTGc                            <  1:207393/149‑1 (MQ=255)
                                                                                                                        tcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCg                           <  1:40762/148‑1 (MQ=255)
                                                                                                                        gcTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCCCCCAACGTTGCg                           >  2:148510/1‑149 (MQ=255)
                                                                                                                             taaaTAAAATCCTGCCCCATATCTCCATCGGGCAGTTGTTCATTCTTTTAGTGGGGTAATTCACACGCCAGCAAAAACTCTGCCGTGCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATa                      <  2:611175/147‑1 (MQ=255)
                                                                                                                             tCAATAAAATCCTGCCCCATATCTACTTGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAACACTCTGCCGTTCATTCATCAACAATAAGGTCCGGGACATATCCTCACAGCAGGGCACCCAACGTTGCGTCATa                      >  1:582576/1‑149 (MQ=255)
                                                                                                                               aaaaaaaTCCAGCCCCATATCTACATGGGTCAGATGATCAATCTTTTAGTGTGGTAAATCACACCCCAGAAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCAAGCAGGGAACCCCACGATGCGACATAGc                    >  1:587419/4‑149 (MQ=255)
                                                                                                                               aaTAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACATTCAGGTCCGTGACATATCCTCCCAGCAGGGCCCCCAACGTGGCGTCATAGc                    >  1:530500/1‑149 (MQ=255)
                                                                                                                               aaTAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACACTCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGc                    <  1:9907/149‑1 (MQ=255)
                                                                                                                               aaTAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGc                    <  1:865156/149‑1 (MQ=255)
                                                                                                                                  aaaaTCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCctc                 <  2:820014/149‑1 (MQ=255)
                                                                                                                                     aaTCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCtt               >  2:998393/1‑149 (MQ=255)
                                                                                                                                     aaTCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCtt               <  2:455480/149‑1 (MQ=255)
                                                                                                                                       tCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTc             <  1:393096/149‑1 (MQ=255)
                                                                                                                                       tCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTc             <  1:179329/149‑1 (MQ=255)
                                                                                                                                       tCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTc             >  1:999512/1‑149 (MQ=255)
                                                                                                                                         cTGCCCCATGTAGATATGGGGCAGTTTTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTTCCGTTCCTTAAACAACAAACAGGACCTTGAAATATACTCCCATCAGGGCACCCAACGTTGCCTCATAGCCTCTTTacc           >  1:246563/1‑149 (MQ=255)
                                                                                                                                         cTGCCCCATGTAGATATGGGGCAGTAGTTCATTCTTTTAGTCTGTTAATTAACACGCCAGAAAAAACTCTGCAGTTCCTTCATCAACAATCAGGTCAGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGGCATAGACTCTTTccc           >  1:386571/1‑149 (MQ=255)
                                                                                                                                          tGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCCACCAACGTTGCGTCATAGCCTCTTTcccc          >  1:329198/1‑149 (MQ=255)
                                                                                                                                          tGCCCCATGTAGATATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcccc          >  1:812067/1‑149 (MQ=255)
                                                                                                                                          tGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCCGCAAAAACTCTCGCGTTCCTTAATCAGCAATCAGGTCCGTGACATATCCTCCCAGCAGGGCGCCCAACGATGCGTCATAGCCTCTATcccc          >  1:495521/1‑149 (MQ=255)
                                                                                                                                            ccccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTcccccc        >  1:250272/1‑147 (MQ=255)
                                                                                                                                              ccATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCccacca      <  1:499083/149‑1 (MQ=255)
                                                                                                                                                  atCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCaa  <  1:657415/149‑1 (MQ=255)
                                                                                                                                                  atCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCaa  <  1:953333/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATAAAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAA  >  NC_000913/4519166‑4519459

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: