Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,922,436	A→T	6.0%	L684H (CTC→CAC)	gcvP‑II ←	glycine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,922,436	0	A	T	6.0%	200.8 / 2.5	66	L684H (CTC→CAC)	gcvP‑II	glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (25/37);  new base T (2/2);  total (27/39)
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.69e-01

GACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGGGTAGGTAATCATCAGGCACGACAGACGCTCCCCGGCCGCATGCGCCTTGG  >  NC_002947/5922288‑5922580                                                                                                                                                     |                                                                                                                                                 gacgaACGGCTTGAGGTGCGCGCGGATACCGATCGG‑CCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAgg                                                                                                                                                <  2:907990/150‑1 (MQ=255)    gaACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGtt                                                                                                                                              <  2:1633462/150‑1 (MQ=255)          cTTGAGGTGCTAGAGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCAcc                                                                                                                                        <  1:744628/150‑1 (MQ=255)                gTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCAt                                                                                                                                  >  2:1044911/1‑150 (MQ=255)                          aTACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTgg                                                                                                                        >  2:1052625/1‑150 (MQ=255)                              cGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGgccgc                                                                                                                    >  1:224707/1‑150 (MQ=255)                                aTCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGccgccg                                                                                                                  >  2:1369244/1‑150 (MQ=255)                                   ggCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTAc                                                                                                               >  2:604889/1‑150 (MQ=255)                                    gCCCCATCCCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGGGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGACGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACt                                                                                                              <  2:222544/150‑1 (MQ=255)                                     ccccATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGCCCGCCGTACTg                                                                                                             <  2:1511455/150‑1 (MQ=255)                                       ccATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGt                                                                                                           >  1:1304070/1‑150 (MQ=255)                                       ccATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGt                                                                                                           >  2:1134093/1‑150 (MQ=255)                                       ccATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGt                                                                                                           >  2:1133802/1‑150 (MQ=255)                                       ccATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGt                                                                                                           >  2:150225/1‑150 (MQ=255)                                        cATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTg                                                                                                          <  2:1499760/150‑1 (MQ=255)                                             ccGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTTCTGGtccgggc                                                                                                     <  2:1202557/150‑7 (MQ=255)                                               gggCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACcaca                                                                                                   <  1:1052625/150‑1 (MQ=255)                                                    gccgccgccgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGcacca                                                                                                     <  2:1322771/143‑1 (MQ=255)                                                    gccgccgccgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGcacca                                                                                                     >  1:1322771/1‑143 (MQ=255)                                                    gccgccgccgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCa                                                                                              <  2:484612/150‑1 (MQ=255)                                                    gccgccgccgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCa                                                                                              <  2:10453/150‑1 (MQ=255)                                                      cgccgccgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGa                                                                                            <  1:948103/150‑1 (MQ=255)                                                        ccgccgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGTGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGAtt                                                                                          >  2:608401/1‑150 (MQ=255)                                                        ccgccgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGTGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGAtt                                                                                          >  2:608370/1‑150 (MQ=255)                                                          gccgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGTGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTc                                                                                        <  1:608401/150‑1 (MQ=255)                                                          gccgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGTGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTc                                                                                        <  1:608370/150‑1 (MQ=255)                                                            cgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAg                                                                                                                                                 <  1:80186/91‑1 (MQ=255)                                                            cgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAg                                                                                                                                                 >  2:80186/1‑91 (MQ=255)                                                            cgTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTcgc                                                                                      <  2:64965/150‑1 (MQ=255)                                                                gggAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTAt                                                                                  >  2:200985/1‑150 (MQ=255)                                                                  gAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATcc                                                                                >  1:1440761/1‑150 (MQ=255)                                                                      gCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCctc                                                                            <  1:1369244/150‑1 (MQ=255)                                                                        aGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTcgc                                                                                      >  2:376275/1‑138 (MQ=255)                                                                        aGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTcgc                                                                                      <  1:376275/138‑1 (MQ=255)                                                                          aaCGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCtt                                                                          <  2:480844/148‑1 (MQ=255)                                                                          aaCGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCg                                                                        <  2:771313/150‑1 (MQ=255)                                                                          aaCGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTAAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCtt                                                                          >  1:480844/1‑148 (MQ=255)                                                                              tCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTg                                                                                                             <  2:801964/109‑1 (MQ=255)                                                                              tCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTg                                                                                                             >  1:801964/1‑109 (MQ=255)                                                                               cTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTacac                                                                   <  2:1304070/150‑1 (MQ=255)                                                                                  gtgCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGcc                                                                >  1:1090388/1‑150 (MQ=255)                                                                                    gCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGca                                                                                                        <  1:1408878/108‑1 (MQ=255)                                                                                    gCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGca                                                                                                        >  2:1408878/1‑108 (MQ=255)                                                                                     cAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTagc                                                                                 >  2:839460/1‑128 (MQ=255)                                                                                     cAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATc                                                                                 <  1:839460/130‑1 (MQ=255)                                                                                     cAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTg                                                             <  1:790526/150‑1 (MQ=255)                                                                                      aGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTgg                                                            <  2:906801/150‑1 (MQ=255)                                                                                         ttCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTc                                                         <  1:1133802/150‑1 (MQ=255)                                                                                         ttCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTc                                                         <  1:1134093/150‑1 (MQ=255)                                                                                            aTGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAg                                                      <  1:804112/150‑1 (MQ=255)                                                                                              gtgtgAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGTTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCATTGCATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAggg                                                    <  1:1117360/150‑1 (MQ=255)                                                                                                gtgAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCa                                                                                              <  1:152438/106‑1 (MQ=255)                                                                                                gtgAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCa                                                                                              >  2:152438/1‑106 (MQ=255)                                                                                                 tgAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGggta                                                 >  2:1166380/1‑150 (MQ=255)                                                                                                   aaaCGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATc                                                                                 <  1:1472658/116‑1 (MQ=255)                                                                                                   aaaCGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATc                                                                                 >  2:1472658/1‑116 (MQ=255)                                                                                                   aaaCGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGggtagg                                               <  1:604889/150‑1 (MQ=255)                                                                                                           gCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGGGTAGGTAatcatc                                       <  1:150225/150‑1 (MQ=255)                                                                                                             cccAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGGGTAGGTAATCATCAg                                     <  2:1440761/150‑1 (MQ=255)                                                                                                               cAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGGGTAGGTAATCATCAGGc                                   <  1:1216299/150‑1 (MQ=255)                                                                                                                    tCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGGGTAGGTAATCATCAGGCACGAc                              >  2:800714/1‑150 (MQ=255)                                                                                                                      tGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGGGTAGGTAATCATCAGGCACGACAg                            <  2:1090388/150‑1 (MQ=255)                                                                                                                             cgcgcCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGGGTAGGTAATCATCAGGCACGACAGACGCTcc                     <  1:1166380/150‑1 (MQ=255)                                                                                                                                gcCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGGGTAGGTAATCATCAGGCACGACAGACGCTCCCCg                  >  2:573558/1‑150 (MQ=255)                                                                                                                                        cACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGGGTAGGTAATCATCAGGCACGACAGACGCTCCCCGGCCGCATg          <  1:573558/150‑1 (MQ=255)                                                                                                                                                cGTTGAGGTTGGCACCGTCCATGTAAACCTCGCCGTCGTCCTGGTACAAACCAGCCCTGATTTCGCGTATCGCCGCTTCGTACACGCCGTGGGTCGAGGGGTAGGTAATCATCAGGCACGACAGACGCTCCCCGGCCGCATGCGCCTTgg  >  1:30875/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GACGAACGGCTTGAGGTGCGCGCGGATACCGATCGGCCCCATACCCGGGCCGCCGCCGCCGTGGGGAATGCAGAACGTCTTGTGCAGGTTCATGTGTGAAACGTCGGCCCCAATGTCTGCGGGGCGCGCCAGGCCCACCTGGGCGTTGAGGTTGGCACCGTCCATGTAAACCTGGCCGCCGTACTGGTGCACCACATCGCAGATTTCGCGTATCCCCTCTTCGTACACGCCGTGGGTCGAGGGGTAGGTAATCATCAGGCACGACAGACGCTCCCCGGCCGCATGCGCCTTGG  >  NC_002947/5922288‑5922580

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