Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,177,171 T→A 10.2% intergenic (+43/‑187) tufB → / → secE translation elongation factor Tu 2/Sec translocon subunit SecE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,177,1710TA10.2% 169.3 / 10.3 68intergenic (+43/‑187)tufB/secEtranslation elongation factor Tu 2/Sec translocon subunit SecE
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base A (5/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

GATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGA  >  NC_000913/4177032‑4177312
                                                                                                                                           |                                                                                                                                             
gATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTACAAGGGCATCATTTg                                                                                                                                            >  2:555238/1‑143 (MQ=255)
gATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTg                                                                                                                                            >  2:547587/1‑143 (MQ=255)
 aTCCCCCCGATCGCGATGGACAACGGAATGCATTTCGTAAGTCTTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGa                                                                                                                                           <  1:167165/143‑1 (MQ=25)
 aTCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACGTTTGACGCAATGCGCACTAAAAGGGCATCGCTTGa                                                                                                                                           >  2:136312/1‑143 (MQ=255)
    cACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTGCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTCAAAGGGCAGCAACTGATGc                                                                                                                                        >  2:200688/1‑143 (MQ=21)
       gcGATCGGGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCt                                                                                                                                     <  2:569270/142‑1 (MQ=255)
        cGATCGGGAGGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtt                                                                                                                                    <  1:109552/143‑1 (MQ=255)
           tCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttta                                                                                                                                 <  2:157007/143‑2 (MQ=255)
            cgcgATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTg                                                                                                                                <  2:106934/143‑1 (MQ=255)
               gATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGTCGGGCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGCGCTACTTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCAc                                                                                                                             <  1:10984/143‑1 (MQ=255)
                   gacgacGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCtt                                                                                                                         <  1:75415/143‑1 (MQ=255)
                         ggTCGGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAc                                                                                                                   <  1:488791/143‑1 (MQ=255)
                               cGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCGGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATCTGATACCCGTGTTGCACGCCTTCGCACCAGAAc                                                                                                             >  2:210302/1‑143 (MQ=37)
                               cGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAACAAAGGTTCTGAGCTAATTGGCGATGGCATTTGACGCAATGCGCACGAAAAGGGCATCATTTGATGCCCTTTTGGCACGCTTGCGTAACAGaag                                                                                                             >  2:347451/1‑142 (MQ=25)
                                 tttCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCt                                                                                                           >  1:725032/1‑143 (MQ=255)
                                   tCGGAATCCGTGGAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTgg                                                                                                         <  2:53584/143‑1 (MQ=255)
                                    cGCAAGCCGGGAGGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATACCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGc                                                                                                        <  1:891405/143‑1 (MQ=255)
                                     gCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTATGAGCTAATTACCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAACGGAACCTGGCt                                                                                                       >  2:555524/1‑143 (MQ=255)
                                     gCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGAGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGCTGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCt                                                                                                       >  1:743020/1‑143 (MQ=255)
                                      cAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCCTTCGTACCAGAACCTGGCtc                                                                                                      >  2:235513/1‑143 (MQ=255)
                                       aaTCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGCGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtca                                                                                                     >  2:474274/1‑143 (MQ=255)
                                       aaTCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTAGACGCAATGCGCACTAAAAGGGCATCATTCGATGCCCTTTTTGCCCGCTCTCGTACCCGAACCTGGcacc                                                                                                     >  2:761877/1‑140 (MQ=37)
                                       aaTCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATGGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtca                                                                                                     >  1:897528/1‑143 (MQ=255)
                                         tCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCccatc                                                                                                   >  1:416925/1‑143 (MQ=255)
                                          ccGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatca                                                                                                  <  2:378997/143‑1 (MQ=255)
                                          ccGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatca                                                                                                  <  2:668616/143‑1 (MQ=255)
                                             tGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCATTg                                                                                               >  1:556919/1‑143 (MQ=255)
                                                  gcggcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATACCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAtttt                                                                                          >  2:615559/1‑143 (MQ=255)
                                                  gcggcCGTACCGGTGGCGCGGGCGTTGTAGCAACAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGAtttt                                                                                          <  2:519616/143‑1 (MQ=255)
                                                    ggcGGTACCGTGGGCGCGGGCGTTGGAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCt                                                                                        <  2:895665/143‑1 (MQ=255)
                                                     gcCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCtt                                                                                       <  1:53876/143‑1 (MQ=255)
                                                      ccGTGCCGTTGGCGCGGGCGTTGTAGCAAAGGTTCTGAGCTGATTGCCGATACTATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCttt                                                                                      <  2:643967/143‑1 (MQ=255)
                                                         tACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTc                                                                                   >  1:693097/1‑143 (MQ=255)
                                                                  cgcgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAtt                                                                          <  1:313085/143‑1 (MQ=255)
                                                                     gggCGTTGCAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCCtttcg                                                                       <  2:855840/143‑4 (MQ=255)
                                                                         tttttAGCAAAAGTTCTGTGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTgaga                                                                   <  1:687035/139‑1 (MQ=255)
                                                                            gTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAGTCATTGCTGAGACAg                                                                >  1:568818/1‑143 (MQ=255)
                                                                                cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCCGTGATTTTCTTTGTCATAATCATTGCCGAGACCGGccc                                                            >  1:852709/1‑141 (MQ=255)
                                                                                cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTTTCATAATCATTGATGAGACAggacc                                                            >  1:303584/1‑140 (MQ=255)
                                                                                cAAAAGTTCGGAGCGAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGTTCATCATTTGATGCCCTGTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGctc                                                            <  1:171433/143‑1 (MQ=255)
                                                                                cAAAAGGTCTGAGCTAATTGCCGATAACATGTGACGCAATGCGCACTAAAAGTGCATCATTTGATGCACTTTTTGCACGCTTTCGTACGAGAACCTGGATAAGCAGTGATTTTCTCTGTCATAAGCAGTGCTGAGGCAggagg                                                            >  2:813534/1‑140 (MQ=255)
                                                                                  aaaGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGGTGAGACCGGCTCTg                                                          >  1:639106/1‑143 (MQ=255)
                                                                                    aGTTCGGAGCTAATTGCAGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGtt                                                        <  2:839138/143‑1 (MQ=255)
                                                                                      ttCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa                                                      >  1:281424/1‑143 (MQ=255)
                                                                                       tCTGAGCTAATTGCCGATAACATTGGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAg                                                     <  1:863388/143‑1 (MQ=255)
                                                                                         tGAGCTAATTGCCGATAACATTTGACGCAATGCGCCCTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATCTTCTTTGTCATACTCATTGCTGAGACCGGCGCTGATGAggg                                                   >  2:858075/1‑143 (MQ=255)
                                                                                          gAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGc                                                  <  2:750269/143‑1 (MQ=255)
                                                                                                  tGCCGATAACATTTGACGCAATGCGCACTACAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGCCAGGCTCTGTTGAGGGCGTATAATc                                          <  1:438952/143‑1 (MQ=255)
                                                                                                       aTAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATAAGTGATCTTCTTTGTCATAATCATTGCTGACACAGGCTCTGTTGGGGGCGTATCACCCGacg                                     >  2:661077/1‑141 (MQ=255)
                                                                                                        tAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTGTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGACCGTTTGAGGGCGTATAATCCGacaa                                    >  2:222198/1‑143 (MQ=255)
                                                                                                             tttGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTaa                               <  1:547587/143‑1 (MQ=255)
                                                                                                             tttGACGCAATCCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTaa                               <  1:555238/143‑1 (MQ=255)
                                                                                                             tcggccGCTATGGGCACTACAAGGGCATCATTTGATGCCCGTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTaa                               <  2:588987/138‑1 (MQ=255)
                                                                                                                    cAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGt                        >  1:450700/1‑143 (MQ=255)
                                                                                                                     aaTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGGGTGTAATCCGAAAAGCTAATACGCGtt                       >  1:863086/1‑143 (MQ=255)
                                                                                                                     aaTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGtt                       <  1:801817/143‑1 (MQ=255)
                                                                                                                       cgcgcACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTc                     <  1:353373/142‑1 (MQ=255)
                                                                                                                           tgcaaaaaGGGCATCAAATGATGCCCTTTTAGTGCGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGGTCCGAtt                 >  2:382078/5‑143 (MQ=38)
                                                                                                                           cACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGAtt                 <  1:879926/143‑1 (MQ=255)
                                                                                                                              aaaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTgg              >  1:392253/2‑143 (MQ=255)
                                                                                                                                aaaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCAGAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGtt            <  1:439189/143‑1 (MQ=255)
                                                                                                                                aaaGGGCACCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCAGAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGtt            <  1:231549/143‑1 (MQ=255)
                                                                                                                                 aaGGGCATGATTTGATGCCCTTTCTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGttt           <  1:18845/143‑1 (MQ=255)
                                                                                                                                 aaGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCAGTGCTGAGGCACGCGCTGTTGAGGGCGTCAAAGCACACAAGAAAGTACTCGTGTTGATTTGGttt           >  2:359122/1‑143 (MQ=255)
                                                                                                                                  aGGGCATCAAATGATGCCCTTTTAGTGCGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGGTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTg          >  1:460747/1‑143 (MQ=255)
                                                                                                                                  aGGGCATCAAATGATGCCCTTTTAGTGCGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGCGGGCGTATAATCCGAAAAGCTAATACGCGGTTCGATTTGGTTTg          >  1:338542/1‑143 (MQ=255)
                                                                                                                                    ggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGCGACAGGCTCCGTTGAGGGCGTATAATCCGAAACGCTAATCCGCGTCTCGATTTGGTTTGcc        >  1:419855/1‑143 (MQ=255)
                                                                                                                                     gCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCGTCAGTGGTTTTCTTTGTCATAATCATTGCGGAGCCCGGCGCTGTTGAGGGCGTATAATCCGCAAAGATAATACGCGTTTAGCATTGGTGTGGCt       >  2:741007/1‑143 (MQ=255)
                                                                                                                                        tcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACCGGCTCTGTTGCGGGCGTATACTCCGAAAAGCTAGTACGAGTGTCGATTTGGTTTGCCTcga    >  1:165872/1‑142 (MQ=255)
                                                                                                                                          atTTGATGCCCTGTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGa  <  2:555739/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
GATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCGA  >  NC_000913/4177032‑4177312

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: