Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,294,591	C→T	5.2%	G105D (GGC→GAC)	ccmF ←	heme lyase, CcmF subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,294,591	0	C	T	5.2%	244.9 / 2.2	76	G105D (GGC→GAC)	ccmF	heme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base C (40/32);  new base T (2/2);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

CAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCA  >  NC_000913/2294452‑2294738                                                                                                                                            |                                                                                                                                                    aaaCGTGCGAGATAACGGGTTTTAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCAGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCtcatca                                                                                                                                            <  2:346479/148‑1 (MQ=255)  aaCGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACAATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATcag                                                                                                                                           >  7:345434/1‑149 (MQ=255)    cGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATcagca                                                                                                                                         <  1:254511/149‑1 (MQ=255)    cGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATcagca                                                                                                                                         <  6:77083/149‑1 (MQ=255)    cGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATcagca                                                                                                                                         <  7:363113/149‑1 (MQ=255)         gagagaACGGGTTAGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGCCCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCAcc                                                                                                                                    >  1:260797/1‑149 (MQ=255)              aaCGGGTTAGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCAcag                                                                                                                               >  7:235599/1‑149 (MQ=255)                cGGGTTAGAGGTAAAGAGAATGAACCGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCAcag                                                                                                                               >  8:88773/1‑147 (MQ=255)                        agGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAg                                                                                                                     >  7:279288/1‑149 (MQ=255)                          gTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCCTCAGCAGCACCCACAGCAGTAGCGAGcc                                                                                                                   >  2:271501/1‑149 (MQ=255)                              agagAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCa                                                                                                               <  1:6435/149‑1 (MQ=255)                                 gaATTAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGGACACGGGCCACAAAATACAGCGGAATACGCTTACTACAAATCGCCACCGCAAAGGTCAAGCCGCTCATCAGCAGAACCCACAGCAGTAGCGAGCCTTCATgc                                                                                                            >  6:304284/1‑149 (MQ=255)                                    tGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATgcgcg                                                                                                         >  7:205774/1‑149 (MQ=255)                                       aCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGccc                                                                                                      >  1:175293/1‑149 (MQ=255)                                        cAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGccc                                                                                                      >  2:28574/1‑148 (MQ=255)                                         aGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGTAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCa                                                                                                    <  2:332467/149‑1 (MQ=255)                                                 gCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAgccg                                                                                            >  6:150092/1‑149 (MQ=255)                                                     acacTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCAc                                                                                        >  7:186338/1‑149 (MQ=255)                                                       acTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCGCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCAcgc                                                                                      <  1:271501/149‑1 (MQ=255)                                                           aCCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGa                                                                                    <  8:279288/147‑1 (MQ=255)                                                             cATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCCTCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATAc                                                                                 >  4:225104/1‑148 (MQ=255)                                                             cATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCCCCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATAc                                                                                 >  6:201382/1‑148 (MQ=255)                                                             cATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATAcc                                                                                >  6:196335/1‑149 (MQ=255)                                                              aTCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACca                                                                               <  2:260797/149‑1 (MQ=255)                                                                ccccATTATCGCCAGTACACGGGCCACAATATCCTGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACcaca                                                                             >  8:205184/1‑149 (MQ=255)                                                                  ccATTATCGCCAGTACACGGGCCATAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCGCATCAGCAGCACCCACAGCAGTAGCGCGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATAccccccc                                                                           >  5:111787/1‑144 (MQ=255)                                                                  ccATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTc                                                                                                                <  6:252386/112‑1 (MQ=255)                                                                  ccATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTc                                                                                                                >  5:252386/1‑112 (MQ=255)                                                                   cATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCg                                                                          >  4:97502/1‑149 (MQ=255)                                                                     ttATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCgg                                                                         >  8:35632/1‑148 (MQ=255)                                                                        tCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGc                                                                     >  4:41521/1‑149 (MQ=255)                                                                           ccAGTACACTGGACACAATATCCAGCGGAATACGCTGACTAAAAATCGCCCCCGCAAAGGTCCAGCCGCTGCTCAGCAGCTCCCACAGCTGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTgg                                                                  <  4:177961/149‑1 (MQ=255)                                                                           ccAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTgg                                                                  <  5:150092/149‑1 (MQ=255)                                                                               tACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCCCCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTAt                                                              >  5:99473/1‑149 (MQ=255)                                                                               tACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTAt                                                              >  1:11134/1‑149 (MQ=255)                                                                                 cacGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTg                                                            <  6:218776/149‑1 (MQ=255)                                                                                 cacGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTg                                                            <  6:99473/149‑1 (MQ=255)                                                                                   cGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGa                                                          <  3:41521/149‑1 (MQ=255)                                                                                    gggCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCCCCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTgg                                                                  >  3:220860/1‑140 (MQ=255)                                                                                    gggCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTgg                                                                  <  4:220860/140‑1 (MQ=255)                                                                                      gCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCCCCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGtt                                                       >  2:248995/1‑149 (MQ=255)                                                                                       ccACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTg                                                      >  3:168541/1‑149 (MQ=255)                                                                                           aaTATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTgg                                                  >  8:216862/1‑149 (MQ=255)                                                                                               tCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGAAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCaa                                               >  5:336726/1‑148 (MQ=255)                                                                                                 cAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTg                                                            >  2:298493/1‑133 (MQ=255)                                                                                                  aGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACa                                             >  2:17038/1‑147 (MQ=255)                                                                                                  aGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATa                                           <  8:186338/149‑1 (MQ=255)                                                                                                   gCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAATCTTGGTATTGGAGTTGCTGGCAACATAg                                          <  2:175293/149‑1 (MQ=255)                                                                                                   gCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAg                                          <  7:88773/149‑1 (MQ=255)                                                                                                    cGGAATACGCTGACTAAAAATCGCCACCCCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATTCGAGCCCCAGGTAGCCTCCACGCGATACCACAACGGAAGCTGGGTATTGTAGTTGCTGGCAACATAgg                                         <  5:304284/149‑1 (MQ=255)                                                                                                       aaTACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGa                                      <  7:205184/149‑1 (MQ=255)                                                                                                         tACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGAGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACg                                    <  8:345434/149‑1 (MQ=255)                                                                                                              tGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTCGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCCCGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTg                                 >  4:252145/1‑147 (MQ=255)                                                                                                               gACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAg                              >  2:4621/1‑149 (MQ=255)                                                                                                                  tAAAAATCGCCACCGCAAAGGTCCAGTCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTgg                                                                  <  5:84692/110‑1 (MQ=255)                                                                                                                  tAAAAATCGCCACCGCAAAGGTCCAGTCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTgg                                                                  <  5:323607/110‑1 (MQ=255)                                                                                                                  tAAAAATCGCCACCGCAAAGGTCCAGTCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTgg                                                                  >  6:84692/1‑110 (MQ=255)                                                                                                                  tAAAAATCGCCACCGCAAAGGTCCAGTCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGAGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTgg                                                                  >  6:323607/1‑110 (MQ=255)                                                                                                                  tAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCCCACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCg                           >  2:193937/1‑149 (MQ=255)                                                                                                                  tAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCg                           >  4:156248/1‑149 (MQ=255)                                                                                                                    aaaaTCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGtt                         <  8:205774/149‑1 (MQ=255)                                                                                                                      aaTCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGa                       <  1:4621/149‑1 (MQ=255)                                                                                                                        tCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACTGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTgg                                                                  <  6:224520/104‑1 (MQ=255)                                                                                                                        tCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTgg                                                                  >  5:224520/1‑104 (MQ=255)                                                                                                                          gCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCa                    <  1:28574/148‑1 (MQ=255)                                                                                                                           ccACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCAAGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTCAAGGTGAAGTCGTTGACCAc                   >  1:364927/1‑148 (MQ=255)                                                                                                                             aCCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGaa                <  7:216862/149‑1 (MQ=255)                                                                                                                             aCCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGa                 >  4:174392/1‑148 (MQ=255)                                                                                                                             aCCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGa                 >  4:295777/1‑148 (MQ=255)                                                                                                                              ccGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGaaa               <  5:201382/149‑1 (MQ=255)                                                                                                                              ccGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGaaa               <  5:196335/149‑1 (MQ=255)                                                                                                                              ccGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGaaa               <  3:225104/149‑1 (MQ=255)                                                                                                                              ccGCAAAGCTCCAGCCGCTCCTCAGCAGCACCCACAGCATTAGCGAGCCTTCATGCTCTCCCTATGTAGCGGCCACGCGATAACACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGaaa               <  6:305374/149‑1 (MQ=255)                                                                                                                                  aaaGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCAt           >  7:266153/1‑149 (MQ=255)                                                                                                                                         cAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAg    <  1:17038/149‑1 (MQ=255)                                                                                                                                           gCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCa  >  5:198557/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                    CAACGTGCGAGAGAACGGGTTAGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGCACCCACAGCAGTAGCGAGCCTTCATGCGCGCCCCAGGTAGCCGCCACGCGATACCACACCGGAAGCTGGGTATTGGAGTTGCTGGCAACATAGGTGACGGTGAAGTCGTTGACCACGAAAGCATTGACCAGCA  >  NC_000913/2294452‑2294738

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