Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,263,324	T→A	6.5%	E361D (GAA→GAT)	PP_1105 ←	ATP‑dependent DNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,263,324	0	T	A	6.5%	182.3 / 2.8	61	E361D (GAA→GAT)	PP_1105	ATP‑dependent DNA ligase
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base A (2/2);  total (31/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 = 5.96e-01

GCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATG  >  NC_002947/1263182‑1263471                                                                                                                                                   |                                                                                                                                                    gcCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATcc                                                                                                                                               >  2:299493/1‑149 (MQ=255)   cccACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATccac                                                                                                                                             <  5:386224/149‑1 (MQ=255)   cccACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATccac                                                                                                                                             <  6:144327/149‑1 (MQ=255)    ccACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATccac                                                                                                                                             >  1:361306/1‑148 (MQ=255)    ccACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGAACATccact                                                                                                                                            >  4:261785/1‑148 (MQ=255)     cacaCCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATccacca                                                                                                                                           <  1:58978/149‑1 (MQ=255)     cacaCCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATccacca                                                                                                                                           <  6:335256/149‑1 (MQ=255)     cacaCCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATccacca                                                                                                                                           <  2:211330/149‑1 (MQ=255)     cacaCCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGAACATccactt                                                                                                                                           <  1:275737/149‑3 (MQ=255)      acacCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAg                                                                                                                                          <  2:361306/149‑1 (MQ=255)       cacCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGAACATccacttgc                                                                                                                                         <  3:261785/149‑5 (MQ=255)         ccGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGgtg                                                                                                                                       <  5:243603/149‑1 (MQ=255)                   ccTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTgg                                                                                                                             >  3:318421/1‑149 (MQ=255)                    cTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGc                                                                                                                            >  4:132553/1‑149 (MQ=255)                    cTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGc                                                                                                                            >  8:92961/1‑149 (MQ=255)                      cgcgTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATccacc                                                                                                                                            >  6:323721/1‑131 (MQ=255)                      cgcgTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATccacc                                                                                                                                            <  5:323721/131‑1 (MQ=255)                              ttCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATc                                                                                                                                                <  2:389106/119‑1 (MQ=255)                              ttCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATc                                                                                                                                                >  1:389106/1‑119 (MQ=255)                              ttCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTc                                                                                                                  >  5:172311/1‑149 (MQ=255)                                 aGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGt                                                                                                                <  7:112885/148‑1 (MQ=255)                                 aGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGt                                                                                                                <  8:25099/148‑1 (MQ=255)                                   cATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGGAGAGCCTTGCA‑GCAGAGGCGACAACAGCCCGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGt                                                                                                                >  1:116786/1‑149 (MQ=255)                                     tCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGCGTTCCAGCTGGCTTCGGCGTTCGtgctgc                                                                                                           >  5:62506/1‑149 (MQ=255)                                       aGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGAGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCgg                                                                                                         >  5:125044/1‑149 (MQ=255)                                         cccTTCCACACCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTc                                                                                                       >  4:67175/1‑149 (MQ=255)                                            ttCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGtg                                                                                                               >  4:356913/1‑138 (MQ=255)                                            ttCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGtg                                                                                                               <  3:356913/138‑1 (MQ=255)                                                   cccAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCa                                                                                             <  3:113363/149‑1 (MQ=255)                                                          tGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCAGGAAGAGGCGACAACAGCACGGGGGCAATTGGATTTTCAACCACCAGGTGTACAAGATGGATTCGCCGTACGTCATGCGTACGGATCCGCCATTcaac                                                                                       >  1:291055/1‑146 (MQ=255)                                                            cGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCg                                                                                      <  6:302522/147‑1 (MQ=255)                                                               gACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTg                                                                                 >  7:124145/1‑149 (MQ=255)                                                                 cTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGTCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGc                                                                                <  8:120351/148‑1 (MQ=255)                                                                    ttCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGTTGCGCCAGTCCTCGCCCTGCCAc                                                                             <  8:115516/148‑1 (MQ=255)                                                                       gcgcTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTcc                                                                          <  4:318421/148‑1 (MQ=255)                                                                         gcTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCcagc                                                                       <  1:299493/149‑1 (MQ=255)                                                                         gcTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCcagc                                                                       <  3:132553/149‑1 (MQ=255)                                                                                   aaGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAgg                                                             >  7:143168/1‑149 (MQ=255)                                                                                   aaGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAgg                                                             >  4:190487/1‑149 (MQ=255)                                                                                      gTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCt                                                          >  5:170790/1‑149 (MQ=255)                                                                                       tCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTA‑GCCTgc                                                        <  3:67175/149‑1 (MQ=255)                                                                                          gCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTgcag                                                      >  2:315890/1‑149 (MQ=255)                                                                                             cAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTGCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGcac                                                   >  8:111227/1‑149 (MQ=255)                                                                                              aGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGcac                                                   >  2:112922/1‑148 (MQ=255)                                                                                                    g‑‑‑gCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGgc                                            <  2:291055/149‑1 (MQ=255)                                                                                                             gca‑gcaGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTgca                                                       <  5:227326/132‑1 (MQ=255)                                                                                                             gca‑gcaGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTgca                                                       >  6:227326/1‑132 (MQ=255)                                                                                                             gca‑gcaGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCg                                      >  8:218614/1‑149 (MQ=255)                                                                                                             gca‑gcaGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCg                                      >  3:99141/1‑149 (MQ=255)                                                                                                               a‑gcagAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGtt                                    <  8:143168/149‑1 (MQ=255)                                                                                                               a‑gcagAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGtt                                    >  6:76399/1‑149 (MQ=255)                                                                                                                 gcagAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGtt                                    >  6:21002/1‑148 (MQ=255)                                                                                                                        cGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGtt                            <  1:112922/149‑1 (MQ=255)                                                                                                                             acaGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTg                       <  8:290430/149‑1 (MQ=255)                                                                                                                             acaGCACGGGGGCAAGTGGATGTTCATCCACCAGGCGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTg                       <  6:62506/149‑1 (MQ=255)                                                                                                                                    gggggCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCg                >  4:229608/1‑149 (MQ=255)                                                                                                                                       ggCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGttt             <  4:99141/149‑1 (MQ=255)                                                                                                                                          tggTGGATGAACATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGc          >  2:363263/3‑149 (MQ=255)                                                                                                                                            gTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGa        <  7:92961/149‑1 (MQ=255)                                                                                                                                             tGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGAc       >  5:258551/1‑149 (MQ=255)                                                                                                                                             tGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGAc       >  5:252999/1‑149 (MQ=255)                                                                                                                                                   ttCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATg  <  6:145089/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GCCCCACACCGTACAGGGCCTCGCGTTGCTTCAGCATCAGCCCTTCCACGCCCAGGCTGCGTGACTGTTCGCGCTGGCGGGCAAGGTCGGCCCAGTCCGG‑‑‑GCCTTGCA‑GCAGAGGCGACAACAGCACGGGGGCAAGTGGATGTTCATCCACCAGGTGTTCCAGCTGGCTTCGGCGTTCGTGCTGCGGTCGGCTGCGCCAGTCCTCGCCCTGCCACTCCAGCAGGTCGTAGGCCTGCAGCACCACCGGCGCATCGTTCAGCAGTTTTTTGCCCAGCGTTTTGCGACCGATG  >  NC_002947/1263182‑1263471

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