Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	668,266	G→T	5.5%	F87L (TTC→TTA) ‡	PP_0571 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	668,266	0	G	T	5.5%	294.9 / 3.0	91	F87L (TTC→TTA) ‡	PP_0571	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (45/41);  new base T (2/3);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACC  >  NC_002947/668119‑668405                                                                                                                                                    |                                                                                                                                            cgAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGaa                                                                                                                                           <  1:575913/150‑1 (MQ=255)  gAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAAt                                                                                                                                          >  1:94198/1‑150 (MQ=255)    aGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATgc                                                                                                                                        >  1:728224/1‑150 (MQ=255)        gAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATTCATGcgcgc                                                                                                                                    <  1:1389502/150‑2 (MQ=255)         aaaGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGa                                                                                                                                   >  1:585643/1‑150 (MQ=255)         aaaGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGa                                                                                                                                   >  2:533599/1‑150 (MQ=255)              cTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATTCATGcgcgc                                                                                                                                    >  2:1391380/1‑143 (MQ=255)              cTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATTCATGcgcgc                                                                                                                                    <  1:1391380/144‑2 (MQ=255)              cTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGcc                                                                                                                              >  2:585176/1‑150 (MQ=255)              cTGGGGCTGCTGTGCAGCACTGTCGCGCCGCTGGGCAGCAGCTGACCTGACCGTACTTATGAGTGGATCATGCAGCTCATCACCCACAGCCCCAGCACCAGCCAGCTGATTCCGCGGCTGATCGCGGCGCGTCTGTATGCGCTTATTgtt                                                                                                                              <  2:766319/150‑3 (MQ=255)              cTGGCGCGGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCACACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGcc                                                                                                                              <  2:1432463/150‑1 (MQ=255)                ggggCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGa                                                                                                                            <  1:404045/150‑1 (MQ=255)                   gctgctGTGCAGACCAGTGAGAACTCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTa                                                                                                                         <  2:94198/150‑1 (MQ=255)                      gctgTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTAcag                                                                                                                      <  2:1516652/150‑1 (MQ=255)                          tgCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGAtt                                                                                                                                 <  1:1159776/135‑1 (MQ=255)                          tgCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGAtt                                                                                                                                 >  2:1159776/1‑135 (MQ=255)                           gCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAgcagca                                                                                                                 >  1:117018/1‑150 (MQ=255)                              gCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCa                                                                                                              <  2:308531/150‑1 (MQ=255)                                           gCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAg                                                                                                 >  1:1501070/1‑150 (MQ=255)                                            cAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAgg                                                                                                <  1:1004703/150‑1 (MQ=255)                                             aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCCGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTAca                                                                                                                       <  2:32720/126‑1 (MQ=255)                                             aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTAca                                                                                                                       >  1:32720/1‑126 (MQ=255)                                              aGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGc                                                                                              <  2:117018/150‑1 (MQ=255)                                                       ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAAt                                                                                                                                          <  2:294910/97‑1 (MQ=255)                                                       ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAAt                                                                                                                                          >  1:294910/1‑97 (MQ=255)                                                       ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGcgcg                                                                                                                                     <  1:295874/102‑1 (MQ=255)                                                       ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGcgcg                                                                                                                                     >  2:295874/1‑102 (MQ=255)                                                       ccACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCt                                                                                     >  1:924190/1‑150 (MQ=255)                                                        cacaAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGCTGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTc                                                                                    >  2:532479/1‑150 (MQ=255)                                                        cacaAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTc                                                                                    >  2:572447/1‑150 (MQ=255)                                                        cacaAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTc                                                                                    >  1:634010/1‑150 (MQ=255)                                                           aaaGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTcagc                                                                                 >  1:199581/1‑150 (MQ=255)                                                             aGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCgatg                                                                                                      >  2:1026970/1‑127 (MQ=255)                                                             aGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCgatg                                                                                                      <  1:1026970/127‑1 (MQ=255)                                                              gACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCagcagc                                                                              <  2:924190/150‑1 (MQ=255)                                                                 cgcgCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGaa                                                                           <  1:533599/150‑1 (MQ=255)                                                                 cgcgCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGaa                                                                           <  2:199581/150‑1 (MQ=255)                                                                  gcgcAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATc                                                                                                               <  2:765530/113‑1 (MQ=255)                                                                  gcgcAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATc                                                                                                               >  1:765530/1‑113 (MQ=255)                                                                   cgcAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACagcagc                                                                                                                  >  2:764936/1‑109 (MQ=255)                                                                   cgcAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACagcagc                                                                                                                  <  1:764936/109‑1 (MQ=255)                                                                   cgcAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAgt                                                                         <  1:585176/150‑1 (MQ=255)                                                                    gcAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAgtg                                                                        >  1:839127/1‑150 (MQ=255)                                                                        aTCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCa                                                                    >  2:160748/1‑150 (MQ=255)                                                                        aTCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCa                                                                    >  2:672934/1‑150 (MQ=255)                                                                           aGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAgg                                                                                          >  2:1145566/1‑125 (MQ=255)                                                                           aGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAgg                                                                                          <  1:1145566/125‑1 (MQ=255)                                                                            gTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGcc                                                                <  1:324587/150‑1 (MQ=255)                                                                               gATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCaa                                                             >  2:725592/1‑150 (MQ=255)                                                                                aTCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAg                                                            <  1:672934/150‑1 (MQ=255)                                                                                    cccAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGGTGATCGATGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGAAGATCAGGGCCAGGGTGCTCAGCAGCGAAGTGTCCATGCGCAAGGTgc                                                        >  1:378686/1‑150 (MQ=255)                                                                                     ccAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAg                                                                                                 <  1:1183736/108‑1 (MQ=255)                                                                                     ccAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAg                                                                                                 >  2:1183736/1‑108 (MQ=255)                                                                                     ccAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGaa                                                                           >  2:969217/1‑130 (MQ=255)                                                                                     ccAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGaa                                                                           <  1:969217/130‑1 (MQ=255)                                                                                     ccAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTg                                                         >  1:398549/1‑148 (MQ=255)                                                                                     ccAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTg                                                         <  2:398549/148‑1 (MQ=255)                                                                                      cAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGCCCATGCCCAAGGTgcgc                                                      <  2:839127/150‑1 (MQ=255)                                                                                          tcatcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACa                                                  <  1:52966/150‑1 (MQ=255)                                                                                           catcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAg                                                 <  2:1501070/150‑1 (MQ=255)                                                                                            atcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGTTTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTca                                                                                   <  2:638881/115‑1 (MQ=255)                                                                                            atcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTca                                                                                   >  1:638881/1‑115 (MQ=255)                                                                                             tcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGcc                                                                                                                              <  2:1453471/71‑1 (MQ=255)                                                                                             tcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGcc                                                                                                                              >  1:1453471/1‑71 (MQ=255)                                                                                             tcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGc                                               >  2:1044217/1‑150 (MQ=255)                                                                                               aCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCa                                             >  2:1286532/1‑150 (MQ=255)                                                                                                 ccACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTAAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATc                                           >  1:582978/1‑150 (MQ=255)                                                                                                      gCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGa                                      >  1:763439/1‑150 (MQ=255)                                                                                                       ccccAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGaa                                     <  1:572447/150‑1 (MQ=255)                                                                                                       ccccAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGaa                                     >  1:1166233/1‑150 (MQ=255)                                                                                                        cccAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTcag                                                                                  <  1:573703/104‑1 (MQ=255)                                                                                                        cccAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTcag                                                                                  >  2:573703/1‑104 (MQ=255)                                                                                                        cccAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAg                                    >  2:261885/1‑150 (MQ=255)                                                                                                                  gccagATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCgg                          >  2:1048946/1‑150 (MQ=255)                                                                                                                  gccagATGATTCCGAGGATGATCGAGGCGCGCATGAACCACCGGCCAGAACTATTTAGCAACATCAGGCCGAGTATCAGGGCCAGGATGCTCAGGAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCgg                          >  1:445068/1‑150 (MQ=255)                                                                                                                   ccagATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGaa                                                                           <  2:1614330/100‑1 (MQ=255)                                                                                                                   ccagATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGaa                                                                           >  1:1614330/1‑100 (MQ=255)                                                                                                                    cagATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTgcgc                                                      <  2:1577545/120‑1 (MQ=255)                                                                                                                    cagATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTgcgc                                                      >  1:1577545/1‑120 (MQ=255)                                                                                                                    cagATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCa                        >  1:1210228/1‑150 (MQ=255)                                                                                                                        tgatTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTgg                    <  1:1048946/150‑1 (MQ=255)                                                                                                                          atTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGcc                  <  1:532479/150‑1 (MQ=255)                                                                                                                           ttCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCca                 <  2:585643/150‑1 (MQ=255)                                                                                                                           ttCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTAAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCca                 <  2:582978/150‑1 (MQ=255)                                                                                                                                         gAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGc                   >  1:1476320/1‑134 (MQ=255)                                                                                                                                         gAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGc                   <  2:1476320/134‑1 (MQ=255)                                                                                                                                         aatccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGaac      <  2:852438/143‑1 (MQ=255)                                                                                                                                         aatccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGaac      >  1:852438/5‑147 (MQ=255)                                                                                                                                          aGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAgtgt                                                                       <  1:937427/81‑1 (MQ=255)                                                                                                                                          aGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAgtgt                                                                       >  2:937427/1‑81 (MQ=255)                                                                                                                                          aGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAAcc  >  1:93458/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                            CGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACC  >  NC_002947/668119‑668405

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑