Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,496,512 C→T 6.2% intergenic (+23/‑238) leuX → / → intB tRNA‑Leu/KpLE2 phage‑like element; putative integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,496,5120CT 6.2% 179.2 / 2.8 65intergenic (+23/‑238)leuX/intBtRNA‑Leu/KpLE2 phage‑like element; putative integrase
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base T (2/2);  total (35/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 = 5.20e-01

ACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGC  >  NC_000913/4496378‑4496650
                                                                                                                                       |                                                                                                                                          
acaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGCGGTCt                                                                                                                                    >  1:342670/1‑143 (MQ=255)
  aaCGTTTTCCGCATCCCTCTTCAGTGCCGAAGTGGCGAACCCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAAGACGTGCCAGCTCCAGTCCAGCCTTCGGCCCCAACAGTTTGTAAATAGACCTCAACTGTGGTCttt                                                                                                                                  >  1:428794/1‑143 (MQ=255)
  aaCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGGTCGAGGCCGGCCTTCGGCACCCACAGTATGTAAATAGACCTCAACTGAGGTCttt                                                                                                                                  >  2:629582/1‑143 (MQ=255)
   aCGTTTTCCGCATGCCGCTTCAGTGCCGAAGGGGCGAAATCGGTAGACGCAGTGGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGTTGAGGTCtatt                                                                                                                                 <  1:614160/143‑1 (MQ=255)
   aCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGCCCGGCCGTCAGCACCAAACGCATGTAAATAGACCTCAACTGAggcccata                                                                                                                                 >  1:907653/1‑137 (MQ=255)
    cGTTTTCCGCATGCCTCTTCAGTGCCGTAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAGCTGAGGTCttttt                                                                                                                                <  2:559046/143‑1 (MQ=255)
      ttttCCGCATACCTCTTCAGTGCCGAAGTGGTGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttttttt                                                                                                                              <  1:426525/143‑1 (MQ=255)
      ttttCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGACATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTCAGGTCttttttt                                                                                                                              >  2:44661/1‑143 (MQ=255)
      ttttCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCttttttt                                                                                                                              <  2:756630/143‑1 (MQ=255)
         tCCACACACCTTTTCAGTGCCGAAGTGGCGAAACCGGGAGACGCGGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAAAAGACCTCAGTTGAGGTCtatttacata                                                                                                                           <  2:468718/143‑6 (MQ=255)
          ccGCCTGCCTCTCCAGCGCCGTCGTGGGGACATCGGTCGCCGCAGTTGATTCAAAATCACCCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGc                                                                                                                          <  1:80043/143‑1 (MQ=255)
            gCATACCTCTTCAGTGGCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCCAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCt                                                                                                                        <  2:715335/143‑1 (MQ=255)
               tACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGaa                                                                                                                     <  2:868661/143‑1 (MQ=255)
                        aGTGCCGAATTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAgtgt                                                                                                            <  1:317194/143‑1 (MQ=255)
                         gTGCCGAGGGGGCGAGAGCGGTGGGCGCAGTTGGTTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGtt                                                                                                           <  2:146942/143‑1 (MQ=255)
                          tGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGttt                                                                                                          >  1:98624/1‑143 (MQ=255)
                          tGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGttt                                                                                                          <  1:590366/143‑1 (MQ=255)
                          tGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTACGGCACCAACAGTATGGAAATAGACCTCCACCGAGGTCTTGTTTTATGCCTGAACTCCAgtcgtt                                                                                                          >  2:441423/1‑139 (MQ=255)
                                    gCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCttt                                                                                                >  1:555537/1‑143 (MQ=255)
                                        aaTCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGGCTTTTTTTGTGCCTGAAATCAAGTGATTATGTATCTTTCCCg                                                                                            >  2:463572/1‑143 (MQ=255)
                                           cGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCta                                                                                         >  2:49336/1‑143 (MQ=255)
                                              tAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTatat                                                                                      >  1:603142/1‑143 (MQ=255)
                                               gggcccAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATAtt                                                                                     <  2:870196/138‑1 (MQ=255)
                                                 tccctgTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTaa                                                                                   >  2:415526/6‑143 (MQ=255)
                                                 acccAGGTGGTTCAAAATCCACCCTAGAAACACGGGCCGGTTCGAGTCCGGCCCTCCGCACCCAACGTATGTAAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTa                                                                                    <  1:17854/140‑1 (MQ=255)
                                                   gCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAAct                                                                                 >  2:474306/1‑143 (MQ=255)
                                                   gCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAAct                                                                                 >  2:883303/1‑143 (MQ=255)
                                                      gTTGAGTCAACATCAGCCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAActctc                                                                              <  2:907653/143‑1 (MQ=255)
                                                          aTTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCaa                                                                          <  2:661984/143‑1 (MQ=255)
                                                           ttCAAAATCATCCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAg                                                                         <  2:112843/143‑1 (MQ=255)
                                                            tCAAAATCAACCGTAGAAATACGTGCCGGTTCGGGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATTCCTGAGATCCAGTGGTTAGGTATCTTTCACGCTGTATTAAGTCTCGCaatt                                                                        >  1:876781/1‑141 (MQ=255)
                                                              aaaaTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCCCTCTATTAACTCTCTCAAGGTc                                                                      >  2:584341/1‑143 (MQ=255)
                                                                aaTCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCaa                                                                    >  2:735485/1‑143 (MQ=255)
                                                                  tCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAAcc                                                                  >  1:581108/1‑143 (MQ=255)
                                                                  tCAACCGCATAAATGCGTGCCGGTTCGAGGCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAAcc                                                                  <  2:528658/143‑1 (MQ=255)
                                                                       cGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTAAAGGTCAACCGatat                                                             >  1:749855/1‑143 (MQ=255)
                                                                         tAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCa                                                           <  2:630907/143‑1 (MQ=255)
                                                                         tAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCAAGTGTTTATATATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGataaac                                                           >  2:632801/1‑140 (MQ=255)
                                                                             aaTACGTGCCGTTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGt                                                       <  2:667504/143‑1 (MQ=255)
                                                                              aTACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTa                                                      <  1:735485/143‑1 (MQ=255)
                                                                                aCGTGCCGGTTCGAGTACGGCCTTCGGCACCAAAAGTATGTAAATAGACTTCAACTGAGGTCTTTTTTTATGGCTGAAATCAAGTGTTTATCTATCTTTTACGCTATATTAACTCTCTCACGGTCAAACGATATCAACGTACa                                                    >  2:770525/1‑143 (MQ=255)
                                                                                   tGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCt                                                 >  2:375593/1‑143 (MQ=255)
                                                                                        gTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTCTGCCTGGAAACCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTATCAAAGAAAACCGATAACCAAGTCCATCTACCaa                                            >  2:848021/1‑143 (MQ=255)
                                                                                        gTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGCCGACCGATATCAACGTACATCTACCaa                                            >  2:547094/1‑143 (MQ=255)
                                                                                                 cgggcgtCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTg                                   <  2:6818/137‑1 (MQ=255)
                                                                                                 cGGCATTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAGATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACACCTACCAACACATCttc                                   >  2:237535/1‑142 (MQ=255)
                                                                                                    ccgtCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTa                                <  1:707142/140‑1 (MQ=255)
                                                                                                    ccTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTCACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCCACATCTGTTgtga                                >  2:732601/1‑140 (MQ=255)
                                                                                                     cTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTAc                               >  2:677287/1‑143 (MQ=255)
                                                                                                        cGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAga                            <  1:435323/143‑1 (MQ=255)
                                                                                                         ggCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAgat                           <  2:729076/143‑1 (MQ=255)
                                                                                                          gCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAgatg                          >  1:451167/1‑143 (MQ=255)
                                                                                                          gCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGCTGGTACAgatg                          >  2:724447/1‑143 (MQ=255)
                                                                                                          gCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCCTCCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCACAATATGTTGTTCCAgatg                          >  2:553776/1‑143 (MQ=255)
                                                                                                             ccAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGatgatg                       <  1:490234/143‑1 (MQ=255)
                                                                                                             ccAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGatgatg                       <  1:375593/143‑1 (MQ=255)
                                                                                                                aaaGTAGGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACCCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTa                    <  1:845306/143‑1 (MQ=255)
                                                                                                                   ttATGTAAATAGACCTAGACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATAAACTCTCTCTAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTAttt                 <  2:126497/142‑1 (MQ=255)
                                                                                                                          aaaaAGACCTCAGTTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGtt           >  1:148361/1‑143 (MQ=255)
                                                                                                                           actAGACCTCAACTGAGGTCTTTTGTTATGGCTGAAATCCAGTGTTTATCTAGCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTc          <  2:207810/141‑1 (MQ=255)
                                                                                                                           aaaaGACCTCAGTTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTc          >  2:388442/4‑143 (MQ=255)
                                                                                                                           aaTAGACCTCAACTGAGGTCTTTTTTCATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTc          <  1:113094/143‑1 (MQ=255)
                                                                                                                              aGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGAt       >  1:216407/1‑143 (MQ=255)
                                                                                                                                  cTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCCACGTACATCTACCAACATATGTTGGATCAGATGATGGTATTTCCGGTTCGATAATg   >  2:654874/1‑143 (MQ=255)
                                                                                                                                   tCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGCTCGATAATGc  >  2:461982/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                          
ACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCAACATATGTTGGTACAGATGATGGTATTTCCGGTTCGATAATGC  >  NC_000913/4496378‑4496650

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: