Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	792,038	G→A	5.7%	F338F (TTC→TTT)	PP_0678 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	792,038	0	G	A	5.7%	209.3 / 2.3	70	F338F (TTC→TTT)	PP_0678	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/30);  new base A (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GCGGGTTGCTGTTGATCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGAGAGGGGCGGTAAGTTTCGCTGGCGACCAGCCCGGACTGCTTC  >  NC_002947/791896‑792185                                                                                                                                               |                                                                                                                                                    gcggGTTGCTGTTGATCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTg                                                                                                                                              <  1:1430261/150‑1 (MQ=255) gcggGTTGCTGTTGATCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTg                                                                                                                                              <  1:520300/150‑1 (MQ=255)   gggTTGCTGTTGATCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTggg                                                                                                                                            <  1:865053/150‑1 (MQ=255)   gggTTGCTGTTGATCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTggg                                                                                                                                            <  2:1624418/150‑1 (MQ=255)       tGCTGTTGATCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCttt                                                                                                                                        <  1:1363497/150‑1 (MQ=255)           gtgggtCGCCGGGAGCAGCTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCaa                                                                                                                                    <  2:100327/145‑1 (MQ=255)           gTTGATCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCaa                                                                                                                                    >  2:161226/1‑150 (MQ=255)             tGATCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAAtt                                                                                                                                  <  1:1429505/150‑1 (MQ=255)               aTCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGTAATTgc                                                                                                                                <  1:1712827/150‑1 (MQ=255)               aTCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTgc                                                                                                                                <  2:1007555/150‑1 (MQ=255)               aTCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTgc                                                                                                                                >  2:317147/1‑150 (MQ=255)                    ccGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAgt                                                                                                                           >  1:445908/1‑150 (MQ=255)                       cagcaACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTc                                                                                                                        <  2:545119/150‑1 (MQ=255)                              tCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCtt                                                                                                                 >  2:1697369/1‑150 (MQ=255)                                gAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATAGCGCTCGAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCGTAc                                                                                                               >  2:544910/1‑150 (MQ=255)                                  aGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTAcag                                                                                                             >  1:316000/1‑150 (MQ=255)                                   gCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTAcagc                                                                                                            <  1:1283195/150‑1 (MQ=255)                                      gCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTAcagccag                                                                                                         <  1:317147/150‑1 (MQ=255)                                             gCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCCAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCATGGCCTGGAATTCTGGGCTTTGCACTTGCGTAGTGTCGGCGATTACTGCCAGCACCTgg                                                                                                  >  2:852922/1‑150 (MQ=255)                                             gCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTgg                                                                                                  >  2:366243/1‑150 (MQ=255)                                             gCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTgg                                                                                                  >  1:477058/1‑150 (MQ=255)                                              cTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCANGGCGNTACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGgc                                                                                                 >  1:208045/1‑150 (MQ=255)                                                  cccTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGATGGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCtt                                                                                             >  2:153721/1‑150 (MQ=255)                                                   ccTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCttt                                                                                            >  2:336614/1‑150 (MQ=255)                                                    cTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCCATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCtttt                                                                                           >  2:1342933/1‑150 (MQ=255)                                                     tGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTtc                                                                                          <  1:161226/150‑1 (MQ=255)                                                        gTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTtcgtc                                                                                       <  1:1697369/150‑1 (MQ=255)                                                         tCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTgg                                                                                                                                             <  2:1547043/95‑1 (MQ=255)                                                         tCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTgg                                                                                                                                             >  1:1547043/1‑95 (MQ=255)                                                           gcgTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTg                                                                                    <  2:1845797/150‑1 (MQ=255)                                                            ggTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTgg                                                                                   <  2:154503/149‑1 (MQ=255)                                                             gTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTgc                                                                                                                                >  1:1666575/1‑104 (MQ=255)                                                             gTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTgc                                                                                                                                <  2:1666575/104‑1 (MQ=255)                                                                 cTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCtttt                                                                                           <  1:1428335/137‑1 (MQ=255)                                                                 cTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCtttt                                                                                           >  2:1428335/1‑137 (MQ=255)                                                                 cTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAg                                                                              >  2:9010/1‑150 (MQ=255)                                                                 cTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAg                                                                              >  2:1077730/1‑150 (MQ=255)                                                                  tGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGt                                                                             >  2:1868645/1‑150 (MQ=255)                                                                    gCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTgc                                                                                                                                <  2:1071082/97‑1 (MQ=255)                                                                    gCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTgc                                                                                                                                >  1:1071082/1‑97 (MQ=255)                                                                             ccGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCttt                                                                                            >  2:1407617/1‑124 (MQ=255)                                                                             ccGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCttt                                                                                            <  1:1407617/124‑1 (MQ=255)                                                                                aTACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGc                                                               >  2:1044808/1‑150 (MQ=255)                                                                                    cGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCt                                                                  >  1:456793/1‑143 (MQ=255)                                                                                    cGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCt                                                                  <  2:456793/143‑1 (MQ=255)                                                                                      cTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGAAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTgg                                                                                                  <  2:1460959/109‑1 (MQ=255)                                                                                      cTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGAAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTgg                                                                                                  >  1:1460959/1‑109 (MQ=255)                                                                                      cTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGAAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTgg                                                                                                  >  1:1461134/1‑109 (MQ=255)                                                                                      cTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGAAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTgg                                                                                                  <  2:1461134/109‑1 (MQ=255)                                                                                       tCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTgg                                                                                                                                             <  2:1497380/65‑1 (MQ=255)                                                                                       tCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTgg                                                                                                                                             >  1:1497380/1‑65 (MQ=255)                                                                                          aGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCa                                                     <  2:445908/150‑1 (MQ=255)                                                                                            gCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACg                                                   <  1:366243/150‑1 (MQ=255)                                                                                              gATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGtt                                                 >  2:1861593/1‑150 (MQ=255)                                                                                                   gcgGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTc                                                                                 >  2:1495850/1‑113 (MQ=255)                                                                                                   gcgGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTc                                                                                 <  1:1495850/113‑1 (MQ=255)                                                                                                    cgGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCgag                                           <  1:9010/150‑1 (MQ=255)                                                                                                     ggCATCTTGCACCCACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCgaga                                          <  1:25187/150‑1 (MQ=255)                                                                                                      gCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGagag                                         >  2:1836546/1‑150 (MQ=255)                                                                                                         tCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGAGAgggg                                      <  1:734531/150‑1 (MQ=255)                                                                                                               aCCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGAGAGGGGCGGTaa                                >  1:1739299/1‑150 (MQ=255)                                                                                                                ccAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGAGAGGGGCGGTAAg                               >  2:686800/1‑150 (MQ=255)                                                                                                                 cAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGAGAGGGGCGGTAAGt                              >  1:1882840/1‑150 (MQ=255)                                                                                                                            gtcccccaGGCAGGGCCTGGAATTCTGGGATTTGCTATTGTGCAGTGTCGGGGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGAGAGGGGCGGTAAGTTTCGCTGGCGa                   <  2:312082/143‑1 (MQ=255)                                                                                                                             ttCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAAGTCACGTTCGCGAGAGGGGCGGTAAGTTTCGCTGGCGAc                  >  1:312082/1‑150 (MQ=255)                                                                                                                                gcgAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGAGAGGGGCGGTGAGGTTCGCTGGCGATCAg               >  2:1535137/1‑150 (MQ=255)                                                                                                                                gcgAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGAGAGGGGCGGTAAGTTTCGCTGGCGACCAg               >  1:743431/1‑150 (MQ=255)                                                                                                                                       aGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGAGAGGGGCGGTAAGTTTCGCTGGCGACCAGCCCGGAc        >  2:1806252/1‑150 (MQ=255)                                                                                                                                             tGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAACGATTCACGTTCGCGAGACGGGCGGTAAGTTTCGCTGGCGACCCGCCCGGACTGCTTc  >  1:1756035/1‑150 (MQ=255)                                                                                                                                             tGGAATTCTGGGCGTCGCAATTGCTCAGGGTCGGCGCTTACAGCCAGCACCTTGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGAGAGGGGCGGTACGTTGCGCTGTCGACCAGCCCGGACTGCttt  >  2:77451/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                    GCGGGTTGCTGTTGATCGCCCGCAGCAACTCGAAGCTGCGTTTGGCTTGCCCTGGGTCGCGTTCCTGGCCGTTGGCCCGATACCGCTCAAGGCGATACGCGGCATCTTGCACCAACGCTTGGCGTTCGCGAGGCAGGGCCTGGAATTCTGGGCTTTGCAATTGCGCAGTGTCGGCGCTTACAGCCAGCACCTGGCGCTTTTCGTCGTGGTCGAGTGGCTCGGCCCTGGCCAGCAATTCACGTTCGCGAGAGGGGCGGTAAGTTTCGCTGGCGACCAGCCCGGACTGCTTC  >  NC_002947/791896‑792185

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑