Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,519,157	A→G	6.4%	K136K (AAA→AAG)	PP_2212 →	leucine‑rich repeat‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,519,157	0	A	G	6.4%	179.4 / 2.1	63	K136K (AAA→AAG)	PP_2212	leucine‑rich repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base G (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

GCTTCGGCACACTGCCAGGTACCTTGCAACAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCGGGGTTGAAGCTGGCGA  >  NC_002947/2519012‑2519305                                                                                                                                                  |                                                                                                                                                     gCTTCGGCACACTGCCAGGTACCTTGCAACAGACATACGGGTACAGCCAGCCGCGTGAGAAAGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGCTTTTGTCAGGCATCGATAAAcga                                                                                                                                                   >  6:351367/1‑147 (MQ=255)   ttCGGCACACTGCCAGGTACCTTGCAACAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGACTACCCGCTACTACGAGGGCACGGGGCTGGGTTTGTCAGGCAGCGATAAACCAtt                                                                                                                                                 >  8:82591/1‑149 (MQ=255)      ggCACACTGCCAGGTACCTTGCAACAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGaa                                                                                                                                              <  5:238464/149‑1 (MQ=255)          cacTGCCAGGTACCTTGCAACAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGAGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGc                                                                                                                                           >  5:73666/1‑148 (MQ=255)            cTGCCAGGTACCTTGCAACAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATc                                                                                                                                        <  3:6457/149‑1 (MQ=255)               ccAGGTACCTTGCAACAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCt                                                                                                                                     <  5:263844/149‑1 (MQ=255)                  ggTACCTTGCAACAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGa                                                                                                                                  <  2:108190/149‑1 (MQ=255)                       cTTGCAACAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGGTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCAt                                                                                                                             >  8:136108/1‑149 (MQ=255)                        ttGCAACAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCAtt                                                                                                                            >  2:277487/1‑149 (MQ=255)                            aaCAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCAtt                                                                                                                                                 >  3:117738/1‑124 (MQ=255)                            aaCAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCAtt                                                                                                                                                 <  4:117738/124‑1 (MQ=255)                            aaCAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGAt                                                                                                                         >  3:144966/1‑148 (MQ=255)                             aCAGACATACGGGTACAGCCAGACGCGAGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTa                                                                                                                       <  3:294256/149‑1 (MQ=255)                               aGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAGCCATTGAATGCATCGCTGGAAGCATTGATTatg                                                                                                                     >  7:58914/1‑147 (MQ=255)                               aGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAGCCATTGAATGCATCGCTGGAAGCATTGATTatg                                                                                                                     >  7:64123/1‑147 (MQ=255)                               aGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACg                                                                                                                     <  2:127602/149‑1 (MQ=255)                                     aCGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGt                                                                                                               >  1:306002/1‑149 (MQ=255)                                               ccAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTgg                                                                                                     <  2:188916/149‑1 (MQ=255)                                               ccAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTgg                                                                                                     <  3:96679/149‑1 (MQ=255)                                                 aGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTg                                                                                                                                                <  6:165264/104‑1 (MQ=255)                                                 aGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTg                                                                                                                                                >  5:165264/1‑104 (MQ=255)                                                     gcgTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTg                                                                                               <  8:39619/149‑1 (MQ=255)                                                           gaACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGc                                                                                         >  2:194159/1‑149 (MQ=255)                                                                 gTTACTGCGGCTTATGCGGAACATTGCTGCGAAAACCGGCGACTAGGAGGCCACGGGGCAGGATTTGTCAGGCAGGGAAACACCATTGCATGCATCGCTGGAAGCATTGATTACGGCCAGACGCAATCTCGACGTTGGTCCGCGTTAt                                                                                    >  4:227710/1‑148 (MQ=255)                                                                    aCTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAg                                                                                 <  7:136108/148‑1 (MQ=255)                                                                      tGCGGCTTATCCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGCCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTTGTCAGCGTTATCAGGAc                                                                              >  3:71206/1‑149 (MQ=255)                                                                          gCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGTGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATCCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGt                                                                          <  5:351367/149‑1 (MQ=255)                                                                          gCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGt                                                                          >  4:89698/1‑149 (MQ=255)                                                                              aTGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTg                                                                        <  4:82861/147‑1 (MQ=255)                                                                                 cgcAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAg                                                                   <  1:194159/149‑1 (MQ=255)                                                                                   cAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGa                                                                 >  7:159411/1‑149 (MQ=255)                                                                                      cTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCt                                                              >  4:164018/1‑149 (MQ=255)                                                                                       tttGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGt                                                                          >  5:106705/1‑136 (MQ=255)                                                                                       tttGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGt                                                                          <  6:106705/136‑1 (MQ=255)                                                                                         tgctgcGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTc                                                            <  6:15899/148‑1 (MQ=255)                                                                                         tgctgcGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTc                                                            <  6:73666/148‑1 (MQ=255)                                                                                              cGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGc                                                        <  8:159411/147‑1 (MQ=255)                                                                                                aaTACCCGCTAATACTATGGCACGGGGCTGGTTTTGTCAGGCATCTATAAACCATTGAATGCATCGCTGGAAGCATTGATTAAGGCCTGTCGCAATCTGGAAGTTGGTCATCGTTATCAGGACGCGTTGAAAAGTATCTTCCCTCCAg                                                     <  2:306002/148‑1 (MQ=255)                                                                                                     ccGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGAcgc                                                                            <  8:35030/120‑1 (MQ=255)                                                                                                     ccGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGAcgc                                                                            >  7:35030/1‑120 (MQ=255)                                                                                                     ccGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTAcc                                                >  3:45333/1‑148 (MQ=255)                                                                                                     ccGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACAATTGAATGCATAGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCACGCCAGCTAcc                                                >  6:377395/1‑148 (MQ=255)                                                                                                     ccGCTACTACGAGGGCACGGGGCTGGGTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGAGGTTGGTCGGCGTTATCAGGACGCGTGGGAATGGATCTTCCCGCCAGCTAcc                                                >  8:177634/1‑148 (MQ=255)                                                                                                        ctactaCGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAgg                                            <  1:277487/149‑1 (MQ=255)                                                                                                              cGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTg                                       <  4:144966/148‑1 (MQ=255)                                                                                                                 gggCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTgctgg                                   >  1:254365/1‑149 (MQ=255)                                                                                                                 gggCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTgctgg                                   >  5:32369/1‑149 (MQ=255)                                                                                                                   gCACGCGGCGGGTTTTGTCAGGCCGCGATAAGCCATTGAATGCATCGCTGGAATCATTGATTATGGCCTGTCGCAATCTGGACGTTGGTCAGCTTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCATGCGGTgctggct                                 <  8:58914/149‑1 (MQ=255)                                                                                                                   ccACGGGGCTGGTTTTCGCAGGCAGCGATAAGCCATTGAATGCATCGCTGGAAGCACTGATAAAGGCCTGTCGCAACCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTgctggct                                 <  8:64123/148‑1 (MQ=255)                                                                                                                     aCGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTgctggctg                                >  5:322038/1‑148 (MQ=255)                                                                                                                      cGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACTGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGaa                              <  3:89698/149‑1 (MQ=255)                                                                                                                                   gTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGaa                              <  7:15321/136‑1 (MQ=255)                                                                                                                                   gTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGaa                              >  8:15321/1‑136 (MQ=255)                                                                                                                                   gTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCgg                 >  1:85586/1‑149 (MQ=255)                                                                                                                                    tCAGGCAGCGATAACCGATTGAATGCATCGCTGGAAGAATTGATTACGGCCTGGCGCATGCTGGAGGTTGGTCAGCGGTATCAGGTCGCGTTGCAAAGGATGTTCCCGCCAGGTACGCCGGCGGTGCTGGCTGAAGACAAACGTTCgg                 <  3:227710/148‑1 (MQ=255)                                                                                                                                    tCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCggg                >  8:148531/1‑149 (MQ=255)                                                                                                                                      aGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCGGGGt              <  7:82591/149‑1 (MQ=255)                                                                                                                                        gCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCGGGGtt             >  1:36995/1‑148 (MQ=255)                                                                                                                                        gCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCGGGGtt             >  8:65687/1‑148 (MQ=255)                                                                                                                                          aGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTcc                                                           >  5:104270/1‑100 (MQ=255)                                                                                                                                          aGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTcc                                                           <  6:104270/100‑1 (MQ=255)                                                                                                                                            cGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCGGGGTTGAAGc        >  4:224000/1‑149 (MQ=255)                                                                                                                                              aTAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCGGGGTTGAAGCTg      >  1:354410/1‑149 (MQ=255)                                                                                                                                                 aaCCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCGGGGTTGAAGCTGGc    <  2:254365/148‑1 (MQ=255)                                                                                                                                                 aaCCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCGGGGTTGAAGCTGGc    <  6:32369/148‑1 (MQ=255)                                                                                                                                                  aCCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCGGGGTTGAAGCTGGCGa  <  2:354410/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GCTTCGGCACACTGCCAGGTACCTTGCAACAGACATACGGGTACAGCCAGACGCGTGAGAACGGGTTACTGCGGCTTATGCGCAACTTTGCTGCGAATACCCGCTACTACGAGGGCACGGGGCTGGTTTTGTCAGGCAGCGATAAACCATTGAATGCATCGCTGGAAGCATTGATTACGGCCTGTCGCAATCTGGACGTTGGTCAGCGTTATCAGGACGCGTTGCAAAGGATCTTCCCGCCAGCTACCCAGGCGGTGCTGGCTGAAGACAAACGTTCGGGGTTGAAGCTGGCGA  >  NC_002947/2519012‑2519305

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