Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,405,321	G→A	10.0%	H80H (CAC→CAT)	PP_3886 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,405,321	0	G	A	10.0%	164.1 / 6.5	61	H80H (CAC→CAT)	PP_3886	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/25);  new base A (4/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTG  >  NC_002947/4405173‑4405465                                                                                                                                                       |                                                                                                                                                 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                                                                                   <  4:95050/140‑1 (MQ=255)                                                                         aCTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                                   <  4:89469/140‑1 (MQ=255)                                                                         aCTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                                   >  3:95050/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                                          <  1:401252/149‑1 (MQ=255)                                                                                                           aCGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCgg                                          >  7:88332/1‑149 (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                                  >  8:304621/1‑135 (MQ=255)                                                                                                                                 tttAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGat                                  <  7:304621/135‑1 (MQ=255)                                                                                                                                     aCGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACa                 >  1:146288/1‑148 (MQ=255)                                                                                                                                       gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg              >  6:404080/1‑149 (MQ=255)                                                                                                                                       gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg              >  3:191196/1‑149 (MQ=255)                                                                                                                                       gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg              >  5:3883/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                                                >  8:129917/1‑103 (MQ=255)                                                                                                                                                   gCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTgc                                                <  7:129917/103‑1 (MQ=255)                                                                                                                                                   gCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTg  <  8:50161/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 CGATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCT‑C‑GACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTG  >  NC_002947/4405173‑4405465

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