Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	448,547	A→T	6.0%	I376F (ATC→TTC)	PP_0369 →	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	448,547	0	A	T	6.0%	203.6 / 2.5	67	I376F (ATC→TTC)	PP_0369	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base T (2/2);  total (34/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 = 8.58e-01

GCGGAATGCGTGCCGCATGTGGAAGCCGTGCGTGAAGTGATCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATGGTGCATGGGCGACAATCACC  >  NC_002947/448400‑448680                                                                                                                                                    |                                                                                                                                      gcgGAATGCGTGCCGCATGTGGAAGCCGTGCGTGAAGTGATCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCAt                                                                                                                                      <  1:363988/149‑1 (MQ=255)   ggAATGCGTGCCGCATGTGGAAGCCGTGCGTGAAGTGATCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCa                                                                                                                                    <  5:12339/149‑1 (MQ=255)       tGCGTGCCGCATGTGGAAGCCGTGCGTGAAGTGATCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGTGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTAt                                                                                                                                <  7:233310/149‑1 (MQ=255)               gCATGTGGAAGCCGTGCGTGAAGTGATCGCCAATTACGTGATGCACCTGCGCGCTCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTgg                                                                                                                        <  2:341045/149‑1 (MQ=255)                     ggAAGCCGTGCGTGAAGTGATCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAgc                                                                                                                  >  6:54242/1‑149 (MQ=255)                      gAAGCCGTGCGTGAAGTGATCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAgcg                                                                                                                 >  8:178691/1‑149 (MQ=255)                               cgtgAAGTGATCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGa                                                                                                        <  7:294992/149‑1 (MQ=255)                                  gAAGTGATCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGatca                                                                                                     >  7:273753/1‑149 (MQ=255)                                  gAAGTGATCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGatca                                                                                                     >  6:11784/1‑149 (MQ=255)                                      tGATCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGc                                                                                                 <  5:164260/149‑1 (MQ=255)                                        aTCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCaa                                                                                               >  6:235142/1‑149 (MQ=255)                                            ccAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAAcccc                                                                                           >  8:158465/1‑149 (MQ=255)                                            ccAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAAcccc                                                                                           >  3:353727/1‑149 (MQ=255)                                              aaTTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAAcccc                                                                                           >  2:235764/1‑147 (MQ=255)                                               aTTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAg                                                                                        <  7:158465/149‑1 (MQ=255)                                                  aCGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGACGATCATCGCAACCCCGAGGCg                                                                                     >  2:266524/1‑149 (MQ=255)                                                           aCCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGatcatc                                                                                                   <  1:296463/126‑1 (MQ=255)                                                           aCCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGatcatc                                                                                                   >  2:296463/1‑126 (MQ=255)                                                              tGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATcc                                                                         >  6:389523/1‑149 (MQ=255)                                                               gcgcGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAAcc                                                                                             <  4:339878/128‑1 (MQ=255)                                                               gcgcGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAAcc                                                                                             >  3:339878/1‑128 (MQ=255)                                                                   gATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTg                                                                               <  2:218312/138‑1 (MQ=255)                                                                   gATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTg                                                                               >  1:218312/1‑138 (MQ=255)                                                                     tCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGAc                                                                   >  8:85932/1‑148 (MQ=255)                                                                       aGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAg                                                                <  4:216256/149‑1 (MQ=255)                                                                        gCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAgg                                                               <  5:54242/149‑1 (MQ=255)                                                                         cacaACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGgc                                                              >  8:191742/1‑149 (MQ=255)                                                                          acaACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGgcg                                                             >  3:350093/1‑149 (MQ=255)                                                                           caACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGgcgc                                                            >  6:165730/1‑149 (MQ=255)                                                                            aaCCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGc                                                                                 <  1:55244/127‑1 (MQ=255)                                                                            aaCCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGc                                                                                 >  2:55244/1‑127 (MQ=255)                                                                            aaCCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGa                                                                    >  1:21015/1‑140 (MQ=255)                                                                            aaCCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGa                                                                    <  2:21015/140‑1 (MQ=255)                                                                            aaCCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCt                                                           >  1:58860/1‑149 (MQ=255)                                                                            aaCCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCt                                                           >  2:224543/1‑149 (MQ=255)                                                                                 cccgcAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTa                                                        >  6:146148/1‑147 (MQ=255)                                                                                   cgcAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGt                                                    <  7:178691/149‑1 (MQ=255)                                                                                   cgcAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGt                                                    <  4:353727/149‑1 (MQ=255)                                                                                    gcAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGTCGATCATCGCTACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTg                                                   <  1:266524/149‑1 (MQ=255)                                                                                       ggACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGAGcc                                                 >  3:194615/1‑148 (MQ=255)                                                                                        gACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGACACGGTCTCGGTCACCTTCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCa                                                >  8:113269/1‑148 (MQ=255)                                                                                        gACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGACACGGTCTCGGTCACCTTCAGTATCGGCGTGGCCGAGCGCAAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCa                                                >  8:107949/1‑148 (MQ=255)                                                                                             tACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCTGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGgcg                                          >  8:261955/1‑149 (MQ=255)                                                                                             tACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGgc                                           <  7:322054/148‑1 (MQ=255)                                                                                             tACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGgc                                           >  8:322054/1‑148 (MQ=255)                                                                                                   cgGGCGCCAGCGCCGGGCTGGCAGCGCCTGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGgc                                    <  7:261955/149‑1 (MQ=255)                                                                                                         ccAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAAtt                              <  5:389523/149‑1 (MQ=255)                                                                                                            gcgcCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATtgtg                           <  7:89363/149‑1 (MQ=255)                                                                                                               ccGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACACTGCCAAAGGCGCGGGGCGCAATTGTGTCa                        <  1:235764/149‑1 (MQ=255)                                                                                                               ccGGGCTGGCAGCGCCGGCGACACGGTCTCGGTCACCTTCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCa                        <  7:113269/149‑1 (MQ=255)                                                                                                               ccGGGCTGGCAGCGCCGGCGACACGGTCTCGGTCACCTTCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCa                        <  7:107949/149‑1 (MQ=255)                                                                                                                  ggCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATgg                     <  2:58860/149‑1 (MQ=255)                                                                                                                   gCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATGGt                    <  1:85655/149‑1 (MQ=255)                                                                                                                   gCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATGGt                    <  8:273753/149‑1 (MQ=255)                                                                                                                   gCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATGGt                    <  5:235142/149‑1 (MQ=255)                                                                                                                    cTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATGGTg                   <  5:11784/149‑1 (MQ=255)                                                                                                                     tGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATGGTg                   >  2:222072/1‑148 (MQ=255)                                                                                                                      ggCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATGGTGCa                 >  7:1988/1‑149 (MQ=255)                                                                                                                       gCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATGGTGCAt                >  5:372650/1‑149 (MQ=255)                                                                                                                             ccGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATGGTGCATGGGCGa          <  1:224543/149‑1 (MQ=255)                                                                                                                              cggcggCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATtgt                            <  1:276953/130‑1 (MQ=255)                                                                                                                              cggcggCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATtgt                            >  2:276953/1‑130 (MQ=255)                                                                                                                              cggcggCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATgg                     <  8:153267/137‑1 (MQ=255)                                                                                                                              cggcggCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATgg                     >  7:153267/1‑137 (MQ=255)                                                                                                                                     aCGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATGGTGCATGGGCGACAATCAcc  >  3:269216/1‑149 (MQ=255)                                                                                                                                       ggTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCaa                                >  6:66052/1‑117 (MQ=255)                                                                                                                                       ggTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCaa                                <  5:66052/117‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                      GCGGAATGCGTGCCGCATGTGGAAGCCGTGCGTGAAGTGATCGCCAATTACGTGATGCACCTGCGCGATCAGCACAACCGCCCGCAGGACGATACCGCCGGGCGCCAGCGCCGGGCTGGCAGCGCCGGCGGCACGGTCTCGGTCACCATCAGTATCGGCGTGGCCGAGCGCCAGGCCGATCATCGCAACCCCGAGGCGGTGCTGAAATCCGCCGACCAGGCGCTTTACAGTGCCAAAGGCGCGGGGCGCAATTGTGTCATGGTGCATGGGCGACAATCACC  >  NC_002947/448400‑448680

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