Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,407,076	A→T	9.6%	T240T (ACT→ACA)	dgcM ←	diguanylate cyclase DgcM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,407,076	0	A	T	9.6%	155.6 / 9.5	63	T240T (ACT→ACA)	dgcM	diguanylate cyclase DgcM
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (3/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-01

TAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGA  >  NC_000913/1406942‑1407218                                                                                                                                       |                                                                                                                                                gaaTCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATc                                                                                                                                         >  2:329631/2‑143 (MQ=255)    tCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATg                                                                                                                                      <  2:24644/143‑1 (MQ=255)    tCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATg                                                                                                                                      >  1:713622/1‑143 (MQ=255)    tCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTCCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATg                                                                                                                                      >  1:694872/1‑143 (MQ=255)    tCCGGAGCGATGGCGATATGCTGCATTGGGCCTGGTCCCGTAATGTGATAAAACTGTCGCCGATTCAGCAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATg                                                                                                                                      <  1:713511/143‑1 (MQ=255)       tGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATGCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCGGCCGTTCTTGCCCAGCAATGCCATGcac                                                                                                                                   >  2:436960/1‑143 (MQ=255)        gAGCGATGGCGAGATGCTGCATGTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGcaca                                                                                                                                  <  2:495220/143‑1 (MQ=255)         aGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACaa                                                                                                                                 <  1:127924/143‑1 (MQ=255)           cggtGGCGGGATGCTGCATTTGGTCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATa                                                                                                                               <  2:365250/140‑1 (MQ=255)               ggCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATa                                                                                                                           >  1:270442/1‑143 (MQ=255)                  gagaTGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGCCTCAGTAATCCGCTCATCGCGTCCTGGTGAGCAGCCTGTTCCCGCTGCTCCTCCAGCCGTTCTTGCACAGTAATATAAAGCACAATACATaaaa                                                                                                                        >  1:932100/1‑141 (MQ=255)                  gagaTGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACGATACATAACa                                                                                                                        >  1:454560/1‑143 (MQ=255)                       gctgcATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTGATATCATGCACAATACATAACAT‑GAGc                                                                                                                  >  1:132956/1‑143 (MQ=255)                         tgcATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCACGTTCCAGCTGCTCCTCCAGCCGTTTTTGCGCAGTAATATCATGCACAATACATAACAT‑GCGCtt                                                                                                                >  2:312856/1‑143 (MQ=255)                           cATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGt                                                                                                              >  1:39032/1‑143 (MQ=255)                              ttGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTGAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGc                                                                                                           >  1:832103/1‑143 (MQ=255)                              ttGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCTGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGc                                                                                                           >  2:524820/1‑143 (MQ=255)                              ttGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCTGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGc                                                                                                           >  2:519467/1‑143 (MQ=255)                              ttGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCTGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGc                                                                                                           >  2:24682/1‑143 (MQ=255)                                  ccTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATa                                                                                                       >  1:258446/1‑143 (MQ=255)                                    tGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAAGCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCGCCAGCCGTTTTTGCGCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATaaa                                                                                                     <  1:64746/143‑1 (MQ=255)                                          cGTAATGTGATAAAACTGTCGCCGCTCCAGTACCCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGCCATATCACCCACACTACATTTTAT‑GAGCTTGTCGCCATAAATTTCaa                                                                                               <  1:500901/143‑1 (MQ=255)                                           gTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAAt                                                                                              <  2:283326/143‑1 (MQ=255)                                                    tAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGc                                                                                     >  1:447541/1‑143 (MQ=255)                                                        aCTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTCTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAg                                                                                 <  1:548158/143‑1 (MQ=255)                                                            tCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCt                                                                             >  1:578580/1‑143 (MQ=255)                                                              gCCGATTCGGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGc                                                                           >  1:112076/1‑143 (MQ=255)                                                               ccGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCCCAATACATAACAT‑GAGCTCGTACCCACAAACTCCAATCGGTCCGCGCTAGGACTGca                                                                          >  1:1037977/1‑143 (MQ=255)                                                                cGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGcac                                                                         >  1:226028/1‑143 (MQ=255)                                                                  aTTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACAt                                                                       <  1:329631/143‑1 (MQ=255)                                                                     gcgTACTCCGGTCATGGCGTCATGGGGGGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATTCATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGAc                                                                    <  2:932100/141‑1 (MQ=255)                                                                      agcactCCGGTCATCGCGTCATGGTGTGCACCATGTTCCATCTGCTCCTCCAGCCGTTTTTGCTCAGTAACATCATGCACAATACACACCAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGCCATAGGTCTGCACATGACg                                                                   <  2:554753/138‑1 (MQ=255)                                                                      aGTAATCCGGTCATCGCGTCATGGTTAGCAGCATGTTCTAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATACCAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACg                                                                   >  1:477719/1‑143 (MQ=255)                                                                        tAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAg                                                                 <  2:39032/143‑1 (MQ=255)                                                                         aaTCCGGTCCTCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGt                                                                <  1:343416/143‑1 (MQ=255)                                                                             cGGTCGTCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATGACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGaa                                                            <  2:1069513/143‑1 (MQ=255)                                                                             cGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCCTAGGTCTGCACATGACGCGTCGaa                                                            >  2:461608/1‑143 (MQ=255)                                                                             cGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGaa                                                            <  1:929200/143‑1 (MQ=255)                                                                                tCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGCGTCGAACCa                                                         >  1:429379/1‑143 (MQ=255)                                                                                  aTCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACACTACATAACAT‑GAGCTTGTGGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACCTGACGAGCCGAACAATc                                                       >  1:175293/1‑143 (MQ=255)                                                                                    cgcgTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATACATATCCATCGGTCCGGCGTACGTCTGCACCGGGCGAGTCGAACCATACg                                                     >  1:954379/1‑143 (MQ=255)                                                                                    cgcgTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCg                                                     >  1:76743/1‑143 (MQ=255)                                                                                        tcttGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATCTCATGCACAATACATACCAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAg                                                 <  2:1004774/140‑1 (MQ=255)                                                                                           tGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGttt                                              >  2:880951/1‑143 (MQ=255)                                                                                               gAGCAGCATGTTCCGGCTGCTCCTCCAGCCGTTCTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGa                                          <  2:947586/143‑1 (MQ=255)                                                                                                    gcaTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACaaa                                     <  2:112076/143‑1 (MQ=255)                                                                                                      aTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATCCCGAGTCGCACCACCCGCCAGTTTATGAACacaaa                                   >  2:112943/1‑142 (MQ=255)                                                                                                       tGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAAtt                                  <  2:578580/143‑1 (MQ=255)                                                                                                         ttCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCa                                <  2:226028/143‑1 (MQ=255)                                                                                                           ccAGCTGCTCCTCCAGCCGTTTTTGCTCTGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaaa                              <  1:524820/143‑1 (MQ=255)                                                                                                           ccAGCTGCTCCTCCAGCCGTTTTTGCTCTGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaaa                              <  1:24682/143‑1 (MQ=255)                                                                                                           ccAGCTGCTCCTCCAGCCGTTTTTGCTCTGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaaa                              <  1:519467/143‑1 (MQ=255)                                                                                                           ccAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaaa                              <  1:599349/143‑1 (MQ=255)                                                                                                            cAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAACTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATAACGAGTCGAACCATCCACCAGTATACGAACAAAATTCaaac                             >  2:552405/1‑142 (MQ=255)                                                                                                                  ctcctcCAGCCGTTTTTGCTCAGTAATTTCATGCAGAATACATAACAT‑GAGCTTGTCGCCAGAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTAt                       <  2:171048/143‑1 (MQ=255)                                                                                                                    cctcctGCCGCTTTTACTCAGTACTACCACGCACCACAC‑TTACATCGATCTGTTCGCCATAAATTTCAATCGGCCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGa                     <  1:1015052/143‑1 (MQ=255)                                                                                                                     ctccAGCCGTTTTTGCTCAGTCAGATCATGCACACTACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGac                    <  2:833765/143‑1 (MQ=255)                                                                                                                       ccAGCCGTTGTTGCTCAGTAATATCAGGCACAATACATAACAT‑GAGCTTGTCGCCATAACTTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGacca                  <  2:122133/143‑1 (MQ=255)                                                                                                                       ccAGCCGGTTTGGCACAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGacca                  <  1:1006445/143‑1 (MQ=255)                                                                                                                           ccGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGc                                                    <  2:23531/105‑1 (MQ=255)                                                                                                                           ccGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGc                                                    >  1:23531/1‑105 (MQ=255)                                                                                                                                    cTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCg     <  1:625024/143‑1 (MQ=255)                                                                                                                                     tCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCga    <  2:447541/143‑1 (MQ=255)                                                                                                                                       aGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTACGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCCGTTTATGAACCACATTCAAAGGCTTATGACCACCGGGTAAATGCgagg  >  2:660459/1‑142 (MQ=255)                                                                                                                                       |                                                                                                                                                TAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAT‑GAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGA  >  NC_000913/1406942‑1407218

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