Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,484,724	A→T	5.9%	intergenic (+368/‑101)	PP_3976 → / → PP_3977	isochorismatase family hydrolase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,484,724	0	A	T	5.9%	257.4 / 4.9	84	intergenic (+368/‑101)	PP_3976/PP_3977	isochorismatase family hydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/42);  new base T (2/3);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

TCAAACCCTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAA  >  NC_002947/4484580‑4484862                                                                                                                                                   |                                                                                                                                            tCAAACCCTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGATAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGaataat                                                                                                                                         >  2:221932/1‑149 (MQ=255) tCAAACCCTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGattatt                                                                                                                                         >  4:104161/1‑147 (MQ=255)   aaaCCCTGGCAGTTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGATCAGAGATATGCGGTCGTGATCAGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCa                                                                                                                                       <  5:261321/149‑1 (MQ=255)    aaCCCTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAg                                                                                                                                      <  1:319942/149‑1 (MQ=255)    aaCCCTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAg                                                                                                                                      <  4:161373/149‑1 (MQ=255)     tctcTGTCAGGTCGGAAGGGCCAAGCAGATGATGTTGGAGGCATATCGGTCATAGCATGCGGAAATGAACAGAGAAATTCGGGCGTGATACGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCGCTGAATAATCAg                                                                                                                                      <  7:346823/145‑1 (MQ=255)        cTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCt                                                                                                                                  <  8:127186/149‑1 (MQ=255)         tGGCAGGTAGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAAAAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAAAAATCatgagct                                                                                                                                 >  8:346823/1‑143 (MQ=255)           gCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAACA‑CGCAAAAAGCCCTGAATAATCAGGGCttt                                                                                                                                <  1:221932/149‑1 (MQ=255)           gCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCttt                                                                                                                                <  3:141198/148‑1 (MQ=255)            cAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGATTATTCAGGGCttttt                                                                                                                              <  3:104161/149‑3 (MQ=255)            cAGGTCGCAAAGGCCAAGCAGCAGATGTTGGAGGCATATCGGTCACAGCATGCTTAAATGAACCGAGAAATGCGGGCGGGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGATTATTCAGGGCttttt                                                                                                                              <  6:124854/149‑3 (MQ=255)               gTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAata                                                                                                                           <  7:169281/149‑1 (MQ=255)                tCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATAACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAatat                                                                                                                          >  5:354466/1‑149 (MQ=255)                     aCGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCg                                                                                                                      >  3:248420/1‑148 (MQ=255)                              cAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTaga                                                                                                            <  7:166994/149‑1 (MQ=255)                                gatgatGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTa                                                                                                              >  6:107252/1‑145 (MQ=255)                                gatgatGTTGGAGGCATATCGATCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAga‑ga                                                                                                         <  7:181052/149‑1 (MQ=255)                                   gatgTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GAtt                                                                                                       >  5:169124/1‑148 (MQ=255)                                    atgTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTc                                                                                                      >  5:229435/1‑148 (MQ=255)                                    atgTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCa                                                                                                     <  7:39284/149‑1 (MQ=255)                                     tgTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCaa                                                                                                    >  6:260323/1‑149 (MQ=255)                                       ttGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GAtt                                                                                                       >  4:147714/1‑144 (MQ=255)                                       ttGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GAtt                                                                                                       <  3:147714/144‑1 (MQ=255)                                            ggCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTa                                                                                             >  3:238924/1‑149 (MQ=255)                                             gCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTag                                                                                                             <  5:297600/133‑1 (MQ=255)                                             gCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTag                                                                                                             >  6:297600/1‑133 (MQ=255)                                             gCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTAt                                                                                            >  7:55900/1‑149 (MQ=255)                                                 atCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAAc                                                                                                  <  7:260449/139‑1 (MQ=255)                                                 atCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAAc                                                                                                  >  8:260449/1‑139 (MQ=255)                                                 atCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATcc                                                                                          >  1:186706/1‑147 (MQ=255)                                                  tCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTa                                                                                        <  4:248420/148‑1 (MQ=255)                                                      tCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGAGAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCACTGAATAATCAGGGCTTTGAATATGGCGGAATCGTAGA‑GATTCAAACCCTTATCCTAGGCt                                                                                    >  7:55229/1‑148 (MQ=255)                                                          aGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTc                                                                                                      <  7:298718/126‑1 (MQ=255)                                                          aGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTc                                                                                                      >  8:298718/1‑126 (MQ=255)                                                                    aaTGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACt                                                                      >  4:173768/1‑148 (MQ=255)                                                                          cAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTc                                                               <  5:107252/149‑1 (MQ=255)                                                                          cAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAACCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTc                                                               <  2:122710/149‑1 (MQ=255)                                                                           agagAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCt                                                              >  3:53491/1‑149 (MQ=255)                                                                           agagAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCt                                                              <  2:186706/149‑1 (MQ=255)                                                                           agagAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCt                                                              >  4:362223/1‑149 (MQ=255)                                                                           agagAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCACAAAA‑A‑CGAAAAAAGCCATGAATAATCAGGGCTTTGAATATGGCGGAAGCGTATA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCt                                                              >  7:63531/1‑149 (MQ=255)                                                                            gagaAATGCGGGCGTGATCCGTATGACAAGTGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTATCTACAGACTTAATTTCTg                                                             <  8:55900/149‑1 (MQ=255)                                                                             agaAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGTGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGa                                                            >  1:10036/1‑149 (MQ=255)                                                                               aaaTGCGGGCGTGATCCGTATGACAAGGGTAAAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATc                                                          >  1:225515/1‑149 (MQ=255)                                                                                     gggCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTc                                                               >  8:29069/1‑138 (MQ=255)                                                                                     gggCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTc                                                               <  7:29069/138‑1 (MQ=255)                                                                                     gggCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAg                                                     <  4:238924/148‑1 (MQ=255)                                                                                      ggCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCt                                                   <  6:229435/149‑1 (MQ=255)                                                                                       gCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTc                                                  <  5:130994/149‑1 (MQ=255)                                                                                                 tATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCttt                                        <  5:260323/149‑1 (MQ=255)                                                                                                    gACAAGGGTAGAAAATCGGAGCCCCAAAA‑AGCGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTAt                                      >  1:221636/1‑149 (MQ=255)                                                                                                     aCAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGTGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATg                                     <  2:10036/148‑1 (MQ=255)                                                                                                         gggTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTg                                >  1:357075/1‑149 (MQ=255)                                                                                                          ggTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAATAGCCCTGAATAACCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATg                                     >  1:29110/1‑143 (MQ=255)                                                                                                          ggTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAACCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATg                                     <  2:29110/143‑1 (MQ=255)                                                                                                           gTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATc                                                          <  4:61507/121‑1 (MQ=255)                                                                                                           gTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATc                                                          >  3:61507/1‑121 (MQ=255)                                                                                                            tAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCgg                                                        >  6:243189/1‑122 (MQ=255)                                                                                                            tAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCgg                                                        <  5:243189/122‑1 (MQ=255)                                                                                                            tAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGt                               >  5:133610/1‑147 (MQ=255)                                                                                                              gAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGc                           <  3:362223/149‑1 (MQ=255)                                                                                                              gAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGc                           >  2:191446/1‑149 (MQ=255)                                                                                                                aaaTGGGAGACCCAAAA‑T‑CGAAAAAAGCCCTGAATTATGAGGGCTTTGAATATGGCTGAAGCGTAGA‑GATTCAAACCCTTATCCTTGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTGGAGGTCCTTTATGGACTGTAAGCTa                         <  2:225515/149‑1 (MQ=255)                                                                                                                aaaTCTGATCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTa                         <  8:55229/149‑1 (MQ=255)                                                                                                                aaaTCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTa                         <  6:169124/149‑1 (MQ=255)                                                                                                                 aaTCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAAt                                                                  >  7:2869/1‑107 (MQ=255)                                                                                                                 aaTCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAAt                                                                  <  8:2869/107‑1 (MQ=255)                                                                                                                  aTCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGc                       <  1:191446/149‑1 (MQ=255)                                                                                                                   tCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGTATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGc                       >  1:261843/1‑148 (MQ=255)                                                                                                                   tCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGc                       >  6:175407/1‑148 (MQ=255)                                                                                                                    cGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCa                      >  6:227634/1‑148 (MQ=255)                                                                                                                     ggAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCaaa                    <  3:173768/149‑1 (MQ=255)                                                                                                                          cccAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCgg               <  6:354466/149‑1 (MQ=255)                                                                                                                          cccAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCgg               <  4:53491/149‑1 (MQ=255)                                                                                                                           ccAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCggg              >  7:230104/1‑149 (MQ=255)                                                                                                                           ccAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCgg               <  2:6208/148‑1 (MQ=255)                                                                                                                           ccAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCgg               >  1:6208/1‑148 (MQ=255)                                                                                                                                  catattcaaaGCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTc         >  4:342747/8‑149 (MQ=255)                                                                                                                                         aaaGCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGAGGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCgg               >  1:77020/1‑137 (MQ=255)                                                                                                                                         aaaGCCCTGATTATTCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCgg               <  2:77020/136‑1 (MQ=255)                                                                                                                                         aaaGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCaa  <  5:175407/149‑1 (MQ=255)                                                                                                                                             cccTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGc                       >  4:82615/1‑124 (MQ=255)                                                                                                                                             cccTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGc                       <  3:82615/124‑1 (MQ=255)                                                                                                                                               cTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAAc                  <  4:127193/127‑1 (MQ=255)                                                                                                                                               cTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAAc                  >  3:127193/1‑127 (MQ=255)                                                                                                                                                   |                                                                                                                                            TCAAACCCTGGCAGGTCGCAACGGCCAAGCAGATGATGTTGGAGGCATATCGGTCACAGCATGCGGAAATGAACAGAGAAATGCGGGCGTGATCCGTATGACAAGGGTAGAAAATCGGAGCCCCAAAA‑A‑CGCAAAAAGCCCTGAATAATCAGGGCTTTGAATATGGCGGAAGCGTAGA‑GATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAA  >  NC_002947/4484580‑4484862

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