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

CGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAA  >  NC_000913/4519166‑4519459
                                                                                                                                                    |                                                                                                                                                  
cGGTGGAGTTTCATGACAAGGTCAACGGGCATTATCTGAACATAAAACACAATCAGTAAGTTGGAGTGATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCt                                                                                                                                                   <  2:1027309/149‑1 (MQ=255)
  gTGGAGTTGCATGACAAGGTCATCGGGCATTAACTGAACATAAAACACTATCAGTAAGTTGGAGTAATTACCGCCCATTTTTATTTCATACATTGTGGGTTTTGTTCTTATTATCGCCGATAAGCAATA‑AAATCCTGCACAATATATAc                                                                                                                                                 >  2:1027735/1‑149 (MQ=255)
      aGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATgg                                                                                                                                             <  1:76811/149‑1 (MQ=255)
          gCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCa                                                                                                                                         >  1:182661/1‑149 (MQ=255)
           cATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAg                                                                                                                                        >  2:179329/1‑149 (MQ=255)
           aaTGACAAGGTCATCGGGCATTATATGACCATAAAACACTAACAGTAAGTAGGAGTCATTACCGACCATGTTTATTTCATACATTTTGGGTATTGTTCTTAATATCGCCGCTAATCAATA‑AAATCCTGCCCCATCTAGATATGGGGCAg                                                                                                                                        <  2:261183/148‑1 (MQ=255)
                  aGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCa                                                                                                                                 <  1:157445/149‑1 (MQ=255)
                  aGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATAGCTACATCGGGCAGTTGTTCa                                                                                                                                 <  2:53091/149‑1 (MQ=255)
                    gTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCAtt                                                                                                                               >  2:761853/1‑149 (MQ=255)
                     gcATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTc                                                                                                                              <  1:797819/148‑1 (MQ=255)
                        tCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCttt                                                                                                                           <  2:377625/149‑1 (MQ=255)
                          gggCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTa                                                                                                                         <  1:308867/149‑1 (MQ=255)
                             gaTTATCTGAACATAAATAACTATCAGTAAGTTGTAGTCATTACCGACTATGTTTTTTTTTTAAATTTTGTGTATTGTTCTTATTATTTCAGATAATGAATT‑AAATACTGCCCCATATATACATGGGGCAGTTGTTCATTCTTTTAgtg                                                                                                                      <  2:558772/148‑1 (MQ=255)
                             cATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAgtg                                                                                                                      >  1:546986/1‑149 (MQ=255)
                               ttATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCATTA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGtgtg                                                                                                                    >  1:826023/1‑149 (MQ=255)
                                  tCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTa                                                                                                                 >  1:981229/1‑149 (MQ=255)
                                     gAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAAtt                                                                                                              >  2:25287/1‑149 (MQ=255)
                                      aaCATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTc                                                                                                             <  2:546986/149‑1 (MQ=255)
                                        cATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTcac                                                                                                           >  2:953333/1‑149 (MQ=255)
                                        cATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTATATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCAATATCTACATGGGGCATTTGTTCATTCTTTTAGTGTGGTAATTcac                                                                                                           >  2:657415/1‑149 (MQ=255)
                                         aTAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTcaca                                                                                                          >  2:660001/1‑149 (MQ=255)
                                         aTAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTCGTAATTcaca                                                                                                          >  2:865156/1‑149 (MQ=255)
                                            aaaCACTATCAGTAAGTGGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGc                                                                                                       <  1:904850/149‑1 (MQ=255)
                                             aaCACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGcc                                                                                                      >  1:820014/1‑149 (MQ=255)
                                             aaCACTATCAGTAAGTTGGAGTCATTACCGACAATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGACCCATATCTACATGGGGCAGTTTTTCATTCTTTTAGTTTGGTAATTCACAcccc                                                                                                      >  1:534072/1‑149 (MQ=255)
                                                acTATCAGTAAGTTGGAGTCATTACCGAACATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGGCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGc                                                                                                   <  2:329556/149‑1 (MQ=255)
                                                 cTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCa                                                                                                  <  1:325558/149‑1 (MQ=255)
                                                  tATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCACATATCTACATGTTGCAGTTCTTCATTCTTTTAGTGTGGTAATTCACACGccctgaa                                                                                                 >  2:429867/1‑144 (MQ=255)
                                                       gTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCTTATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAAct                                                                                            >  1:464901/1‑149 (MQ=255)
                                                        tAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTGTTTTAGTGTGGTAATTCACACGCCAGCAAAAActc                                                                                           >  2:29855/1‑149 (MQ=255)
                                                        tAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAActc                                                                                           >  2:959137/1‑149 (MQ=255)
                                                        tAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAActa                                                                                           >  1:653747/1‑148 (MQ=255)
                                                        tAAGTTGGAGTAATTACCGACCATGTTTTTTTCATAGATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAActc                                                                                           <  2:182661/149‑1 (MQ=255)
                                                          aGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTg                                                                                         >  2:143412/1‑149 (MQ=255)
                                                             tGGAGTCAGTACCGACCATGTTTAATTCATACAATCTGTGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCAGATCGACATGCTCCACTAGTTAATTCTTTTACTGTGGTAATTCACACGCCAGAAAAAACTCTCCCg                                                                                      >  2:182340/1‑149 (MQ=255)
                                                              ggAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGt                                                                                     >  1:659351/1‑149 (MQ=255)
                                                                 gTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTcc                                                                                  >  2:40762/1‑149 (MQ=255)
                                                                   cATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCtt                                                                                <  1:270805/149‑1 (MQ=255)
                                                                     ttACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTTATTAACACGCCAGCAAAAGCTCTGCCGTTCCTtca                                                                              >  1:80690/1‑149 (MQ=255)
                                                                       aCCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTtcatc                                                                            <  1:111625/149‑1 (MQ=255)
                                                                        ccGACCATGTTTATTTCATACATAGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCAGCCCAATATATACATGGGGCAGTTGTTCATTCTTTAAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCGTtcatca                                                                           >  2:279850/1‑149 (MQ=255)
                                                                         cGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCGATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGACGTTCCTTCATcaa                                                                          >  1:445061/1‑149 (MQ=255)
                                                                         cGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaa                                                                          >  1:1014964/1‑149 (MQ=255)
                                                                            ccATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAAACCTGCCCCATATCTACATGGGGCAGTTGTTCATTCATTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATcaacaa                                                                       >  2:171174/1‑149 (MQ=255)
                                                                              atttttATTTCATAAATTGTGGGTATTGTTATTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATc                                                                     <  2:740697/146‑1 (MQ=255)
                                                                                gTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAg                                                                   <  1:909148/149‑1 (MQ=255)
                                                                                  ttATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGt                                                                 <  2:630646/149‑1 (MQ=255)
                                                                                      ttCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGt                                                             <  2:1027197/149‑1 (MQ=255)
                                                                                         aTACATTGTGGGTATTGTTCTTAATATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGAc                                                          <  2:158470/149‑1 (MQ=255)
                                                                                             aTTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACatat                                                      >  2:554351/1‑149 (MQ=255)
                                                                                                   ggTATTGTTCTTATTATCGCCGCTAATTAATA‑AAATCCTGCCCCTTTTCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCTCACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTccc                                                <  1:850727/149‑1 (MQ=255)
                                                                                                    gTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCca                                               <  1:330558/149‑1 (MQ=255)
                                                                                                     tATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCcag                                              >  1:338796/1‑149 (MQ=255)
                                                                                                       ttGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCcagca                                            >  1:763109/1‑149 (MQ=255)
                                                                                                        tGTTCTTATTATCGTCGCAAATCAATA‑AAATCCTGCCCCATATCTCCATGGGGCAGTTGATCATTCTTTTAGTGTGGTAAATCACACGCCAGCAAAAACTCTGCCGTTCCATCATCAACAATCAGGTCCGTGACATATCCTCCcagcag                                           <  2:241714/149‑1 (MQ=255)
                                                                                                        tGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCcagcag                                           <  1:906719/149‑1 (MQ=255)
                                                                                                         gTTCTTATTTTCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAgg                                          >  2:524360/1‑149 (MQ=255)
                                                                                                         gTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCGTTAATCAACAATCAGGTCCGTGACATATCCTCCCAGCAgg                                          >  1:363020/1‑149 (MQ=255)
                                                                                                              tcttatCGCCGCAAAT‑AATATAAATCCTGCCCCTTCTCTCCATGGTGCAGTGTTTCATTCTATTAGTGTGGTATTTCACACGCAAGCCAAAACTCTGCCGGTCCTTCATCAACAATCAGGTCCGTGACAAATCCTCCCAGCAGGGCAcc                                     <  2:692216/147‑1 (MQ=255)
                                                                                                              tattatCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCATTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCAcc                                     >  1:463814/1‑149 (MQ=255)
                                                                                                               attatCGCAGCTAATCAATA‑AAATCCTTCCCCATATCTACATGTGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCAccc                                    <  1:4221/149‑1 (MQ=255)
                                                                                                                  atCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTTGTAATTCACACGCCAGCAAAACCTCTGCCGTTCCTTCATCAACAATCATGTCCGTGCCATATCCTCCCAGCAGGGCACCCAAc                                 >  2:829941/1‑149 (MQ=255)
                                                                                                                   acgccgcTAATCAATA‑AAATCCTGCCCCATAACTACATGGGGCACTTGTTCATTCTTTTAGTGAGGTAATTCTCACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCTGATCCGAGACATATCCACCCAGCAGGGCACCCATCg                                >  1:441808/2‑149 (MQ=255)
                                                                                                                       cgcTCATCAATA‑CAATCCATCCCCATAGCAACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTACTTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACACTCAGGTCCGTGACAAATCCTCCCAGCAGGGGACCCAACGTTGc                            <  1:207393/149‑1 (MQ=255)
                                                                                                                        tcTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCg                           <  1:40762/148‑1 (MQ=255)
                                                                                                                        gcTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCCCCCAACGTTGCg                           >  2:148510/1‑149 (MQ=255)
                                                                                                                             taaaTA‑AAATCCTGCCCCATATCTCCATCGGGCAGTTGTTCATTCTTTTAGTGGGGTAATTCACACGCCAGCAAAAACTCTGCCGTGCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATa                      <  2:611175/147‑1 (MQ=255)
                                                                                                                             tCAATA‑AAATCCTGCCCCATATCTACTTGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAACACTCTGCCGTTCATTCATCAACAATAAGGTCCGGGACATATCCTCACAGCAGGGCACCCAACGTTGCGTCATa                      >  1:582576/1‑149 (MQ=255)
                                                                                                                               aaaa‑aaaTCCAGCCCCATATCTACATGGGTCAGATGATCAATCTTTTAGTGTGGTAAATCACACCCCAGAAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCAAGCAGGGAACCCCACGATGCGACATAGc                    >  1:587419/4‑149 (MQ=255)
                                                                                                                               aaTA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACATTCAGGTCCGTGACATATCCTCCCAGCAGGGCCCCCAACGTGGCGTCATAGc                    >  1:530500/1‑149 (MQ=255)
                                                                                                                               aaTA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACACTCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGc                    <  1:9907/149‑1 (MQ=255)
                                                                                                                               aaTA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGc                    <  1:865156/149‑1 (MQ=255)
                                                                                                                                  a‑aaaTCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCctc                 <  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)
                                                                                                                                                    |                                                                                                                                                  
CGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTATTATCGCCGCTAATCAATA‑AAATCCTGCCCCATATCTACATGGGGCAGTTGTTCATTCTTTTAGTGTGGTAATTCACACGCCAGCAAAAACTCTGCCGTTCCTTCATCAACAATCAGGTCCGTGACATATCCTCCCAGCAGGGCACCCAACGTTGCGTCATAGCCTCTTTCCCCACCAGCAA  >  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: ‑