Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,284,555	A→G	6.0%	S165G (AGC→GGC)	PP_1121 →	OmpA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,284,555	0	A	G	6.0%	197.5 / 2.4	67	S165G (AGC→GGC)	PP_1121	OmpA family protein
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base G (2/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

AAGAAGTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTGGCAGACAACG  >  NC_002947/1284409‑1284699                                                                                                                                                   |                                                                                                                                                 aagaagTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGTGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGc                                                                                                                                                >  4:360280/1‑149 (MQ=255) aagaagTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGc                                                                                                                                                >  2:199057/1‑149 (MQ=255)   gaagTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGt                                                                                                                                              <  1:35997/149‑1 (MQ=255)     agTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGTGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCg                                                                                                                                            <  3:360280/149‑1 (MQ=255)     agTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCGCGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCg                                                                                                                                            <  4:279877/149‑1 (MQ=255)       tCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGt                                                                                                                                          <  5:275703/149‑1 (MQ=255)         aTTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCGTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGtt                                                                                                                                         <  2:58592/148‑1 (MQ=255)               aTCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAg                                                                                                                                  <  1:163071/149‑1 (MQ=255)               aTCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTCCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTc                                                                                                                                        >  2:230142/1‑143 (MQ=255)               aTCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTCCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTc                                                                                                                                        <  1:230142/143‑1 (MQ=255)               aTCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGAAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACa                                                                                                                                   <  1:141573/148‑1 (MQ=255)                 ccGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGc                                                                                                                                 >  7:337152/1‑148 (MQ=255)                                 gAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTc                                                                                                                <  8:337152/149‑1 (MQ=255)                                    ttCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTcc                                                                                                               >  3:359388/1‑147 (MQ=255)                                     tCGATTCAGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGa                                                                                                             <  6:308492/148‑1 (MQ=255)                                       gATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGa                                                                                                                       <  4:277212/136‑1 (MQ=255)                                       gATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGa                                                                                                                       >  3:277212/1‑136 (MQ=255)                                       gATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAgcg                                                                                                          >  4:146087/1‑149 (MQ=255)                                                   cTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTc                                                                                              >  2:226111/1‑149 (MQ=255)                                                      aCGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCg                                                                                             <  7:107228/147‑1 (MQ=255)                                                      aCGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGt                                                                                           <  3:146087/149‑1 (MQ=255)                                                         gCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGAc                                                                                        <  1:18744/149‑1 (MQ=255)                                                           cAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGAcc                                                                                       >  1:150727/1‑148 (MQ=255)                                                              tGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACGCTGTCGGtatgg                                                                                                                                      >  7:86171/1‑94 (MQ=255)                                                              tGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACGCTGTCGGtatgg                                                                                                                                      <  8:86171/98‑5 (MQ=255)                                                                  aaaGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGa                                                                                                                                     >  2:357207/1‑95 (MQ=255)                                                                  aaaGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGa                                                                                                                                     <  1:357207/95‑1 (MQ=255)                                                                   aaGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTg                                                                              >  2:212280/1‑149 (MQ=255)                                                                   aaGAGCGCCTCAATACCACTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTg                                                                              <  5:183911/149‑1 (MQ=255)                                                                    agagCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTgg                                                                             >  5:265277/1‑149 (MQ=255)                                                                         gcCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTAcc                                                                                <  4:76860/141‑1 (MQ=255)                                                                         gcCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTAcc                                                                                >  3:76860/1‑141 (MQ=255)                                                                          ccTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCTTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCga                                                                         >  6:34872/1‑147 (MQ=255)                                                                           cTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACGGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGagag                                                                      >  5:123654/1‑149 (MQ=255)                                                                              aaTACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCACTCGGTGACCGATTACCTGGTCGAGAGCg                                                                    <  4:359388/148‑1 (MQ=255)                                                                               aTACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGgt                                                                  <  6:265277/149‑1 (MQ=255)                                                                                tACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTAc                                                                                 >  6:182911/1‑133 (MQ=255)                                                                                tACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTAc                                                                                 <  5:182911/133‑1 (MQ=255)                                                                                   cATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTAc                                                              >  4:43983/1‑149 (MQ=255)                                                                                     ttGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGAc                                                                                                                                    >  2:164034/1‑77 (MQ=255)                                                                                        ccACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGACACTAACCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACcgcg                                                          >  1:118967/1‑148 (MQ=255)                                                                                         cACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCgca                                                        <  1:361478/149‑1 (MQ=255)                                                                                          aCGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCgcag                                                       <  3:43983/149‑1 (MQ=255)                                                                                              gcCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCagcagc                                                   <  1:199057/149‑1 (MQ=255)                                                                                                 tGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTc                                                >  6:68520/1‑149 (MQ=255)                                                                                                   aaGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTc                                                <  3:355370/147‑1 (MQ=255)                                                                                                   aaGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTc                                                <  1:226111/147‑1 (MQ=255)                                                                                                     gCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCga                                                                         <  6:200415/120‑1 (MQ=255)                                                                                                     gCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGttt                                            >  2:268303/1‑149 (MQ=255)                                                                                                       aGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCg                                          >  2:178175/1‑149 (MQ=255)                                                                                                           aGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACGGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCggtgg                                      <  6:123654/149‑1 (MQ=255)                                                                                                             cgcCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGgt                                                                  >  4:363809/1‑119 (MQ=255)                                                                                                             cgcCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGgt                                                                  <  3:363809/119‑1 (MQ=255)                                                                                                                   ctGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGcgcg                              >  2:211616/1‑149 (MQ=255)                                                                                                                      cccGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCgg                             <  5:68520/147‑1 (MQ=255)                                                                                                                      cccGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCgg                             <  1:178175/147‑1 (MQ=255)                                                                                                                       ccGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCGCATTCGGTGa                                                                                         <  4:11504/86‑1 (MQ=255)                                                                                                                                 gTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGcc                >  4:216740/1‑149 (MQ=255)                                                                                                                                  tCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTc                                    <  7:256628/128‑1 (MQ=255)                                                                                                                                  tCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTc                                    >  8:256628/1‑128 (MQ=255)                                                                                                                                    aGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGt             <  1:211616/149‑1 (MQ=255)                                                                                                                                        gCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTgg           >  6:110566/1‑147 (MQ=255)                                                                                                                                         ccATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTg            <  8:316438/145‑1 (MQ=255)                                                                                                                                         ccATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTg            >  7:316438/1‑145 (MQ=255)                                                                                                                                             aCCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTGGCAGCCaa    >  3:9949/1‑149 (MQ=255)                                                                                                                                             aCCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACAGATTACCTGGTCGAGAGAGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTGGCAGACaa    >  5:207813/1‑149 (MQ=255)                                                                                                                                               cGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTCGCAGACAACg  >  8:36112/1‑149 (MQ=255)                                                                                                                                                   aGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCggg                            >  4:109203/1‑119 (MQ=255)                                                                                                                                                   aGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCggg                            <  3:109203/119‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AAGAAGTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGACAGCGTCGGTTCCGACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTGGCAGACAACG  >  NC_002947/1284409‑1284699

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