Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,189,646	G→T	5.4%	V171V (GTG→GTT)	PP_1937 →	helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,189,646	0	G	T	5.4%	274.0 / 3.3	85	V171V (GTG→GTT)	PP_1937	helicase
Reads supporting (aligned to +/- strand):  ref base G (38/42);  new base T (2/3);  total (40/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 = 9.99e-01

GCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTT  >  NC_002947/2189499‑2189789                                                                                                                                                     |                                                                                                                                                gCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTgc                                                                                                                                                <  6:69774/149‑1 (MQ=255)  cTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGGACGTgcg                                                                                                                                               >  4:331602/1‑149 (MQ=255)      ttCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCtt                                                                                                                                            >  5:372872/1‑148 (MQ=255)      ttCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACAGCTCCCCTGGAGCAGTGACGCTTTGGTCCAGCAACTACACTGGACCAGGTT‑TACAAACTCAGCATGAATCGAACGTGCGCtt                                                                                                                                            >  1:278986/1‑148 (MQ=255)         gcgcCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTATACATACTCAGCATGAAGCGCACGTTCGCTTCat                                                                                                                                         <  1:50054/149‑2 (MQ=255)         gcgcCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGCACTTTCGCTTCatg                                                                                                                                        <  8:17957/149‑3 (MQ=255)         gcgcCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGCACGTTCGCTTCat                                                                                                                                         >  7:17957/1‑147 (MQ=255)         gcgcCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCAcc                                                                                                                                        >  3:325988/1‑149 (MQ=255)              ggTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTgg                                                                                                                                   <  6:131826/149‑1 (MQ=255)              ggTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTgg                                                                                                                                   <  8:158559/149‑1 (MQ=255)                  aCGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGa                                                                                                                                >  5:274760/1‑148 (MQ=255)                    gTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACTTGCGCTTCACCTCTGGGGa                                                                                                                                >  2:290855/1‑146 (MQ=255)                     tGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAaga                                                                                                                            <  8:335040/149‑1 (MQ=255)                       cGAATTCGTTACCCGACAGGGTACAGACGGTTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAgc                                                                                                                          <  3:263111/149‑1 (MQ=255)                       cGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAagag                                                                                                                           <  7:127709/148‑1 (MQ=255)                        gAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAgcg                                                                                                                         >  4:21750/1‑149 (MQ=255)                         aaTTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCAc                                                                                                                                         <  8:291200/132‑1 (MQ=255)                         aaTTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCAc                                                                                                                                         >  7:291200/1‑132 (MQ=255)                              gTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACTCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCttt                                                                                                                   <  8:300105/149‑1 (MQ=255)                              gTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAaga                                                                                                                            <  3:325808/140‑1 (MQ=255)                              gTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAaga                                                                                                                            >  4:325808/1‑140 (MQ=255)                              gTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCtt                                                                                                                    <  1:317431/148‑1 (MQ=255)                                   ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgtt                                                                                                                >  6:219336/1‑147 (MQ=255)                                   ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgtt                                                                                                                >  2:247629/1‑147 (MQ=255)                                   ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgtt                                                                                                                <  5:219336/147‑1 (MQ=255)                                   ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgtt                                                                                                                >  5:137708/1‑147 (MQ=255)                                   ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgtt                                                                                                                <  1:247629/147‑1 (MQ=255)                                       cAGGGTACAGACTGCTACATGCAAAAATGTACTGCTCCGCTTGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGAcc                                                                                                          <  2:278986/149‑1 (MQ=255)                                       cAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGAcc                                                                                                          <  6:372872/149‑1 (MQ=255)                                          ggTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAAt                                                                                                       <  6:274760/149‑1 (MQ=255)                                           gTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCT‑GGGAAAGAGCGCGCTTTGTTGAGACCAAtct                                                                                                     >  4:127307/1‑149 (MQ=255)                                                 aCGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAgc                                                                                                                          <  4:144532/123‑1 (MQ=255)                                                         aTGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTa                                                                                                 >  7:348573/1‑140 (MQ=255)                                                         aTGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTa                                                                                                 <  8:348573/140‑1 (MQ=255)                                                          tGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCt                                                                                       >  1:236470/1‑149 (MQ=255)                                                             aaaaaTGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCtttt                                                                                    >  3:155542/1‑149 (MQ=255)                                                              aataTGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTg                                                                                   <  1:290855/146‑1 (MQ=255)                                                              aaaaTGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTg                                                                                   <  3:21750/149‑1 (MQ=255)                                                              aaaaTGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTg                                                                                   <  4:325988/149‑1 (MQ=255)                                                                gcgtcaCTGCTCCAGCGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGGGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTg                                                                                   >  5:207606/6‑147 (MQ=255)                                                                  gtcaCTGCTCCAGCGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACGAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACGTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGc                                                                               >  3:234025/4‑149 (MQ=255)                                                                  gtcaCTGCTCCAGCGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTa                                                                                                 >  3:356869/4‑131 (MQ=255)                                                                  gtcaCTGCTCCAGCGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTa                                                                                                 <  4:356869/128‑1 (MQ=255)                                                                   gTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCt                                                                              >  7:10417/1‑149 (MQ=255)                                                                     aCTGCTCCAGCGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTg                                                                            <  4:234025/149‑1 (MQ=255)                                                                         cTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGGACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTggg                                                                          <  3:331602/147‑1 (MQ=255)                                                                         cTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGa                                                                                            >  2:3517/1‑129 (MQ=255)                                                                         cTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGa                                                                                            <  1:3517/129‑1 (MQ=255)                                                                         cTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTggg                                                                          <  4:157958/147‑1 (MQ=255)                                                                         cTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTggg                                                                          >  2:79703/1‑147 (MQ=255)                                                                          tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTggg                                                                          >  1:221584/1‑146 (MQ=255)                                                                          tCCGCTGGAGCAGTCACGCTTGTGTAAAGCAACTACACTGGACCACGTT‑TACATACTCAGCATGAAGAGAACGTGCGCTACACCTCTGGGGAAAGACCGCGCTTTGTTGAGACCAATCTCCTAGATCAGTGCTTTTGGATCTTggg                                                                          >  8:279761/1‑146 (MQ=255)                                                                              cTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGATCGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAAc                                                                   <  8:10417/149‑1 (MQ=255)                                                                                 gAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGcc                                                                <  5:253665/149‑1 (MQ=255)                                                                                      gTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGa                                                           >  3:362344/1‑149 (MQ=255)                                                                                         aCGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGc                                                        >  4:8243/1‑149 (MQ=255)                                                                                         aCGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGc                                                        >  4:205644/1‑149 (MQ=255)                                                                                          cGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttg                                                                                                                  >  6:60183/1‑90 (MQ=255)                                                                                          cGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttg                                                                                                                  <  5:60183/90‑1 (MQ=255)                                                                                           gCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCtt                                                      >  8:89483/1‑149 (MQ=255)                                                                                                ggTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCAtt                                                 <  2:85341/149‑1 (MQ=255)                                                                                                ggTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCAtt                                                 <  7:89483/149‑1 (MQ=255)                                                                                                  tCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGa                                               <  1:79703/149‑1 (MQ=255)                                                                                                   ccAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGa                                                              >  2:61981/1‑133 (MQ=255)                                                                                                   ccAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGa                                                              <  1:61981/133‑1 (MQ=255)                                                                                                                  gACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACg                                                                  >  4:197811/1‑114 (MQ=255)                                                                                                                  gACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACg                                                                  <  3:197811/114‑1 (MQ=255)                                                                                                                   aCCATGTT‑TACATTCTCAGGATGACGCGAACGTGCTCTTCACCTCGGGGTAAAGAGAGGGCTTTGTGGAGACGAATCTGCTCGGTGAGTGCTTTTGGAGCTTGGGGAGTAAGTCCGATGACGCTTTCATTTTGCTGCTAGAGCAGCTtc                              <  4:356921/149‑1 (MQ=255)                                                                                                                   aCCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTCTAGCAGCtca                              <  6:137708/149‑3 (MQ=255)                                                                                                                     cAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTtctc                            <  2:236470/149‑1 (MQ=255)                                                                                                                             aCATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTc                      <  3:8243/148‑1 (MQ=255)                                                                                                                              cATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGt                    >  5:354827/1‑149 (MQ=255)                                                                                                                                 aCTCAGCATGAAGCGAACGTTCGCTTCTCCTCTGTGTAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTTTAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGATCTGCTa                                        <  2:321652/126‑1 (MQ=255)                                                                                                                                 aCTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTa                                        >  1:321652/1‑126 (MQ=255)                                                                                                                                     gaggtGAAGCGCACGTTCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGCCGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTgg             >  8:16510/5‑149 (MQ=255)                                                                                                                                     aGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCaa                          >  1:72854/1‑136 (MQ=255)                                                                                                                                     aGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCaa                          <  2:72854/136‑1 (MQ=255)                                                                                                                                     aGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTgg             >  7:289825/1‑149 (MQ=255)                                                                                                                                     aGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTgg             <  3:205644/149‑1 (MQ=255)                                                                                                                                     aGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTg              <  6:207606/148‑1 (MQ=255)                                                                                                                                     aGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTAACCTgg             >  7:284589/1‑149 (MQ=255)                                                                                                                                      gCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTgg             >  2:272775/1‑148 (MQ=255)                                                                                                                                         tGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGc            >  3:354544/1‑146 (MQ=255)                                                                                                                                          gAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCaa          >  1:280151/1‑147 (MQ=255)                                                                                                                                           aaGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAaga        <  4:354544/148‑1 (MQ=255)                                                                                                                                               gAACGTGCGCTTCACCTCT‑GGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCttt  <  3:127307/149‑1 (MQ=255)                                                                                                                                                aaCGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCttt  <  7:279761/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                GCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTT  >  NC_002947/2189499‑2189789

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