Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,627,312	C→T	6.9%	V64I (GTC→ATC)	PP_4944 ←	carbamoyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,627,312	0	C	T	6.9%	171.2 / 2.6	59	V64I (GTC→ATC)	PP_4944	carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/28);  new base T (2/2);  total (29/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 = 9.12e-01

CTTGCGATAGCGCTTGTAGCGACGATTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGG  >  NC_002947/5627169‑5627457                                                                                                                                                |                                                                                                                                                  cttGCGATAGCGCTTGTAGCGACGATTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCacgtcgacgtcg                                                                                                                                               <  2:6487/148‑1 (MQ=255) cttGCGATAGCGCTTGTAGCGACGATTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGAGCGAACGGGATGGCTACCacgtcgacgtc                                                                                                                                                <  4:141370/147‑1 (MQ=255)    gCGATAGCGCTTGTAGCGACGATTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACAACGTCGATGTAGGACg                                                                                                                                           >  4:199052/1‑149 (MQ=255)    gCGATAGCGCTTGTAGCGACGATTGCCCATTAGGATGGAATCGAGTGCACGAACCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGCAGAGGCTCATCTGGGCGAACGGGATGGATACAACGTCGATGTCGGACg                                                                                                                                           >  4:203611/1‑149 (MQ=255)       aTAGCGCTTGTAGCGACGATTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCt                                                                                                                                        <  5:202747/149‑1 (MQ=255)                  aGCGACGATTGCCCATTAGGATGCCATCGAGTGCAAGAACCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGCAGTCGGACGGCTTGAtgcctgcc                                                                                                                             >  8:300842/1‑149 (MQ=255)                  aGCGAAGATTGCCCATTACGATGGCAACGAGTGCACGATCCGGGGCGTACCAATAGCGCTTCGCATAGTGCCACCGCGCCTTGCCGAACAGGCTCATCGGGGCGAACGGGCTGGCGACGACGTCGCCGTCGGACGGCTTGCtgcctgcc                                                                                                                             >  8:295959/1‑149 (MQ=255)                   gCGACGTTTGCCCATTAGGATGGCATCGAGGGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACCGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGTTGTCGGACGGCTTGATgcctgcct                                                                                                                            <  3:203611/149‑1 (MQ=255)                   gCGACGATTGCGCTTTAGGATTGCATAGTGTGCAGGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGGTGATCGGGGCGAACGGGATGGGTACCACGTCGATGTCGGACGGCTTGATgcctgcct                                                                                                                            <  3:199052/149‑1 (MQ=255)                   gCGACGATTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATgcctgcct                                                                                                                            <  6:97495/149‑1 (MQ=255)                        gaTTGCCCATTAGGATGGCATCGAGTGCCCGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCCCGTCGCCGGCGGGCGGCTTGATGCCTGCCTGTTcc                                                                                                                       >  8:93731/1‑149 (MQ=255)                         aTTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCa                                                                                                                      <  5:288431/149‑1 (MQ=255)                            gCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAgg                                                                                                                    <  7:210655/148‑1 (MQ=255)                             cccATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAgg                                                                                                                    >  7:288448/1‑147 (MQ=255)                               cATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAg                                                                                                                 >  3:45635/1‑148 (MQ=255)                                aTTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGaa                                                                                                               <  7:93731/149‑1 (MQ=255)                                       tGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGcc                                                                                                        <  5:304958/149‑1 (MQ=255)                                                   cACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAggg                                                                                             >  3:138752/1‑148 (MQ=255)                                                   cACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAggg                                                                                             >  4:213753/1‑148 (MQ=255)                                                      gATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATg                                                                                         <  1:83500/149‑1 (MQ=255)                                                         ccGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCgg                                                                                      >  6:256534/1‑149 (MQ=255)                                                             ggCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGttct                                                                                  >  2:283107/1‑149 (MQ=255)                                                                  aCCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAgg                                                                                              <  6:112612/132‑1 (MQ=255)                                                                  aCCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAgg                                                                                              >  5:112612/1‑132 (MQ=255)                                                                  aCCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCa                                                                             >  6:54537/1‑149 (MQ=255)                                                                   ccAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCgg                                                                                      >  5:338702/1‑139 (MQ=255)                                                                   ccAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCgg                                                                                      <  6:338702/139‑1 (MQ=255)                                                                   ccAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCAt                                                                            >  1:213195/1‑149 (MQ=255)                                                                    cAATAGCGATTGGCATAGTGCCACCGCGCCTTGCCGAAAAGGCTGATCGGGGCGAACGGGATGGATACCACGTCGACGTCGGACGGCTTGATGCCTGACTGTTCCAGGCAGAACTTCGCCGACTCGTAGGTCATGCGGTTCTTCGCATg                                                                           >  2:225571/1‑149 (MQ=255)                                                                      aTAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCt                                                                         >  4:270337/1‑149 (MQ=255)                                                                       tAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCtt                                                                        <  8:288448/149‑1 (MQ=255)                                                                       tAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTCTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCtt                                                                        <  1:19162/149‑1 (MQ=255)                                                                        aGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTg                                                                       >  1:223926/1‑149 (MQ=255)                                                                               ggCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTcacgca                                                                >  2:217435/1‑149 (MQ=255)                                                                                gCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTcacgcac                                                               >  7:327897/1‑149 (MQ=255)                                                                                gCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTcacgcac                                                               >  5:107465/1‑149 (MQ=255)                                                                                 cATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTcacgcac                                                               >  3:62995/1‑148 (MQ=255)                                                                                  aTAGTGCCACCGCGCCTTGTCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTcacgcacg                                                              >  7:6667/1‑148 (MQ=255)                                                                                  aTAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGaa                                                                                                               >  5:249826/1‑99 (MQ=255)                                                                                  aTAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGaa                                                                                                               <  6:249826/99‑1 (MQ=255)                                                                                  aTAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTcacgcacg                                                              >  5:145989/1‑148 (MQ=255)                                                                                     gTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAg                                                           <  3:213753/148‑1 (MQ=255)                                                                                        ccACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGct                                                       <  1:217435/149‑1 (MQ=255)                                                                                             gcgcCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCtcttct                                                  <  6:107465/149‑1 (MQ=255)                                                                                             gcgcCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCtcttct                                                  <  8:327897/149‑1 (MQ=255)                                                                                                     ccGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTc                                                                                                           <  6:44996/84‑1 (MQ=255)                                                                                                     ccGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTcggcgg                                           <  3:270337/148‑1 (MQ=255)                                                                                                      cgtaCAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCTTGCTTGTCACGCACGAAGCGCTCTTCTTcggcggcg                                         <  7:295959/146‑1 (MQ=255)                                                                                                      cGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTcggcggcg                                         <  7:300842/149‑1 (MQ=255)                                                                                                       gAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGACACGCACGAAGCGCTCTTCTTcggcggcgg                                        <  1:31419/149‑1 (MQ=255)                                                                                                       gAACAGGCTGATCGGGGCGAACGGGAGGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTcggcggcgg                                        <  6:344218/149‑1 (MQ=255)                                                                                                             gCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTACTTCGGCGGCGGCAAt                                    >  1:323535/1‑147 (MQ=255)                                                                                                              cTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAg                                 <  4:138752/149‑1 (MQ=255)                                                                                                                 aTCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCtt                              <  4:62995/149‑1 (MQ=255)                                                                                                                         gAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCaa                      >  4:186976/1‑149 (MQ=255)                                                                                                                             gggATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTa                  <  5:256534/149‑1 (MQ=255)                                                                                                                               gATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACa                <  4:45635/149‑1 (MQ=255)                                                                                                                                aTGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAg               <  5:54537/149‑1 (MQ=255)                                                                                                                                   gCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAggg             <  6:145989/148‑1 (MQ=255)                                                                                                                                     tACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCcg          >  5:221506/1‑149 (MQ=255)                                                                                                                                     tACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCGTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCGTGCCGTCAATGTACAGGGCcg          >  4:239009/1‑149 (MQ=255)                                                                                                                                          cgtcgacgtcgGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCgg       >  6:157487/1‑147 (MQ=255)                                                                                                                                              gacgtcgGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAggg  <  2:223926/148‑1 (MQ=255)                                                                                                                                                cgtcgGACGGCTTGATGCCTGCCTGTTCCAGGCGGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTcgg                                              >  2:148102/1‑102 (MQ=255)                                                                                                                                                cgtcgGACGGCTTGATGCCTGCCTGTTCCAGGCGGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTcgg                                              <  1:148102/102‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  CTTGCGATAGCGCTTGTAGCGACGATTGCCCATTAGGATGGCATCGAGTGCACGATCCGGGGCGTACCAATAGCGCTTGGCATAGTGCCACCGCGCCTTGCCGAACAGGCTGATCGGGGCGAACGGGATGGCTACCACGTCGACGTCGGACGGCTTGATGCCTGCCTGTTCCAGGCAGAACTTCGCCGACTCGTAGGGCATGCGGTTCTTCGCATGCTTGTCACGCACGAAGCGCTCTTCTTCGGCGGCGGCAATCAGCTTGCCGTCAATGTACAGGGCCGCGGAAGGG  >  NC_002947/5627169‑5627457

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