Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,405,327	A→G	7.8%	P78P (CCT→CCC)	PP_3886 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,405,327	0	A	G	7.8%	181.5 / 4.5	64	P78P (CCT→CCC)	PP_3886	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base G (3/2);  total (34/30)
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.00e+00

AGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCT  >  NC_002947/4405181‑4405473                                                                                                                                                     |                                                                                                                                                   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                                                                                           <  4:95050/140‑1 (MQ=255)                                                                 aCTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                                           <  4:89469/140‑1 (MQ=255)                                                                 aCTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                                           >  3:95050/1‑140 (MQ=255)                                                                 aCTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                                           >  3:89469/1‑140 (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                                                                                                    <  1:164587/90‑1 (MQ=255)                                                                                                            aTGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGgc                                                                                                    >  2:164587/1‑90 (MQ=255)                                                                                                                cTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCTCCGGGCTTTGATCAtct                                     <  6:357562/149‑1 (MQ=255)                                                                                                                         tttAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGat                                          >  8:304621/1‑135 (MQ=255)                                                                                                                         tttAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGat                                          <  7:304621/135‑1 (MQ=255)                                                                                                                             aCGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACa                         >  1:146288/1‑148 (MQ=255)                                                                                                                               gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg                      >  5:3883/1‑149 (MQ=255)                                                                                                                               gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg                      >  3:191196/1‑149 (MQ=255)                                                                                                                               gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg                      >  6:404080/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                >  2:89360/1‑149 (MQ=255)                                                                                                                                     gCCTCGGCACATGCCGGGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGc                >  5:32932/1‑149 (MQ=255)                                                                                                                                     gCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGc                >  3:90253/1‑149 (MQ=255)                                                                                                                                           gCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTgc                                                        <  7:129917/103‑1 (MQ=255)                                                                                                                                           gCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTgc                                                        >  8:129917/1‑103 (MQ=255)                                                                                                                                           gCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTg          <  8:50161/149‑1 (MQ=255)                                                                                                                                                tgCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCg      <  7:280845/148‑1 (MQ=255)                                                                                                                                                tgCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCg      <  2:353344/148‑1 (MQ=255)                                                                                                                                                   cGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCt  >  5:17457/1‑149 (MQ=255)                                                                                                                                                   cGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCt  <  4:95060/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGGCT  >  NC_002947/4405181‑4405473

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