Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,683,744	A→G	9.3%	L274P (CTG→CCG)	PP_2352 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,683,744	0	A	G	9.3%	175.8 / 5.4	64	L274P (CTG→CCG)	PP_2352	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base G (4/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCCAGCGAACACGCCGAGAGGTCATCCAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGA  >  NC_002947/2683597‑2683890                                                                                                                                                    |                                                                                                                                                   gccgccAGCGAACACGCCGAGAGGTCATCCAGCCCATCGATCAGCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCAGAGATCGTCAAAGCCGAGCTCCACCCGCCGGCCATCGATGGCCAg                                                                                                                                                   <  7:234182/149‑1 (MQ=255)   cgccAGCGAACACGCCGAGAGGTCATCCAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAgcg                                                                                                                                                 >  2:143546/1‑149 (MQ=255)       aGCGAACACGCCGAGAGGTCATCCAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCg                                                                                                                                             >  7:248981/1‑149 (MQ=255)                gCCGAGAGGTCATCCAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGgcg                                                                                                                                    >  1:298524/1‑149 (MQ=255)                  cGAGAGGTCATCCAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCt                                                                                                                                  >  1:130014/1‑149 (MQ=255)                   gagaGGTCATCCAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCtt                                                                                                                                 <  5:51999/149‑1 (MQ=255)                   gagaGGTCATCCAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCtt                                                                                                                                 <  3:164600/149‑1 (MQ=255)                       ggTCATCCTGCCCATCGATCCCCCACTCGTCATACTCCCTGGCCAGCAGCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTc                                                                                                                             <  7:165526/149‑1 (MQ=255)                         tCATCCAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCAc                                                                                                                           <  8:19925/149‑1 (MQ=255)                              cAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTc                                                                                                                        >  7:83070/1‑147 (MQ=255)                              cAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCa                                                                                                                      <  6:121232/149‑1 (MQ=255)                              cAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCa                                                                                                                      <  8:248981/149‑1 (MQ=255)                              cAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCa                                                                                                                      <  1:143546/149‑1 (MQ=255)                               aGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAg                                                                                                                     <  3:276070/149‑1 (MQ=255)                                    aTCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCtt                                                                                                                                 >  8:47474/1‑132 (MQ=255)                                    aTCCATCACCCACTCGTCATACTCCCTGGCCACCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGACCAGCGCTCGCAGCGGGCGCtt                                                                                                                                 <  7:47474/132‑1 (MQ=255)                                             acACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCAACACACGCCGGCCAACGATGGCCAGCGCTCGCAGCCGGCGCTTGTTCACGTACAGAATCACCCgtgttt                                                                                                       >  3:154103/2‑145 (MQ=255)                                               aCTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCGGCCCTCGCAGCGGGCGCTTCTTCACGTCCAGTATAACCCGTTgttcg                                                                                                     >  4:107190/1‑149 (MQ=255)                                                    tCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCgg                                                                                                <  2:130014/149‑1 (MQ=255)                                                        aCTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCAGAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACg                                                                                            >  7:291309/1‑149 (MQ=255)                                                        aCTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGgttgtt                                                                                                       >  1:10482/1‑138 (MQ=255)                                                        aCTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGgttgtt                                                                                                       <  2:10482/138‑1 (MQ=255)                                                            ccTGGCCAGCACGAGGTAGTCGATGACCGCAGTGGCCTGCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCGGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGtt                                                                                          <  3:107190/147‑1 (MQ=255)                                                            ccTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTAt                                                                                                                  <  3:158525/123‑1 (MQ=255)                                                            ccTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTAt                                                                                                                  >  4:158525/1‑123 (MQ=255)                                                                gCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGACGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGttcgttcg                                                                                                 <  8:245614/136‑1 (MQ=255)                                                                gCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGAAGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGttcgttcg                                                                                                 >  7:245614/1‑136 (MQ=255)                                                                    gCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTgg                                                                                >  3:108294/1‑149 (MQ=255)                                                                      aCCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGgcg                                                                              >  2:218578/1‑149 (MQ=255)                                                                       ccAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGGCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGgcgc                                                                             <  8:258347/149‑1 (MQ=255)                                                                                       cggaGTGGCCTTCTCGCACAGATCCTCAAAGCCCAGCACCACCCGCCGGCCATCGATGCCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCg                                                             <  8:291309/146‑1 (MQ=255)                                                                                              gCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGa                                                       >  5:106341/1‑148 (MQ=255)                                                                                               ccTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGgt                                                                                                           >  7:35643/1‑95 (MQ=255)                                                                                               ccTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGgt                                                                                                           <  8:35643/95‑1 (MQ=255)                                                                                                  tctcGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATaa                                                  >  2:88698/1‑149 (MQ=255)                                                                                                       cacaGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTgcgc                                                                        >  3:231965/1‑122 (MQ=255)                                                                                                       cacaGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTgcgc                                                                        <  4:231965/122‑1 (MQ=255)                                                                                                        acaGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGccc                                            >  4:273350/1‑149 (MQ=255)                                                                                                            aTCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTggg                                        <  6:106341/149‑1 (MQ=255)                                                                                                                tCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCGGCGCTCGCAGCGGGCGCTGGTTCCCGTCCAGGATCACCGGGTGTTCGTTCGGCACGGTCAGCGCCGGGCGCTGCGCCGCATCACCCGGCCCGACATAATGGCCGGGGGTaa                                     >  8:290059/1‑148 (MQ=255)                                                                                                                tCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGgc                                                                               >  8:134958/1‑106 (MQ=255)                                                                                                                tCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTgg                                         >  5:86037/1‑144 (MQ=255)                                                                                                                tCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTgg                                         <  6:86037/144‑1 (MQ=255)                                                                                                                tCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTaa                                     >  3:181428/1‑148 (MQ=255)                                                                                                                tCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTaa                                     >  6:244627/1‑148 (MQ=255)                                                                                                                 cAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTaaaa                                   >  1:210458/1‑149 (MQ=255)                                                                                                                 cAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTaaaa                                   >  1:248802/1‑149 (MQ=255)                                                                                                                  aaaGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTaaaa                                   <  2:298524/148‑1 (MQ=255)                                                                                                                           caccacCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCt                         >  5:226266/1‑149 (MQ=255)                                                                                                                           caccacCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCAGGGTAAAACGCTGATGCt                         >  4:21543/1‑149 (MQ=255)                                                                                                                              caccCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTggg                                        >  7:51880/1‑131 (MQ=255)                                                                                                                               accCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTcgcg                     <  6:226266/149‑1 (MQ=255)                                                                                                                                 gcgctgGCCATCGATGGCCGGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGt                   >  4:125930/2‑149 (MQ=255)                                                                                                                                 gcgctgGCCATCGATGGCCGGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGt                   >  1:37481/2‑149 (MQ=255)                                                                                                                                  cgccgGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGtt                  >  6:175963/1‑149 (MQ=255)                                                                                                                                   gccgGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTtg                 <  8:83070/149‑1 (MQ=255)                                                                                                                                    ccgGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTtgg                <  4:181428/149‑1 (MQ=255)                                                                                                                                       gCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTtggctg             <  3:273350/149‑1 (MQ=255)                                                                                                                                       gCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTtggctg             <  5:244627/149‑1 (MQ=255)                                                                                                                                              aTGGCCGGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACt      <  2:37481/149‑1 (MQ=255)                                                                                                                                              aTGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACt      >  3:132971/1‑149 (MQ=255)                                                                                                                                                 gCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGc    >  6:27892/1‑148 (MQ=255)                                                                                                                                                  ccAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCAGTTTGCCAGCCAACCGCGa  <  1:88698/149‑1 (MQ=255)                                                                                                                                                   cAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGa  >  3:129904/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCCGCCAGCGAACACGCCGAGAGGTCATCCAGCCCATCGATCACCCACTCGTCATACTCCCTGGCCAGCACCAGGTAGTCGATGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCCTGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGA  >  NC_002947/2683597‑2683890

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