Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,518,546	A→T	10.3%	I12N (ATT→AAT)	PP_4850 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,518,546	0	A	T	10.3%	181.7 / 10.4	68	I12N (ATT→AAT)	PP_4850	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base T (4/3);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

TTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGC  >  NC_002947/5518401‑5518689                                                                                                                                                   |                                                                                                                                                ttCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGAAGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGc                                                                                                                                              >  2:183823/1‑149 (MQ=255)   cGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGcc                                                                                                                                             >  2:97693/1‑148 (MQ=255)    gAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGc                                                                                                                                           <  8:226594/149‑1 (MQ=255)     aaCACCAGCAGACCGGCCTGCAGGCTTGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCa                                                                                                                                          <  3:24458/149‑1 (MQ=255)     aaCACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCg                                                                                                                                            >  1:51582/1‑147 (MQ=255)       cacCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGt                                                                                                                                        <  1:325347/149‑1 (MQ=255)        acCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTc                                                                                                                                       >  5:101953/1‑149 (MQ=255)          cagcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTg                                                                                                                                     <  6:101953/149‑1 (MQ=255)          cagcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGCCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGAGTATACAGCCTCAACGCCACACAGCGCGAATGCCGCaaaat                                                                                                                                      >  1:119794/1‑145 (MQ=255)           agcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGc                                                                                                                                    <  5:72283/149‑1 (MQ=255)           agcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGc                                                                                                                                    >  8:28448/1‑149 (MQ=255)             cagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCa                                                                                                                                  <  7:28448/149‑1 (MQ=255)                    tccAGCAGGGTGGTGATCAACTTCAGGGTCAAAGCGAGCAGCAGGGCGGCGATGACGATGATCAGGTTGGCGAGGGTGATCTGTCTGTGGGTAAAGG‑TCATTCAGCCTCGACGGCAGACAGCTGTAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                           <  5:361345/148‑1 (MQ=255)                      cagcagGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTggcg                                                                                                                         <  1:45229/149‑1 (MQ=255)                              tGGGCAGCAACTTCAGGTGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGCCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                           <  5:182186/139‑1 (MQ=255)                              tGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGCCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                           >  6:182186/1‑139 (MQ=255)                               gggCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCa                                                                                                                >  4:297250/1‑149 (MQ=255)                                  cagcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGAGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAg                                                                                                             >  4:293446/1‑149 (MQ=255)                                   agcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAgg                                                                                                            <  4:82951/149‑1 (MQ=255)                                   agcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAgg                                                                                                            >  4:204075/1‑149 (MQ=255)                                      aaCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGccgtc                                                                                                                                <  2:272953/126‑4 (MQ=38)                                      aaCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGccgtc                                                                                                                                >  1:272953/1‑123 (MQ=38)                                       aCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                   <  5:296231/138‑1 (MQ=255)                                       aCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                   >  6:296231/1‑138 (MQ=255)                                           cAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCa                                                                                                        <  5:93106/145‑1 (MQ=255)                                           cAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCa                                                                                                        >  6:93106/1‑145 (MQ=255)                                               ggCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGTAGGTCAttt                                                                                                >  5:285884/1‑149 (MQ=255)                                                gCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttc                                                                                               <  1:183823/149‑1 (MQ=255)                                                gCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttc                                                                                               >  1:149683/1‑149 (MQ=255)                                                   aaGCCAGCAGCAGCCCGGCCATTACGATGATCAGGCTGGCGAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGAGGGCAGAGAGCGGTAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttc                                                                                            <  1:97693/149‑1 (MQ=255)                                                   aaGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttc                                                                                            <  3:11807/149‑1 (MQ=255)                                                     gCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGTCGGCAGACAGCGGGAAGGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttctt                                                                                          <  7:157570/149‑1 (MQ=255)                                                       cagcagcagCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                        <  2:51582/149‑1 (MQ=255)                                                       cagcagcagCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                        <  1:8709/149‑1 (MQ=255)                                                         gcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGTAGGTCATTTCTTCTTCAgg                                                                                      <  6:285884/149‑1 (MQ=255)                                                         gcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAgg                                                                                      <  3:297250/149‑1 (MQ=255)                                                             cagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttctt                                                                                             <  2:34218/138‑1 (MQ=255)                                                             cagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttctt                                                                                             >  1:34218/1‑138 (MQ=255)                                                               gcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCg                                                                                >  7:146843/1‑149 (MQ=255)                                                                cccGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGc                                                                               >  8:70816/1‑149 (MQ=255)                                                                   ggCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCa                                                                            <  2:149683/149‑1 (MQ=255)                                                                    gCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCAt                                                                           >  7:92053/1‑149 (MQ=255)                                                                              tgatCTGGCTGGCCAGGGTGTTCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCTGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTg                                                                 <  2:119794/149‑1 (MQ=255)                                                                               gatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgc                                                                >  3:76914/1‑149 (MQ=255)                                                                                atCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgc                                                                <  8:92053/148‑1 (MQ=255)                                                                                 tCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgcgc                                                              <  8:146843/149‑1 (MQ=255)                                                                                            aGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCtt                                                   >  1:173368/1‑149 (MQ=255)                                                                                               gTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAg                                                >  8:118308/1‑149 (MQ=255)                                                                                                 gATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGc                                                                     >  8:210884/1‑126 (MQ=255)                                                                                                 gATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGc                                                                     <  7:210884/126‑1 (MQ=255)                                                                                                    cTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                        >  8:274989/1‑104 (MQ=255)                                                                                                    cTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                        <  7:274989/104‑1 (MQ=255)                                                                                                     tGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCttt                                          >  3:20985/1‑149 (MQ=255)                                                                                                         ggggggTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtg                                      <  7:118308/149‑1 (MQ=255)                                                                                                           ggggTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggt                                    <  4:76914/149‑1 (MQ=255)                                                                                                              gTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGc                                 <  2:173368/149‑1 (MQ=255)                                                                                                                   gg‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAg                            >  8:88393/1‑149 (MQ=255)                                                                                                                           cAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGc                      <  4:20985/148‑1 (MQ=255)                                                                                                                            aGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggtgg                                  >  4:102302/1‑135 (MQ=255)                                                                                                                            aGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggtgg                                  <  3:102302/135‑1 (MQ=255)                                                                                                                                   ctggCAGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCa             >  8:225392/3‑149 (MQ=255)                                                                                                                                      gCAGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGGGGTGGCGTCAGCATTGGAGGTATTCaata          >  3:26475/1‑146 (MQ=255)                                                                                                                                        aGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtgg                                        >  2:126337/1‑117 (MQ=255)                                                                                                                                        aGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtgg                                        <  1:126337/117‑1 (MQ=255)                                                                                                                                         gACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAtttttt       >  6:95453/1‑149 (MQ=255)                                                                                                                                            aGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCtt    <  3:204075/149‑1 (MQ=255)                                                                                                                                              cGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGc  >  7:16252/1‑149 (MQ=255)                                                                                                                                              cGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGc  >  8:288163/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                TTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGC  >  NC_002947/5518401‑5518689

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