Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,018,260	(A)8→7	coding (315/444 nt)	fliJ →	flagellar protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,018,253	0	A	.	92.3%	194.1 / 6.1	52	coding (308/444 nt)	fliJ	flagellar protein
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base . (34/14);  total (34/18)
Fisher's exact test for biased strand distribution p-value = 1.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGC  >  NC_000913/2018104‑2018401                                                                                                                                                      |                                                                                                                                                     cAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                                                                     >  1:74963/1‑151 (MQ=255)      gCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                                                                                >  1:689014/1‑151 (MQ=255)      gCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                                                                                >  1:583082/1‑151 (MQ=255)       cGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                                                                               <  2:880430/151‑2 (MQ=255)       cGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                                                                                >  2:747935/1‑150 (MQ=255)       cGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                                                                                >  1:573710/1‑150 (MQ=255)                   cGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAgg                                                                                                                                  <  2:850553/151‑1 (MQ=255)                    gggATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAgg                                                                                                                                  >  1:963411/1‑150 (MQ=255)                               cAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGacac                                                                                                                      <  1:1090496/151‑1 (MQ=255)                                  ccGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTg                                                                                                                    >  2:33402/1‑150 (MQ=255)                                       ggATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGaa                                                                                                              >  2:210696/1‑151 (MQ=255)                                        gATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGaa                                                                                                              <  1:1057847/150‑1 (MQ=255)                                         aTCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAAc                                                                                                             <  1:795589/150‑1 (MQ=255)                                         aTCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACg                                                                                                            <  1:602191/151‑1 (MQ=255)                                             aCTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCaa                                                                                                        >  1:405520/1‑151 (MQ=255)                                                aTCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATcc                                                                                                     <  2:806128/151‑1 (MQ=255)                                                aTCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATc                                                                                                      >  2:862167/1‑150 (MQ=255)                                                     cagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgg                                                                                                 >  1:50193/1‑150 (MQ=255)                                                      agTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGCGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcg                                                                                               >  2:783318/1‑151 (MQ=255)                                                       gTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcgg                                                                                              >  2:1007489/1‑151 (MQ=255)                                                           aTCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCa                                                                                            >  1:26807/1‑149 (MQ=255)                                                                       gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaaa                                                                              >  2:288035/1‑151 (MQ=255)                                                                       gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaa                                                                               <  2:758982/150‑1 (MQ=255)                                                                             gCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                                                                         >  1:639727/1‑150 (MQ=255)                                                                              ccATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                                                                         <  1:1260103/149‑1 (MQ=255)                                                                              ccATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCt                                                                        >  2:25200/1‑150 (MQ=255)                                                                                attccccAGCACCGCCAGAAACTTAAACAGGGGACGCAAAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCGGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa                                                                     <  1:288100/145‑1 (MQ=255)                                                                                 ttACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt                                                                    >  1:11604/1‑151 (MQ=255)                                                                                   aCTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa                                                                  >  1:1089288/1‑151 (MQ=255)                                                                                      cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGACAACCGCCTCGATCAGaa                                                               >  2:22624/1‑151 (MQ=255)                                                                                        gCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaa                                                             >  1:1216279/1‑151 (MQ=255)                                                                                         cATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAAAAGTTGGCGAG‑AAAAAAACAACGTTGGCAGGCCGGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaaa                                                            <  1:835477/151‑1 (MQ=255)                                                                                         cATCGCCAGCAACTTAATAAGTGGACGCAGAAAGTGGACATTGCCCTGAACAGTTGGCGA‑AAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGACCGGCAAGCCCCGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaa                                                              <  2:1251180/149‑1 (MQ=255)                                                                                            cGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgat                                                         <  1:989232/151‑1 (MQ=255)                                                                                              ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACTTTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCTGCACTGCTTGCCGAAAACCGCCCCGGTCAGAAAACgatgg                                                       >  2:419905/1‑151 (MQ=255)                                                                                                  cAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACATTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGa                                                   >  2:1251671/1‑151 (MQ=255)                                                                                                    aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGt                                                 >  2:182649/1‑151 (MQ=255)                                                                                                     cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGtt                                                >  1:57188/1‑151 (MQ=255)                                                                                                        aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCg                                              >  2:1100711/1‑150 (MQ=255)                                                                                                              gggACGCAGAAAGTTGACATTGCCCTGAACAGTGGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCCCGGCGGCCCTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAg                                         <  1:23304/148‑1 (MQ=255)                                                                                                                    cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGcgccg                                  >  2:12253/1‑150 (MQ=255)                                                                                                                     aaaaaTTTAACATTGCCCTGAACAGTGGGGGAG‑AAAAAACCACCGTTTGCAGGCCGGGCAGACACTGCAGGAACGGCAATCCCCGGCGGCCCTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCgccgcc                                <  1:395769/149‑1 (MQ=255)                                                                                                                           ttaaCATGGCCCAAAACGGTGGGCGA‑AAAAAAAAAAACGTTTCCGGGCCGGGAAGACACCGCGGAAACGGAAACCCCCGGCGGCACGGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAgg                          <  1:1188509/148‑1 (MQ=255)                                                                                                                           tggACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAgg                          <  1:685501/149‑1 (MQ=255)                                                                                                                              aCATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGaaa                       >  2:1131455/1‑151 (MQ=255)                                                                                                                                   gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGa                  >  2:1122618/1‑151 (MQ=255)                                                                                                                                   gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGa                  >  1:177567/1‑151 (MQ=255)                                                                                                                                   gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGa                  >  1:688362/1‑151 (MQ=255)                                                                                                                                   gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGa                  >  1:224072/1‑151 (MQ=255)                                                                                                                                   gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGa                  >  2:936316/1‑151 (MQ=255)                                                                                                                                   gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGa                  >  2:954013/1‑151 (MQ=255)                                                                                                                                   gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGa                  >  2:990724/1‑151 (MQ=255)                                                                                                                                     ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAAt                >  2:548844/1‑151 (MQ=255)                                                                                                                                     ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGA‑GAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATg               >  1:919142/1‑151 (MQ=255)                                                                                                                                     ccTGAACAGTGGGCGA‑AAAAAAAACAACTTTTGCAGGCCTGGCAGACACGGCAGGAACGGCAACCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAAt                <  2:405404/151‑1 (MQ=255)                                                                                                                                     ccTGAACAGTGGGCGA‑AAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAAAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAAt                <  1:991090/151‑1 (MQ=255)                                                                                                                                        gAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGAt             >  1:1119325/1‑151 (MQ=255)                                                                                                                                             gTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCg          >  1:915370/1‑149 (MQ=255)                                                                                                                                                     gaAAAAAACCACCGTTTGCAGGCCGGACAGACGCGACAGGAACGGAAACCCACGGGGGCCCTTCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCGTAgcgc  <  2:555413/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CAACAGCGATATGAGTGCCGGGATAACCAGCAACCGCTGGATCAACTATCAGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGC  >  NC_000913/2018104‑2018401

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑