Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,189,577	T→C	6.9%	A148A (GCT→GCC) ‡	PP_1937 →	helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,189,577	0	T	C	6.9%	250.9 / 4.2	85	A148A (GCT→GCC) ‡	PP_1937	helicase
Reads supporting (aligned to +/- strand):  ref base T (37/42);  new base C (3/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 = 1.00e+00

AATGCGAGAGTACTGGCTCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAAC  >  NC_002947/2189430‑2189724                                                                                                                                                    |                                                                                                                                                     aaTGCGAGAGTACTGGCTCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTg                                                                                                                                                     >  8:195765/1‑149 (MQ=255)  aTGCGAGAGTACTGGCTCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTTCAGACGGCTACATGCAAAAATGTACTGCTCCGCTg                                                                                                                                                     >  5:265450/1‑148 (MQ=255)       agagTACTGGCTCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAg                                                                                                                                               <  6:95997/149‑1 (MQ=255)        gagTACTGGCTCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGt                                                                                                                                              <  3:284005/149‑1 (MQ=255)          gTACTGGCTCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCAGCGGAGCAGta                                                                                                                                             <  4:158185/148‑2 (MQ=255)             cTGGCTCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGc                                                                                                                                         <  1:56249/149‑1 (MQ=255)             cTGGCTCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGc                                                                                                                                         >  4:294239/1‑149 (MQ=255)              tGGCTCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCt                                                                                                                                        >  5:69774/1‑149 (MQ=255)                gCTCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCtt                                                                                                                                       <  4:149869/148‑1 (MQ=255)                  tCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTggg                                                                                                                                    >  6:124456/1‑149 (MQ=255)                  tCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTggg                                                                                                                                    >  2:72628/1‑149 (MQ=255)                   cAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAg                                                                                                                                                  >  6:122017/1‑134 (MQ=255)                   cAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAg                                                                                                                                                  <  5:122017/134‑1 (MQ=255)                    aaGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGt                                                                                                                                   <  6:265450/148‑1 (MQ=255)                                tcGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTacac                                                                                                                      >  1:85341/1‑149 (MQ=255)                                tcGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTacac                                                                                                                      >  4:318876/1‑149 (MQ=255)                                  gCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTacac                                                                                                                      <  3:59182/147‑1 (MQ=255)                                  gCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTg                                                                                                                    <  7:195765/149‑1 (MQ=255)                                  gCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTg                                                                                                                    <  6:356360/149‑1 (MQ=255)                                  gCTAGCGCTGTGTGGCCAAGTACTACATTGCAGCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTg                                                                                                                    <  1:72628/149‑1 (MQ=255)                                    tAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGa                                                                                                                  >  4:320715/1‑149 (MQ=255)                                         cTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAgg                                                                                                             <  5:124456/149‑1 (MQ=255)                                          tgtgtgGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACg                                                                                                                                          <  1:140403/119‑1 (MQ=255)                                          tgtgtgGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACg                                                                                                                                          >  2:140403/1‑119 (MQ=255)                                              tgGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑Ta                                                                                                        <  3:294239/148‑1 (MQ=255)                                              tgGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑Ta                                                                                                        <  3:318876/148‑1 (MQ=255)                                                gCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGTTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACa                                                                                                      >  4:263111/1‑148 (MQ=255)                                                        ctacATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCa                                                                                             <  4:330749/149‑1 (MQ=255)                                                        ctacATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCa                                                                                             <  3:320715/149‑1 (MQ=255)                                                             ttGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAg                                                                                        >  2:317431/1‑149 (MQ=255)                                                                  tCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAAc                                                                                   >  7:335040/1‑149 (MQ=255)                                                                    gAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGt                                                                                 >  6:253665/1‑149 (MQ=255)                                                                      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                                                                  <  8:158559/149‑1 (MQ=255)                                                                                   ggTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTgg                                                                  <  6:131826/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                                                           >  4:325808/1‑140 (MQ=255)                                                                                                   gTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAaga                                                           <  3:325808/140‑1 (MQ=255)                                                                                                   gTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCtt                                                   <  1:317431/148‑1 (MQ=255)                                                                                                        ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgtt                                               <  5:219336/147‑1 (MQ=255)                                                                                                        ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgtt                                               >  5:137708/1‑147 (MQ=255)                                                                                                        ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgtt                                               >  2:247629/1‑147 (MQ=255)                                                                                                        ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgtt                                               <  1:247629/147‑1 (MQ=255)                                                                                                        ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgtt                                               >  6:219336/1‑147 (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                  <  4:325988/149‑1 (MQ=255)                                                                                                                                   aaaaTGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTg                  <  3:21750/149‑1 (MQ=255)                                                                                                                                     gcgtcaCTGCTCCAGCGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGGGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTg                  >  5:207606/6‑147 (MQ=255)                                                                                                                                       gtcaCTGCTCCAGCGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACGAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACGTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGc              >  3:234025/4‑149 (MQ=255)                                                                                                                                       gtcaCTGCTCCAGCGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTa                                <  4:356869/128‑1 (MQ=255)                                                                                                                                       gtcaCTGCTCCAGCGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTa                                >  3:356869/4‑131 (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)                                                                                                                                                    |                                                                                                                                                     AATGCGAGAGTACTGGCTCAAGCGCGGGATCTCGCTAGCGCTGTGTGGCCAAGTACTACATTGCATCGAGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTT‑TACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAAC  >  NC_002947/2189430‑2189724

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