Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,166,837	G→A	5.9%	A213A (GCC→GCT)	creA ←	creatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,166,837	0	G	A	5.9%	205.9 / 2.2	68	A213A (GCC→GCT)	creA	creatinase
Reads supporting (aligned to +/- strand):  ref base G (35/29);  new base A (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

TCCAGCGCGGTGTAATAACCGGCGATCATGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCT  >  NC_002947/4166691‑4166985                                                                                                                                                    |                                                                                                                                                       tCCAGCGCGGTGTAATAACCGGCGATCATGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACTGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGGCCATCAGCTCGGc                                                                                                                                                      >  2:237173/1‑149 (MQ=255)        cgGTGTAATAACCGGCGATCATGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCgggg                                                                                                                                               <  5:273598/149‑1 (MQ=255)         ggTGTAATAACCGGCGATCATGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCgggg                                                                                                                                               >  4:206840/1‑148 (MQ=255)              aataaCCGGCGATCATGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGt                                                                                                                                         >  8:40545/1‑149 (MQ=255)               ataaCCGGCGATCATGGGGAAACAGTTGAGCCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGtt                                                                                                                                        >  1:28299/1‑149 (MQ=255)                taaCCGGCGATCATGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGttt                                                                                                                                       <  8:88741/149‑1 (MQ=255)                taaCCGGCGATCATGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGttt                                                                                                                                       <  3:206840/149‑1 (MQ=255)                    cGGCGATCATGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGc                                                                                                                                    >  5:22309/1‑148 (MQ=255)                    cGGCGATCATGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGc                                                                                                                                    >  7:95754/1‑148 (MQ=255)                      gCGATCATGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAAt                                                                                                                                 >  6:165556/1‑149 (MQ=255)                            aTGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTcgc                                                                                                                            <  3:119469/148‑1 (MQ=255)                             tGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCgcgc                                                                                                                          <  6:22309/149‑1 (MQ=255)                                 gAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCa                                                                                                                      >  4:228204/1‑149 (MQ=255)                                       gTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCt                                                                                                                <  2:28299/149‑1 (MQ=255)                                           aGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCCGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTGTGCAATTTCGCGCACCATCGCCT‑GGGt                                                                                                           >  3:177118/1‑149 (MQ=255)                                            ggCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTg                                                                                                          <  2:137005/149‑1 (MQ=255)                                            ggCTGAGGATCTCGCCCTTCCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTg                                                                                                          >  7:184935/1‑149 (MQ=255)                                               tGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCa                                                                                                       >  7:149623/1‑149 (MQ=255)                                                    aTCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTCTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAg                                                                                                 >  2:111815/1‑149 (MQ=255)                                                         gCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGc                                                                                              <  3:314199/147‑1 (MQ=255)                                                         gCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCa                                                                                            <  3:118667/149‑1 (MQ=255)                                                           ccTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCACTGCA‑TGCCAcc                                                                                          <  6:146935/149‑1 (MQ=255)                                                               gCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCACTGCA‑TGCCAccca                                                                                        <  5:151763/147‑3 (MQ=255)                                                                 ccACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGCGCCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGt                                                                                      >  1:233317/1‑148 (MQ=255)                                                                   acacGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGtt                                                                                                                                        <  7:101374/97‑1 (MQ=255)                                                                   acacGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGtt                                                                                                                                        >  8:101374/1‑97 (MQ=255)                                                                    cacGGCGGCTGGTGTCCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCg                                                                                 <  8:149623/149‑1 (MQ=255)                                                                     acGGCGGCTGGTGACCGGGTGGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCTGGGAAGCGTTTGGCAATTTCGCGCACCATCGGCT‑GTGTGGCA‑TTCAGTGCCACCTCGTACTCgg                                                                                <  6:233217/149‑1 (MQ=255)                                                                          ggcTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCAcc                                                                           >  5:183676/1‑149 (MQ=255)                                                                          ggcTGGTGACCGGGTTCTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCgacggt                                                                            >  2:268740/1‑143 (MQ=255)                                                                             tGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGc                                                                        >  6:272502/1‑149 (MQ=255)                                                                                   ccGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGc                                                                  <  3:228204/149‑1 (MQ=255)                                                                                    cGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACg                                                                   >  6:10565/1‑147 (MQ=255)                                                                                    cGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCAGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGCGCa                                                                             <  8:147070/137‑1 (MQ=255)                                                                                    cGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCAGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGCGCa                                                                             >  7:147070/1‑137 (MQ=255)                                                                                     gggTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAg                                                                <  7:40545/149‑1 (MQ=255)                                                                                      ggTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAgg                                                               <  5:10565/149‑1 (MQ=255)                                                                                             tgCGCCATCGGTAGTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCg                                                         <  1:111815/148‑1 (MQ=255)                                                                                               cgcCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGAc                                                       <  5:165556/148‑1 (MQ=255)                                                                                                    tCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTgc                                                 >  6:50685/1‑149 (MQ=255)                                                                                                      ggTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTcgcg                                                                                                                           >  7:50024/1‑75 (MQ=255)                                                                                                      ggTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTcgcg                                                                                                                           <  8:50024/75‑1 (MQ=255)                                                                                                             aGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGccgccccgcc                                        >  8:74419/1‑149 (MQ=255)                                                                                                                    gACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTc                                  <  5:50685/149‑1 (MQ=255)                                                                                                                      cTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCgg                                >  3:145927/1‑149 (MQ=255)                                                                                                                       tGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCATTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGc                               >  6:254556/1‑149 (MQ=255)                                                                                                                         gAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCa                             >  2:46795/1‑149 (MQ=255)                                                                                                                          aaCCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCCCCATCGCCT‑GGGTGGCA‑TGCAGTGGCAGCTCGTACTCGGGCCCCTGCTCACGCAGGGCTTCGCGCACTGCCGCCCCGCCGATGTCGGCCAc                            >  8:79371/1‑149 (MQ=255)                                                                                                                          aaCCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTAGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCAc                            >  7:262162/1‑149 (MQ=255)                                                                                                                            ccaggtccaggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCt                          <  5:272502/149‑1 (MQ=255)                                                                                                                             caggtccaggtGTCCATCAGCTCAGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCtg                         >  5:42725/1‑149 (MQ=255)                                                                                                                                 tccaggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCa                     >  3:255756/1‑149 (MQ=255)                                                                                                                                   caggtGTCAAACAGCTCGGCGTATGGGAAGCGATAGGCAAATTCGCCCACCATAGCCA‑GCGATGCA‑TCCAGTGCCAACTCGTACTCAGGCACCTGCTCACGCAGCGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCAcc                   >  2:60563/1‑149 (MQ=255)                                                                                                                                    aggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCAc                                                                    <  4:261261/99‑1 (MQ=255)                                                                                                                                    aggtGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCAc                                                                    >  3:261261/1‑99 (MQ=255)                                                                                                                                      gtgtCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCTAGGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGt                 >  1:89111/1‑149 (MQ=255)                                                                                                                                       tgtCCATCAGCTCAGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGa               <  6:42725/149‑1 (MQ=255)                                                                                                                                         tCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGAc              <  3:222490/148‑1 (MQ=255)                                                                                                                                         tCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACg             >  2:44800/1‑149 (MQ=255)                                                                                                                                           cATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACgg            <  4:145927/148‑1 (MQ=255)                                                                                                                                           cATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGa           >  1:292508/1‑149 (MQ=255)                                                                                                                                            aTCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GTGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCATGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGAt          <  1:46795/149‑1 (MQ=255)                                                                                                                                            aTCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGAt          <  1:237173/149‑1 (MQ=255)                                                                                                                                             tCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGAt          >  1:168653/1‑148 (MQ=255)                                                                                                                                               aGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATaa        <  6:183676/148‑1 (MQ=255)                                                                                                                                                gCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGt      <  1:44800/149‑1 (MQ=255)                                                                                                                                                  tCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGc    >  8:92960/1‑149 (MQ=255)                                                                                                                                                   cGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGcc   >  2:19878/1‑149 (MQ=255)                                                                                                                                                    ggCGTCGGGGAAGCGTTTGGCAATTTCGCGCACGATCGACT‑GGGTGGGA‑TGCAGTGCCACCGCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGc    <  4:177118/147‑1 (MQ=255)                                                                                                                                                    ggCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCt  <  2:233317/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       TCCAGCGCGGTGTAATAACCGGCGATCATGGGGAAACAGTTGAGGCTGAGGATCTCGCCCTTGCCCACACGGCGGCTGGTGACCGGGTTGTGTGCGCCATCGGTGTTGAGC‑CCGGACTGGAACCAGGTCCAGGTGTCCATCAGCTCGGCGTCGGGGAAGCGTTTGGCAATTTCGCGCACCATCGCCT‑GGGTGGCA‑TGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCT  >  NC_002947/4166691‑4166985

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