Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,405,314	C→T	10.0%	G83R (GGG→AGG)	PP_3886 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,405,314	0	C	T	10.0%	157.5 / 7.8	60	G83R (GGG→AGG)	PP_3886	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/23);  new base T (4/2);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

CGATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGC  >  NC_002947/4405173‑4405459                                                                                                                                                |                                                                                                                                                  cGATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCAtg                                                                                                                                            >  2:79803/1‑149 (MQ=255)  gATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCAtgt                                                                                                                                           >  7:217680/1‑149 (MQ=255)        cAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCggg                                                                                                                                     >  4:104177/1‑147 (MQ=255)         aGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAgg                                                                                                                                    >  1:246134/1‑149 (MQ=255)           cGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCC‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGc                                                                                                                                   <  1:110794/149‑1 (MQ=255)           cGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑CTCACGCTATCCATGACTGTCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGc                                                                                                                                   <  1:105491/149‑1 (MQ=255)            gAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCtt                                                                                                                                 <  7:320395/149‑1 (MQ=255)             aaCTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTg                                                                                                                                >  5:40050/1‑149 (MQ=255)              aCTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGGGCtttt                                                                                                                               <  3:104177/149‑3 (MQ=255)                tGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCTGCATGTGCCGAGGCTTGGaa                                                                                                                             >  5:137274/1‑149 (MQ=255)                  ccTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGGGCtttttc                                                                                                                             <  5:221094/147‑5 (MQ=255)                   cTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATgg                                                                                                                          >  6:144731/1‑149 (MQ=255)                        aaCAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGga                                                                                                                     <  4:286983/149‑1 (MQ=255)                           aGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GCCGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTgg                                                                                                                               >  5:13271/1‑136 (MQ=255)                           aGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTgg                                                                                                                               <  6:13271/136‑1 (MQ=255)                                  cccGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAAcc                                                                                                           <  3:163963/149‑1 (MQ=255)                                     gAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTg                                                                                                        >  2:401252/1‑149 (MQ=255)                                     gAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTg                                                                                                        >  6:417818/1‑149 (MQ=255)                                            gccgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTg                                                                                                                       <  2:140702/127‑1 (MQ=255)                                            gccgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTg                                                                                                                       >  1:140702/1‑127 (MQ=255)                                                      cGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTgg                                                                                       <  1:79803/149‑1 (MQ=255)                                                       gTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGgc                                                                                      >  2:346404/1‑149 (MQ=255)                                                            ggAATTTGGAGGGACTGGACAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTc                                                                                 >  7:50161/1‑149 (MQ=255)                                                               atgtGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCt                                                                              <  7:226290/146‑1 (MQ=255)                                                                   ggAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCa                                                                           <  6:40050/148‑1 (MQ=255)                                                                   ggAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAg                                                                          >  1:353344/1‑149 (MQ=255)                                                                   ggAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATCTGCCGAGGCTTGGATTGCGTGGAGATGGCAACCCTGACAGGTCGCTGTCGTGGCGCTTCCCTACAg                                                                          >  5:200732/1‑149 (MQ=255)                                                                         aCTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                             >  3:89469/1‑140 (MQ=255)                                                                         aCTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                             >  3:95050/1‑140 (MQ=255)                                                                         aCTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                             <  4:95050/140‑1 (MQ=255)                                                                         aCTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                             <  4:89469/140‑1 (MQ=255)                                                                          cTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCAc                                                                   >  3:95060/1‑149 (MQ=255)                                                                             ggCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCt                                                                <  8:217680/149‑1 (MQ=255)                                                                                aGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCCGGTCGCTGTCGTGGCGGTTCCCTCCAGTCCCCACGCTGAc                                                             >  7:299678/1‑149 (MQ=255)                                                                                aGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGAc                                                             <  1:346404/149‑1 (MQ=255)                                                                                 gAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACt                                                            <  2:246134/149‑1 (MQ=255)                                                                                      cTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTAc                                                       >  3:200597/1‑149 (MQ=255)                                                                                        gAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTAccc                                                     >  7:118730/1‑149 (MQ=255)                                                                                             tCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCa                                                >  1:275675/1‑149 (MQ=255)                                                                                                  cGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTg                                           <  5:144731/149‑1 (MQ=255)                                                                                                        c‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTgccgcc                                      >  8:280845/1‑149 (MQ=255)                                                                                                           aCGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCgg                                    >  7:88332/1‑149 (MQ=255)                                                                                                           aCGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCgg                                    <  1:401252/149‑1 (MQ=255)                                                                                                                 tCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCttt                               <  2:275675/148‑1 (MQ=255)                                                                                                                   cATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGa                             <  1:150151/148‑1 (MQ=255)                                                                                                                    aTGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGgc                                                                                      >  2:164587/1‑90 (MQ=255)                                                                                                                    aTGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGgc                                                                                      <  1:164587/90‑1 (MQ=255)                                                                                                                        cTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCTCCGGGCTTTGATCAtct                       <  6:357562/149‑1 (MQ=255)                                                                                                                                 tttAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGat                            <  7:304621/135‑1 (MQ=255)                                                                                                                                 tttAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGat                            >  8:304621/1‑135 (MQ=255)                                                                                                                                     aCGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACa           >  1:146288/1‑148 (MQ=255)                                                                                                                                       gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg        >  6:404080/1‑149 (MQ=255)                                                                                                                                       gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg        >  5:3883/1‑149 (MQ=255)                                                                                                                                       gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg        >  3:191196/1‑149 (MQ=255)                                                                                                                                          caaGCCTCGGCACATGCCGGGGCTTGGAATGGGTGGAGATGGAAACCCTGTCAGGTCTCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTCTTACCACTGCACGTTGACGCCGGGCTTTGATCATCTCAAAAAACGACatatcg      >  3:267001/2‑143 (MQ=255)                                                                                                                                          aaaGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGcc      >  3:46128/1‑148 (MQ=255)                                                                                                                                           aaGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCa    <  8:88332/149‑1 (MQ=255)                                                                                                                                             gCCTCGGCACATGCCGGGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGc  >  5:32932/1‑149 (MQ=255)                                                                                                                                             gCCTCGGCACATGCCGGGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGc  >  2:89360/1‑149 (MQ=255)                                                                                                                                             gCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGc  >  3:90253/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                  CGATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGC  >  NC_002947/4405173‑4405459

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