Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,183,003	A→T	14.7%	I244N (ATC→AAC)	PP_4563 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,183,003	0	A	T	14.7%	170.2 / 18.1	68	I244N (ATC→AAC)	PP_4563	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/24);  new base T (3/7);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTT  >  NC_002947/5182855‑5183151                                                                                                                                                     |                                                                                                                                                       cGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                                        <  2:276557/149‑1 (MQ=255)      agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                                   <  7:246038/149‑1 (MQ=255)       ggggtaCTTCCGCAACTTCTTACAGCCGTAAGGCCGTAAGGC‑TTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                                   <  1:132889/142‑1 (MQ=255)        cTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgt                                                                                                                                                 <  2:370074/149‑2 (MQ=255)           gACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgtca                                                                                                                                               <  3:207730/148‑4 (MQ=255)           gACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgtc                                                                                                                                                <  4:64707/147‑3 (MQ=255)            aCTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                                   <  1:366241/143‑1 (MQ=255)            aCTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                                   >  2:366241/1‑143 (MQ=255)            aCTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAggg                                                                                                                                             <  6:362455/149‑1 (MQ=255)                 cGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑ct                                                                                                                                       >  2:67899/1‑149 (MQ=255)                         cGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTc                                                                                                                                 <  3:317918/147‑1 (MQ=255)                                gCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACAcccc                                                                                                                        >  1:332797/1‑149 (MQ=255)                                      ggCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAacac                                                                                                                  >  5:206240/1‑149 (MQ=255)                                        ccGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACAccc                                                                                                                >  8:146582/1‑149 (MQ=255)                                             aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCt                                                                                                           >  2:66381/1‑149 (MQ=255)                                               gCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTc                                                                                                         <  8:242853/149‑1 (MQ=255)                                                 tttCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCgg                                                                                                       >  1:73111/1‑149 (MQ=255)                                                   tctcGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAAACCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTCCTGATCAACGAAGGGTCA‑ATCTTGTCGAACACCCCCAACACCACGCATTctcat                                                                                                     >  7:79413/1‑145 (MQ=255)                                                     tcGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTc                                                                                                   >  7:257690/1‑149 (MQ=255)                                                     tcGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTc                                                                                                   >  3:95368/1‑149 (MQ=255)                                                        aTGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAg                                                                                                <  7:146582/149‑1 (MQ=255)                                                                gTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCa                                                                                        >  2:81930/1‑149 (MQ=255)                                                                 tCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTACCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAt                                                                                       >  1:286638/1‑149 (MQ=255)                                                                   aTGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑Ga                                                                                    <  2:195354/149‑1 (MQ=255)                                                                            gATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCtt                                                                           >  6:228066/1‑149 (MQ=255)                                                                            gATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCtt                                                                           <  6:206240/149‑1 (MQ=255)                                                                              tCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTg                                                                          >  4:271161/1‑148 (MQ=255)                                                                                aGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCa                                                                        <  2:286638/148‑1 (MQ=255)                                                                                   cccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGgt                                                                    <  8:257690/149‑1 (MQ=255)                                                                                    ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGgtg                                                                   >  6:204704/1‑149 (MQ=255)                                                                                        cTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTcc                                                                <  5:228066/148‑1 (MQ=255)                                                                                        cTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTcc                                                                <  1:66381/148‑1 (MQ=255)                                                                                         ttGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGc                                                              >  2:185554/1‑149 (MQ=255)                                                                                           gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑ctc                                                                                                                                      >  6:3378/1‑76 (MQ=255)                                                                                           gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑ctc                                                                                                                                      <  5:3378/76‑1 (MQ=255)                                                                                            aacaacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAg                                                           >  4:186195/1‑149 (MQ=255)                                                                                             acaacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCTCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAgg                                                          <  1:67899/149‑1 (MQ=255)                                                                                             acaacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAgg                                                          <  4:95368/149‑1 (MQ=255)                                                                                              caacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAggg                                                         <  2:73111/149‑1 (MQ=255)                                                                                               aacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGggc                                                        >  5:291708/1‑149 (MQ=255)                                                                                                    gCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGgcgc                                                   >  6:151458/1‑149 (MQ=255)                                                                                                    gCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGgcgc                                                   >  8:187869/1‑149 (MQ=255)                                                                                                        aTGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTAc                                               >  4:258085/1‑149 (MQ=255)                                                                                                         tGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTAcc                                              <  5:151458/149‑1 (MQ=255)                                                                                                             cttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTc                                                                                                         <  3:90768/87‑1 (MQ=255)                                                                                                             cttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTc                                                                                                         >  4:90768/1‑87 (MQ=255)                                                                                                              ttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTg                                         >  2:60499/1‑149 (MQ=255)                                                                                                                   gAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATCGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGc                                        <  1:81930/146‑1 (MQ=255)                                                                                                                   gAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGc                                    >  3:7412/1‑149 (MQ=255)                                                                                                                    aaGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGc                                        >  4:124079/1‑144 (MQ=255)                                                                                                                    aaGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGc                                        <  3:124079/144‑1 (MQ=255)                                                                                                                      gTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCt                                   >  1:283353/1‑147 (MQ=255)                                                                                                                       ttGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAg                                <  7:187869/149‑1 (MQ=255)                                                                                                                           ccTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAg                                      >  2:237450/1‑139 (MQ=255)                                                                                                                           ccTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAg                                      <  1:237450/139‑1 (MQ=255)                                                                                                                            cTTGCTGACGTTGAACAGCTTGCTGATCAACGCTGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGaaa                           >  8:250783/1‑149 (MQ=255)                                                                                                                              tGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATc                         >  2:21347/1‑149 (MQ=255)                                                                                                                                cctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCAGCTCTTGTCGCACCCCCCCAACCCCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCt                        >  1:59912/5‑149 (MQ=255)                                                                                                                                cTGACGTTGTACAGATTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCTGGCTCGCGAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCTGCGCCTACCGGCTGCAGGCTCAGTGAAATCtt                       <  8:79413/149‑1 (MQ=255)                                                                                                                                 tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGTCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg                      >  6:44524/1‑149 (MQ=255)                                                                                                                                 tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg                      >  6:374047/1‑149 (MQ=255)                                                                                                                                 ctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg                      >  2:361452/4‑149 (MQ=255)                                                                                                                                   aCGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg                      <  2:332797/147‑1 (MQ=255)                                                                                                                                    cGTTGATCAGCAAGCTGTTCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGcc                    <  1:361452/148‑1 (MQ=255)                                                                                                                                        gAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCt               >  2:364939/1‑149 (MQ=255)                                                                                                                                        gAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCt               >  1:182824/1‑149 (MQ=255)                                                                                                                                          aCAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTgg             <  2:283353/149‑1 (MQ=255)                                                                                                                                              cTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCa         <  1:21347/149‑1 (MQ=255)                                                                                                                                              cTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGAGCGCCAT‑GATCTGCTCTTGCAGGCTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGATCAGTGAAATCTTGCCGTCATGGCGCa         >  7:336004/1‑149 (MQ=255)                                                                                                                                                     aTCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCtt  >  3:25062/1‑149 (MQ=255)                                                                                                                                                     aTCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCt   <  3:258085/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       CGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCTT  >  NC_002947/5182855‑5183151

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