Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,707,048	A→C	10.0%	intergenic (+17/+30)	PP_2369 → / ← PP_2370	leucine‑rich repeat‑containing protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,707,048	0	A	C	10.0%	168.6 / 10.7	60	intergenic (+17/+30)	PP_2369/PP_2370	leucine‑rich repeat‑containing protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/24);  new base C (2/4);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

CAATCGAAAACGCTGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCAATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCGGCACGGATG  >  NC_002947/2706903‑2707186                                                                                                                                                  |                                                                                                                                           cAATCGAAAACGCTGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGagag                                                                                                                                         >  7:331186/1‑149 (MQ=255)  aaTCGAAAACGCTGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAgg                                                                                                                                        >  6:114495/1‑149 (MQ=255)      gAAAACGCTGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCggg                                                                                                                                    <  8:331186/149‑1 (MQ=255)      gAAAACGCTGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCAATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCggg                                                                                                                                    <  4:161322/149‑1 (MQ=255)        aaaCGCTGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCAATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCTCTCGGCGGGCt                                                                                                                                  <  6:322166/149‑1 (MQ=255)              tGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGagag                                                                                                                                         <  8:7515/136‑1 (MQ=255)              tGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGagag                                                                                                                                         <  8:2275/136‑1 (MQ=255)              tGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGagag                                                                                                                                         >  7:7515/1‑136 (MQ=255)              tGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGagag                                                                                                                                         >  7:2275/1‑136 (MQ=255)              tGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGc                                                                                                                            >  1:18626/1‑149 (MQ=255)                aGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACAGCCTGTGAACGAAACAGCCCGCCTCTCGGCGGGCTGtttcg                                                                                                                            <  4:176835/147‑4 (MQ=255)                  tGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCTCTCGGCGGGCTGtttcgttca                                                                                                                        <  8:116161/149‑8 (MQ=255)                       gCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCgg                                                                                                                    >  4:217641/1‑148 (MQ=255)                         gACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCAATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGAc                                                                                                                 >  5:157505/1‑149 (MQ=255)                         gACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGACAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGAc                                                                                                                 >  1:94637/1‑149 (MQ=255)                            tATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGAt                                                                                                              >  8:138959/1‑149 (MQ=255)                                    ccAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGCTCGCTGGGGCGGGACGATCAgcggcg                                                                                                      <  1:242316/149‑1 (MQ=255)                                           gAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGaa                                                                                                >  4:173757/1‑148 (MQ=255)                                           gAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGaa                                                                                                >  4:140201/1‑148 (MQ=255)                                                tatCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCAATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAAcagcgg                                                                                          >  5:174104/1‑149 (MQ=255)                                                   cTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAAcagcggcag                                                                                       >  4:263637/1‑149 (MQ=255)                                                           cgcgTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTtcgt                                                                               <  3:263637/149‑1 (MQ=255)                                                                   gTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTggcg                                                                       >  4:83720/1‑149 (MQ=255)                                                                   gTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCAATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGt                                                                           >  3:100765/1‑145 (MQ=255)                                                                   gTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCAATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGt                                                                           <  4:100765/145‑1 (MQ=255)                                                                   gTCGCGAAACGGCAGCTGTTGCTTGAGTGGAGCCGGGAGGCTATCCAGCGCCACGGCCTGTTAACGAACCAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTggc                                                                        <  6:96/148‑1 (MQ=255)                                                                      gcgAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGcc                                                                    >  2:119289/1‑149 (MQ=255)                                                                        gAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTg                                                                  <  3:217641/149‑1 (MQ=255)                                                                          aaCGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGt                                                                <  3:173757/149‑1 (MQ=255)                                                                           aCGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTa                                                               <  5:114495/149‑1 (MQ=255)                                                                            cGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAg                                                              <  4:90221/149‑1 (MQ=255)                                                                            cGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAg                                                              <  7:138959/149‑1 (MQ=255)                                                                             ggCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTTGGGCGGGACGATCAGCGGCGGCGGAAcagcggcagcg                                                                                     >  1:218873/1‑125 (MQ=255)                                                                             ggCAGCTTTTGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTTGGGCGGGACGATCAGCGGCGGCGGAAcagcggcagcg                                                                                     <  2:218873/125‑1 (MQ=255)                                                                               cAGCTTTTGCTTGAGTTGAGCCGTGAGGCAATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGt                                                            >  2:189461/1‑148 (MQ=255)                                                                                   ttttGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCTCTCGGCGGGCt                                                                                                                                  <  3:143902/74‑1 (MQ=37)                                                                                   ttttGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCTCTCGGCGGGCt                                                                                                                                  >  4:143902/1‑74 (MQ=37)                                                                                   ttttGCTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCAcc                                                        >  3:27821/1‑148 (MQ=255)                                                                                        cTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAAcag                                                                                             >  2:71213/1‑106 (MQ=255)                                                                                        cTTGAGTTGAGCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAAcag                                                                                             <  1:71213/106‑1 (MQ=255)                                                                                         ttgagttgagCCGTGAGGCTATCCAGCGCCACGCCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGt                                                 >  7:218439/1‑149 (MQ=255)                                                                                             gttgagCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGGCACCGAGAAGTCCtt                                             <  1:119289/149‑1 (MQ=255)                                                                                                 agCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAg                                          <  2:18626/148‑1 (MQ=255)                                                                                                 agCCGTGAGGCAATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAgg                                         >  5:319823/1‑149 (MQ=255)                                                                                                  gCCGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGc                                        >  2:362406/1‑149 (MQ=255)                                                                                                   ccGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCt                                       >  7:131675/1‑149 (MQ=255)                                                                                                    cGTGAGGCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCtt                                      <  3:83720/149‑1 (MQ=255)                                                                                                         ggCTATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAg                                 >  1:356262/1‑149 (MQ=255)                                                                                                          gctaTCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGCCGGGCTGATCGCTGGGGCGGGACGATCAGCGTCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAgc                                >  3:149806/4‑149 (MQ=255)                                                                                                                cAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCt                          >  6:83680/1‑149 (MQ=255)                                                                                                                   cgcCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTcc                       >  6:200531/1‑149 (MQ=255)                                                                                                                     ccACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCgg                     >  1:191698/1‑149 (MQ=255)                                                                                                                        cGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAg                                 <  2:260960/134‑1 (MQ=255)                                                                                                                               tgAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGt             <  8:218439/147‑1 (MQ=255)                                                                                                                               tgAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTc            <  6:157505/148‑1 (MQ=255)                                                                                                                               tgAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCg           >  7:177634/1‑149 (MQ=255)                                                                                                                               tgAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCg           >  3:169049/1‑149 (MQ=255)                                                                                                                                 aaCGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCgg          <  8:177634/148‑1 (MQ=255)                                                                                                                                 aaCGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCGGc         >  3:242209/1‑149 (MQ=255)                                                                                                                                   cGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCGGCAc       <  2:356262/149‑1 (MQ=255)                                                                                                                                    gAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCGGCACg      <  8:131675/149‑1 (MQ=255)                                                                                                                                     gaaCAGCCCGCCTCTCGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCGGCACg      >  5:271766/2‑148 (MQ=255)                                                                                                                                      aaCAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCGGCACGGa    >  6:240380/1‑149 (MQ=255)                                                                                                                                       aCAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCGGCACGGAt   <  3:140201/149‑1 (MQ=255)                                                                                                                                        cAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCGGCACGGATg  <  1:189461/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                           CAATCGAAAACGCTGAGTGATGGCGACTATTTTGCCCAGGTCGAGCATATCTGCGCAGCGCGTGAGGTCGCGAAACGGCAGCTTTTGCTTGAGTTGAGCCGTGAGGCAATCCAGCGCCACGGCCTGTGAACGAAACAGCCCGCCGAGAGGCGGGCTGTTCGCTGGGGCGGGACGATCAGCGGCGGCGGAACAGCGGCAGCGGTTCGTCGGTGGCGGCCTGGTAGGTCACCGAGAAGTCCTTCAGGCTTTGCAGCGCGTCTTCCGGGTCCTTGTCGGCACGGATG  >  NC_002947/2706903‑2707186

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