BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,180,497	Δ1 bp	8.3%	coding (77/1161 nt)	ygiC →	ATP‑Grasp family ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,180,497	0	T	.	8.3%	781.4 / 59.4	230	coding (77/1161 nt)	ygiC	ATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand): ref base T (98/113); new base . (10/9); total (108/122)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

Only 100 of 234 total aligned reads displayed.

TGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAA  >  NC_000913/3180349‑3180644

tGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGGTCAATGGGTGGCTGGTACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAt                                                                                                                                                      <  4:326783/149‑1 (MQ=255)

       aaacaaaGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGc                                                                                                                                               <  2:1554500/149‑1 (MQ=255)

       aaacaaaGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGc                                                                                                                                               >  2:446183/1‑149 (MQ=255)

             aGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCg                                                                                                                                         <  2:492208/149‑1 (MQ=255)

                 cTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAt                                                                                                                                                      <  4:417921/132‑1 (MQ=255)

                  tATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑Tg                                                                                                                                   <  1:56265/149‑1 (MQ=255)

                    tGCAGCGTAGTTCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgt                                                                                                                                 <  3:62075/149‑1 (MQ=255)

                    tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgt                                                                                                                                 <  4:658880/149‑1 (MQ=255)

                    tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTAcatggt                                                                                                                                 <  1:1097156/149‑6 (MQ=255)

                    tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTAcatggt                                                                                                                                 >  2:1097156/1‑144 (MQ=255)

                     gCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgtg                                                                                                                                <  2:314058/149‑1 (MQ=255)

                     gCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTAcatggt                                                                                                                                 <  3:1107774/148‑6 (MQ=255)

                      cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGTGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgtgt                                                                                                                               <  3:365879/149‑1 (MQ=255)

                      cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTACATGgtgt                                                                                                                               <  2:15933/149‑1 (MQ=255)

                              tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                       <  3:638137/149‑1 (MQ=255)

                              tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                       <  1:688175/149‑1 (MQ=255)

                               gCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCt                                                                                                                      <  2:908769/149‑1 (MQ=255)

                               gCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCt                                                                                                                      <  2:876001/149‑1 (MQ=255)

                                cAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCtt                                                                                                                     >  4:1294821/1‑149 (MQ=255)

                                cAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCtt                                                                                                                     >  4:239709/1‑149 (MQ=255)

                                 aaCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTta                                                                                                                    <  4:1498742/149‑1 (MQ=255)

                                 aaCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTta                                                                                                                    <  4:947239/149‑1 (MQ=255)

                                       tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACaa                                                                                                              <  4:564948/149‑1 (MQ=255)

                                       tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACaa                                                                                                              <  2:138850/149‑1 (MQ=255)

                                          ctctTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGtt                                                                                                           >  2:953244/1‑149 (MQ=255)

                                           tcttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTg                                                                                                          >  1:738323/1‑149 (MQ=255)

                                            cttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGa                                                                                                         <  1:1578139/149‑1 (MQ=255)

                                              tctcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTa                                                                                                                                       >  4:428996/1‑118 (MQ=255)

                                              tctcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTa                                                                                                                                       <  3:428996/118‑1 (MQ=255)

                                               ctcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGAccc                                                                                                      <  4:603410/149‑1 (MQ=255)

                                                  gTTCAATGTGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCg                                                                                                   >  1:148632/1‑149 (MQ=255)

                                                  gTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTta                                                                                                                    >  2:539589/1‑132 (MQ=255)

                                                  gTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTta                                                                                                                    <  1:539589/132‑1 (MQ=255)

                                                  gTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTAcatg                                                                                                                                   >  2:65663/1‑114 (MQ=255)

                                                  gTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTAcatg                                                                                                                                   >  2:211970/1‑114 (MQ=255)

                                                  gTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCA‑GTACGGCTCGCCGTAcatg                                                                                                                                   <  1:65663/117‑4 (MQ=255)

                                                   ttCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGc                                                                                                  >  4:1033659/1‑149 (MQ=255)

                                                      aaTGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCa                                                                                               >  3:1260674/1‑149 (MQ=255)

                                                        tGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAgg                                                                                             >  2:104531/1‑149 (MQ=255)

                                                        tGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAgg                                                                                             >  2:638950/1‑149 (MQ=255)

                                                         gggTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGt                                                                                            <  4:117536/149‑1 (MQ=255)

                                                                  aTACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgt                                                                                                                                 >  1:958556/1‑103 (MQ=255)

                                                                  aTACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGgt                                                                                                                                 <  2:958556/103‑1 (MQ=255)

                                                                  aTACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                   <  2:851797/149‑1 (MQ=255)

                                                                  aTACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                   <  2:745411/149‑1 (MQ=255)

                                                                    aCCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTgg                                                                                 >  1:978170/1‑149 (MQ=255)

                                                                          ggAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGaagaag                                                                           <  4:162486/149‑1 (MQ=255)

                                                                           gAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGt                                                                          <  1:150531/149‑1 (MQ=255)

                                                                           gAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGGTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGt                                                                          <  3:86029/149‑1 (MQ=255)

                                                                             aaGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCa                                                                        >  4:480125/1‑149 (MQ=255)

                                                                                 gTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTc                                                                         >  2:996515/1‑144 (MQ=255)

                                                                                  tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgcc                                                                   >  4:88987/1‑149 (MQ=255)

                                                                                  tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgcc                                                                   >  1:225021/1‑149 (MQ=255)

                                                                                     gTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGaa                                                                >  3:965247/1‑149 (MQ=255)

                                                                                      tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                       >  1:343720/1‑93 (MQ=255)

                                                                                      tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGc                                                                                                                       <  2:343720/93‑1 (MQ=255)

                                                                                      tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGccc                                                                                                <  3:843428/116‑1 (MQ=255)

                                                                                      tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAAc                                                               >  1:283572/1‑149 (MQ=255)

                                                                                          aCCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCa                                                           <  3:448671/149‑1 (MQ=255)

                                                                                            cGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAcc                                                         >  1:636319/1‑149 (MQ=255)

                                                                                                cgcCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAt                                                     <  1:1400696/149‑1 (MQ=255)

                                                                                                   ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAg                                                                                                             >  4:1448158/1‑90 (MQ=255)

                                                                                                   ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAg                                                                                                             <  3:1448158/90‑1 (MQ=255)

                                                                                                      gACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAgg                                                                                             <  1:208536/103‑1 (MQ=255)

                                                                                                         tGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGaa                                            <  2:278465/149‑1 (MQ=255)

                                                                                                           gCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGaagg                                          >  1:1315578/1‑147 (MQ=255)

                                                                                                           gCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAg                                          >  3:57771/1‑149 (MQ=255)

                                                                                                             gTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtg                                        >  3:1071543/1‑149 (MQ=255)

                                                                                                              tGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                       <  3:1450578/149‑1 (MQ=255)

                                                                                                               gagaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                      <  1:381918/149‑1 (MQ=255)

                                                                                                               gagaAAGCCCACGAATACGCTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                      <  4:1071543/149‑1 (MQ=255)

                                                                                                                 gaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGtggtgg                                    >  2:1220932/1‑149 (MQ=255)

                                                                                                                 gaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTA‑AAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGa                                   >  1:557834/1‑149 (MQ=255)

                                                                                                                  aaaGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGa                                   >  3:1506614/1‑149 (MQ=255)

                                                                                                                    aGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaa                                 <  2:1363849/149‑1 (MQ=255)

                                                                                                                      cccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                               >  3:1249476/1‑149 (MQ=255)

                                                                                                                      cccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                               <  1:848682/149‑1 (MQ=255)

                                                                                                                      cccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                               <  3:522375/149‑1 (MQ=255)

                                                                                                                       ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAg                              >  2:129833/1‑149 (MQ=255)

                                                                                                                        cACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGt                             <  1:446183/149‑1 (MQ=255)

                                                                                                                           gAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                          <  2:776295/149‑1 (MQ=255)

                                                                                                                            aaTACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATc                         >  1:1195584/1‑149 (MQ=255)

                                                                                                                             aTACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                            <  4:151581/97‑1 (MQ=255)

                                                                                                                             aTACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                            <  2:1029349/97‑1 (MQ=255)

                                                                                                                             aTACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                            >  1:1029349/1‑97 (MQ=255)

                                                                                                                                   tttCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCGTGAAAGTGGTGGAAAAAGTGATCGCCAGCg                  <  1:1172198/149‑1 (MQ=255)

                                                                                                                                   tttCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                            <  4:404538/139‑1 (MQ=255)

                                                                                                                                   tttCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAg                    >  1:948382/1‑147 (MQ=255)

                                                                                                                                   tttCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGAGGAAAAAGTg                            >  3:404538/1‑139 (MQ=255)

                                                                                                                                     tCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt                >  2:247391/1‑149 (MQ=255)

                                                                                                                                       aaTTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGa              <  2:395721/149‑1 (MQ=255)

                                                                                                                                        aTTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGt                             <  2:1184143/133‑1 (MQ=255)

                                                                                                                                           ttCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtg          <  3:1294821/149‑1 (MQ=255)

                                                                                                                                           ttCACACCA‑GTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtg          >  2:768155/1‑149 (MQ=39)

                                                                                                                                             cacaCCA‑GTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgat        >  3:1142044/1‑149 (MQ=255)

                                                                                                                                                acCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                       <  2:1497123/115‑1 (MQ=255)

                                                                                                                                                acCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAc     >  3:786468/1‑149 (MQ=255)

                                                                                                                                                acCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAc     >  2:1135401/1‑149 (MQ=255)

                                                                                                                                                acCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAc     >  1:861341/1‑149 (MQ=255)

                                                                                                                                                  cATGTACGGCGAGCCGTAC‑TGGCGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                   >  3:32051/1‑69 (MQ=255)

TGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAC‑TGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAA  >  NC_000913/3180349‑3180644

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`