Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,513,032	A→G	6.7%	E567G (GAA→GGA)	PP_3108 →	rhs family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,513,032	0	A	G	6.7%	172.2 / 2.1	60	E567G (GAA→GGA)	PP_3108	rhs family protein
Reads supporting (aligned to +/- strand):  ref base A (30/25);  new base G (2/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGGGGAGCGTTACGCCATCGAGCTGGACGAACACGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGC  >  NC_002947/3512885‑3513177                                                                                                                                                    |                                                                                                                                                      gCGGGGAGCGTTACGCCATCGAGCTGGACGAACACGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGaa                                                                                                                                                      <  2:361104/149‑1 (MQ=255)   ggggAGCGTTACGCCATCGAGCTGGACGAACACGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCGCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGCata                                                                                                                                                    >  8:39946/1‑149 (MQ=255)     ggAGCGTTACGCCATCGAGCTGGACGAACACGGCAATAAGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAatat                                                                                                                                                   >  4:186841/1‑148 (MQ=255)       aGCGTTACGCCATCGAGCTGGACGAACACGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATcac                                                                                                                                                <  1:2072/149‑1 (MQ=255)          gTTACGCCATCGAGCTGGACGAACACGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATcaccac                                                                                                                                             <  2:190604/149‑1 (MQ=255)                   tCGAGCTGGACGAACACGGCAATATGGTTGGCCTAACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGGATATCACCACCTGACcaca                                                                                                                                    >  5:138512/1‑149 (MQ=255)                   tCGAGCTGGACGAACACGGCAATATGGTTGGCCTAACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGGATATCACCACCTGACcaca                                                                                                                                    >  5:132695/1‑149 (MQ=255)                          ggACGAACACGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGa                                                                                                                           >  4:226744/1‑149 (MQ=255)                              gAACACGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCa                                                                                                                       <  1:158824/149‑1 (MQ=255)                               aaCACGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAg                                                                                                                      <  4:108408/149‑1 (MQ=255)                                acacGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAatat                                                                                                                                                   >  8:69990/1‑121 (MQ=255)                                acacGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAatat                                                                                                                                                   <  7:69990/121‑1 (MQ=255)                                  acGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAat                                                                                                                                                     >  1:28032/1‑117 (MQ=255)                                  acGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAat                                                                                                                                                     <  2:28032/117‑1 (MQ=255)                                       aaTATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTa                                                                                                            >  2:231807/1‑149 (MQ=255)                                       aaTATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTa                                                                                                            >  8:32781/1‑149 (MQ=255)                                             gTTGGCCTGACCTTGCCGGACGGTAACTGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGaa                                                                                                                                                      >  7:221976/1‑105 (MQ=255)                                             gTTGGCCTGACCTTGCCGGACGGTAACTGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGaa                                                                                                                                                      <  8:221976/105‑1 (MQ=255)                                             gTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCgg                                                                                                       <  3:186841/148‑1 (MQ=255)                                             gTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGa                                                                                                      >  8:261829/1‑149 (MQ=255)                                                 gCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATgg                                                                                                   >  5:263157/1‑148 (MQ=255)                                                   cTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGgcag                                                                                                >  6:134488/1‑149 (MQ=255)                                                      aCCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGGATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGgcagca                                                                                              <  6:138512/148‑1 (MQ=255)                                                      aCCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGGATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGgcagca                                                                                              <  6:132695/148‑1 (MQ=255)                                                      aCCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCa                                                                                                                       <  1:190568/125‑1 (MQ=255)                                                      aCCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCa                                                                                                                       >  2:190568/1‑125 (MQ=255)                                                      aCCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCAc                                                                                             >  8:45162/1‑149 (MQ=255)                                                        cTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACtt                                                                                           >  2:248826/1‑149 (MQ=255)                                                         ttGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTg                                                                                          >  1:103370/1‑149 (MQ=255)                                                           gCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTg                                                                                          <  7:45162/147‑1 (MQ=255)                                                             cGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGa                                                                                      <  6:263157/149‑1 (MQ=255)                                                                  ggTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGcgcg                                                                                 >  6:59379/1‑149 (MQ=255)                                                                      aCCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTc                                                                             <  3:226744/149‑1 (MQ=255)                                                                       ccGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTc                                                                             >  8:175038/1‑148 (MQ=255)                                                                           cTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGt                                                                                                                    <  8:144397/107‑1 (MQ=255)                                                                           cTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACAACCTGACCACACT‑GG‑TGACCCAGGt                                                                                                                    >  7:144397/1‑107 (MQ=255)                                                                            tGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTcgacga                                                                        >  7:279988/1‑148 (MQ=255)                                                                                 ttCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAggg                                                                  >  7:327366/1‑149 (MQ=255)                                                                                     aGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAAc                                                              <  8:296988/149‑1 (MQ=255)                                                                                      gTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAgg                                                                   >  2:249503/1‑143 (MQ=255)                                                                                      gTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAgg                                                                   <  1:249503/143‑1 (MQ=255)                                                                                          gacgaATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCt                                                         >  1:243853/1‑149 (MQ=255)                                                                                          gacgaATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCt                                                         >  6:80204/1‑149 (MQ=255)                                                                                           acgaATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTc                                                        <  7:261829/149‑1 (MQ=255)                                                                                                  ctcCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACTAGGAGGACCCAGGTCAAGACTACACGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGcc                                                     >  1:299671/1‑149 (MQ=17)                                                                                                   tcCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGAGCCTTTCGACGACAAGGGCAACCTGCTCGCCGAAtt                                                >  7:167907/1‑149 (MQ=255)                                                                                                    cccGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGCCCCAGTACGAATATCACCACCCGACCACACT‑GG‑TGCCCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTc                                               >  7:113981/1‑149 (MQ=255)                                                                                                      cGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGa                                             >  5:334710/1‑149 (MQ=255)                                                                                                            cTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGccc                                        <  1:248826/148‑1 (MQ=255)                                                                                                            cTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGcc                                         <  5:134488/147‑1 (MQ=255)                                                                                                             tCGAAGAGATCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTc                                      >  8:297012/1‑149 (MQ=255)                                                                                                             tCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATg                                           <  6:181191/144‑1 (MQ=255)                                                                                                             tCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATg                                           >  5:181191/1‑144 (MQ=255)                                                                                                               gAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCgg                                    <  7:32781/149‑1 (MQ=255)                                                                                                                     accgaccCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGAt                              >  5:360430/1‑149 (MQ=255)                                                                                                                     accgaccCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGAt                              <  7:175038/149‑1 (MQ=255)                                                                                                                                 ggCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTAcc                   <  6:334710/148‑1 (MQ=255)                                                                                                                                 ggCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTAc                    <  5:80204/147‑1 (MQ=255)                                                                                                                                 ggCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCt                  <  1:231807/149‑1 (MQ=255)                                                                                                                                    cGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGaa               >  1:79025/1‑149 (MQ=255)                                                                                                                                      cAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTAc                    <  1:265234/142‑1 (MQ=255)                                                                                                                                      cAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTAc                    >  2:265234/1‑142 (MQ=255)                                                                                                                                                  cGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGc  <  2:103370/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GCGGGGAGCGTTACGCCATCGAGCTGGACGAACACGGCAATATGGTTGGCCTGACCTTGCCGGACGGTAACCGCCTGCAATTCAAGTACGACGAATTCTCCCGCCTGCTCGAAGAGACCGACCCGCTGGGCCGCAAGACCCAGTACGAATATCACCACCTGACCACACT‑GG‑TGACCCAGGTC‑AG‑CTACCCGGATGGCAGCACTTGGCGAGCGCGTTTCGACGACAAGGGCAACCTGCTCGCCGAATTCGATGCCCTCGGCCAGATGACCGAGTACCTGAACAGCGATGACGGC  >  NC_002947/3512885‑3513177

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: ‑