Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	276,021	A→G	pseudogene (53/249 nt)	BW25113_RS01360 ←	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	276,021	0	A	G	93.8%	263.7 / 7.4	81	pseudogene (53/249 nt)	BW25113_RS01360	transposase
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base G (39/37);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

GGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGATATATAC  >  NZ_CP009273/275891‑276152                                                                                                                                   |                                                                                                                                    ggCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATAAACAGCCTGTTTg                                                                                                                             >  2:122209/1‑139 (MQ=255)       aGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCt                                                                                                                       >  1:841516/1‑139 (MQ=255)       aGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCt                                                                                                                       >  2:657064/1‑139 (MQ=255)          aTCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCa                                                                                                                    <  1:574089/139‑1 (MQ=255)                   ttCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTaaa                                                                                                            <  1:199972/138‑1 (MQ=255)                   ttCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAAc                                                                                                           >  1:64911/1‑139 (MQ=255)                   ttCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAAc                                                                                                           >  2:754532/1‑139 (MQ=255)                   ttCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAAc                                                                                                           >  2:152127/1‑139 (MQ=255)                        tattaGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCtt                                                                                                       <  1:494845/138‑1 (MQ=255)                             gTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTTTTTGACTGCTTCAGTTTTAAACTCTTCGGGa                                                                                                  <  1:122209/138‑1 (MQ=255)                             gTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGa                                                                                                  >  1:717358/1‑138 (MQ=255)                             gTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGa                                                                                                  >  2:386019/1‑138 (MQ=255)                             gTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGa                                                                                                  >  1:264914/1‑138 (MQ=255)                             gTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGa                                                                                                  >  2:623041/1‑138 (MQ=255)                                gaagaaTCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAAc                                                                                              >  1:187084/1‑139 (MQ=255)                                gaagaaTCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAAc                                                                                              >  2:319222/1‑139 (MQ=255)                                    aaTCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCtt                                                                                          <  2:195584/139‑1 (MQ=255)                                      tCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTAc                                                                                        <  1:319222/139‑1 (MQ=255)                                        cGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCg                                                                                      >  2:565488/1‑139 (MQ=255)                                        cGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCg                                                                                      >  1:425588/1‑139 (MQ=255)                                              gTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATg                                                                                <  2:64911/139‑1 (MQ=255)                                               tACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATgg                                                                               >  1:334684/1‑139 (MQ=255)                                               tACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGC‑TACCGCTCCTggg                                                                              >  1:428629/1‑139 (MQ=255)                                                    cttTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCAc                                                                           <  1:394282/138‑1 (MQ=255)                                                    cttTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCAc                                                                           <  2:841516/138‑1 (MQ=255)                                                       tATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACct                                                                         <  2:264914/137‑1 (MQ=255)                                                       tATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACctct                                                                       >  1:325504/1‑139 (MQ=255)                                                         tCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACctctct                                                                     >  1:147673/1‑139 (MQ=255)                                                              gCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAg                                                                >  1:233049/1‑139 (MQ=255)                                                                        tgtgGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTaaa                                                      >  2:813216/1‑139 (MQ=255)                                                                        tgtgGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTaaa                                                      >  2:855644/1‑139 (MQ=255)                                                                         gtgggtggTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAAt                                                     >  2:715911/1‑139 (MQ=255)                                                                         gtgggtggTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAAt                                                     >  1:179457/1‑139 (MQ=255)                                                                          tgggtggTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTaaa                                                      <  2:492737/137‑1 (MQ=255)                                                                          tgggtggTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTaaa                                                      <  2:680661/137‑1 (MQ=255)                                                                          tgggtggTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAAt                                                     >  1:44511/1‑138 (MQ=255)                                                                          tgggtggTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTaaa                                                      <  1:623041/137‑1 (MQ=255)                                                                               ggtgATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCg                                               >  1:769521/1‑139 (MQ=255)                                                                                gtgATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCg                                               >  2:349959/1‑138 (MQ=255)                                                                                gtgATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCg                                               >  2:292581/1‑138 (MQ=255)                                                                                 tgATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCg                                               >  2:210937/1‑137 (MQ=255)                                                                                      tCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGgtg                                         <  1:420856/138‑1 (MQ=255)                                                                                      tCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGgtg                                         <  2:234963/138‑1 (MQ=255)                                                                                      tCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGgtg                                         <  1:243515/138‑1 (MQ=255)                                                                                              tgtTGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTaaa                                >  1:393475/1‑139 (MQ=255)                                                                                                ttGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGtt                                   >  1:224017/1‑134 (MQ=255)                                                                                                ttGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGtt                                   <  2:224017/134‑1 (MQ=255)                                                                                                ttGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAAcc                              >  1:530471/1‑139 (MQ=255)                                                                                                ttGCAACGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAAcc                              >  2:271558/1‑139 (MQ=255)                                                                                                      cGCTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGc                        >  1:748155/1‑139 (MQ=255)                                                                                                        cTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACct                                                                         <  1:404200/88‑1 (MQ=255)                                                                                                        cTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACct                                                                         >  2:404200/1‑88 (MQ=255)                                                                                                        cTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGa                      >  2:704193/1‑139 (MQ=255)                                                                                                        cTGGCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGa                      >  1:684405/1‑139 (MQ=255)                                                                                                          ggCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAg                                            <  2:167278/115‑1 (MQ=255)                                                                                                          ggCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAg                                            >  1:167278/1‑115 (MQ=255)                                                                                                           gCAACAGAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGAtt                    >  2:16838/1‑138 (MQ=255)                                                                                                                 gAATAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAAtt             <  2:44511/139‑1 (MQ=255)                                                                                                                  aaTAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTg            <  1:130615/139‑1 (MQ=255)                                                                                                                  aaTAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTg            >  1:626453/1‑139 (MQ=255)                                                                                                                    tAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCg          <  1:152127/139‑1 (MQ=255)                                                                                                                    tAACCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCg          <  1:210937/139‑1 (MQ=255)                                                                                                                      aCCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGat        <  2:334684/139‑1 (MQ=255)                                                                                                                      aCCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGat        <  1:406712/139‑1 (MQ=255)                                                                                                                      aCCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGat        <  2:670238/139‑1 (MQ=255)                                                                                                                      aCCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGat        <  2:393475/139‑1 (MQ=255)                                                                                                                      aCCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGa         <  2:425588/138‑1 (MQ=255)                                                                                                                      aCCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGa         <  2:717358/138‑1 (MQ=255)                                                                                                                      aCCGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGa         <  1:565488/138‑1 (MQ=255)                                                                                                                       ccGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGc                                                               >  1:374492/1‑83 (MQ=255)                                                                                                                       ccGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGc                                                               <  2:374492/83‑1 (MQ=255)                                                                                                                        cGAGATCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGatat      <  1:813216/139‑1 (MQ=255)                                                                                                                         gagaTCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTaaa                                                      >  2:582493/1‑90 (MQ=255)                                                                                                                         gagaTCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTaaa                                                      <  1:582493/90‑1 (MQ=255)                                                                                                                            aTCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCtt                                                                                                       >  1:432356/1‑38 (MQ=255)                                                                                                                            aTCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCtt                                                                                                       <  2:432356/38‑1 (MQ=255)                                                                                                                            aTCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGAtatata   >  1:529090/1‑138 (MQ=255)                                                                                                                            aTCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGATATatac  <  1:450147/139‑1 (MQ=255)                                                                                                                            aTCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGATATatac  <  1:233792/139‑1 (MQ=255)                                                                                                                            aTCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGATATatac  <  2:530471/139‑1 (MQ=255)                                                                                                                            aTCAACAGCCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGATATatac  <  2:529090/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                    GGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTTGACTGCTTCAGTTTTAAACTCTTCGGGATAACGCTTACCGCTCATGGGCACCTCTCTTTAAGCCATTTTAAATGACTCGGAGGTGTCTGTTAAACCCGTGGCGATTCAAAATTGCGATATATAC  >  NZ_CP009273/275891‑276152

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