BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	566,531	3 bp→TGT	100%	pseudogene (156‑158/312 nt)	peaD →	pseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	566,531	0	G	T	97.2%	6137.1 / 72.0	2121	pseudogene (156/312 nt)	peaD	pseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand): ref base G (37/29); new base T (1105/948); total (1143/978)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
*	NC_000913	566,532	0	A	G	97.9%	5949.6 / 35.2	2121	pseudogene (157/312 nt)	peaD	pseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand): ref base A (25/21); new base G (1104/951); total (1143/978)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
*	NC_000913	566,533	0	C	T	97.2%	5920.9 / 56.8	2121	pseudogene (158/312 nt)	peaD	pseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand): ref base C (37/24); new base T (1105/947); total (1143/978)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 2123 total aligned reads displayed.

CGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGT  >  NC_000913/566388‑566681

|||

                 tCCGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAAtt                                                                                                                                   <  6:257923/148‑1 (MQ=255)

                   cGGATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAAc                                                                                                                                >  6:44180/1‑149 (MQ=255)

                    ggATGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAAcc                                                                                                                               <  6:66717/149‑1 (MQ=255)

                       tGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCgtg                                                                                                                            <  4:141647/149‑1 (MQ=255)

                       tGCGGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCgtg                                                                                                                            <  6:148563/149‑1 (MQ=255)

                         cggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCgtggt                                                                                                                          >  7:200370/3‑149 (MQ=255)

                          ggagTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACAgtggt                                                                                                                          >  1:156073/2‑148 (MQ=255)

                          agagTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCgtggt                                                                                                                          >  2:144418/1‑148 (MQ=255)

                             gTCTTATCCGTGGAAATCAAACGCGCGCTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTgagg                                                                                                                      <  3:178519/149‑2 (MQ=255)

                             gTCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCCGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTgagg                                                                                                                      <  3:43350/149‑2 (MQ=255)

                              tCTTATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTgaggc                                                                                                                     >  7:145278/1‑147 (MQ=255)

                                ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTgaggc                                                                                                                     >  2:38336/1‑145 (MQ=255)

                                ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                                                                                                                   >  1:235666/1‑149 (MQ=255)

                                ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                                                                                                                   >  2:237879/1‑149 (MQ=255)

                                ttATCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGa                                                                                                                   >  7:313485/1‑149 (MQ=255)

                                  aTCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGAtc                                                                                                                 >  2:169748/1‑148 (MQ=255)

                                   tCCGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGGGAGGCGagcc                                                                                                                >  1:92675/1‑146 (MQ=255)

                                     cGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCt                                                                                                              >  4:111038/1‑149 (MQ=255)

                                     cGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCt                                                                                                              <  6:202849/149‑1 (MQ=255)

                                       tGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGa                                                                                                            <  2:313527/149‑1 (MQ=255)

                                         gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                                                                                          >  8:343071/1‑149 (MQ=255)

                                         gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                                                                                          <  2:56612/149‑1 (MQ=255)

                                         gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                                                                                          <  7:5476/149‑1 (MQ=255)

                                         caaaTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTAACTGATCCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                                                                                          <  2:8883/148‑1 (MQ=255)

                                           aaTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACAGTGGTGAGGCGATCCCTGAAc                                                                                                          >  7:388381/1‑147 (MQ=255)

                                            aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAg                                                                                                       >  1:59385/1‑149 (MQ=255)

                                            aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAg                                                                                                       >  8:315196/1‑149 (MQ=255)

                                            aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAg                                                                                                       >  1:369681/1‑149 (MQ=255)

                                            aTCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACAGTGGTGAGGCGATCCCTGAACCAg                                                                                                       <  3:265848/149‑1 (MQ=255)

                                             tcaaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAAt                                                                                                                                    >  8:366077/4‑119 (MQ=255)

                                             tcaaaCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGt                                                                                                      <  5:280305/146‑1 (MQ=255)

                                               aaaCGCGCACTACTGGCTGGTTACCAACCTGTATGAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTa                                                                                                     >  5:114815/2‑148 (MQ=255)

                                                 aCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaaaa                                                                                                  >  8:234286/1‑149 (MQ=255)

                                                 aCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaaaa                                                                                                  >  3:328151/1‑149 (MQ=255)

                                                 aCGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTaaaa                                                                                                  >  1:210653/1‑149 (MQ=255)

                                                     gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGCTCCCTGAACCAGTAAaacaac                                                                                              >  1:101471/1‑149 (MQ=255)

                                                     gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGaa                                                                                                           <  8:75621/136‑1 (MQ=255)

                                                      cactACTGGCTGGTTACCAACCTGTATCAGAACATGCTGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACt                                                                                             <  5:365927/146‑1 (MQ=255)

                                                           ctggctggTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCt                                                                                         <  8:143788/148‑1 (MQ=255)

                                                           ctggctggTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTg                                                                                        <  1:125119/149‑1 (MQ=255)

                                                            tggttggtTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑T‑GACCTGTATGACAGCGCGaa                                                                                                                                     >  4:25234/1‑103 (MQ=255)

                                                              gctggttACCAACGTGTATCAGAACTTGCGGGCCATTGCGCTTACTGATGCGGTATTCCGCCGTAAGGCCAAAGATGAGC‑T‑TGTCCATATGACTGCGAGATTTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCATGtca                                                                                     <  8:498/147‑3 (MQ=255)

                                                              gctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGtca                                                                                     >  5:307784/3‑147 (MQ=255)

                                                              gctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGtca                                                                                     <  1:184602/147‑3 (MQ=255)

                                                              gctggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGtca                                                                                     >  6:87169/3‑147 (MQ=255)

                                                               ctggttACCAACCTGTATCAGAACATGCTGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCAt                                                                                    <  2:257233/148‑1 (MQ=255)

                                                                tggttACCAACCTGTATCAGACCATGCGGGCCAATGCGCTTACTGAAGCGGAATTACGCCGTAAGGCCGCAGATGAGCTT‑GTCCAG‑ATGACTGCGAGAATTAACCGAGGAGAGGCGATCCCTGAACCAGTAAAACAACTTCCTCTCATg                                                                                   >  2:124300/1‑149 (MQ=255)

                                                                tggttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATg                                                                                   >  7:95580/1‑149 (MQ=255)

                                                                  gttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                                                                                          >  5:86195/1‑124 (MQ=255)

                                                                  gttACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATggg                                                                                 <  3:208531/149‑1 (MQ=255)

                                                                       cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTa                                                                            <  7:73218/149‑1 (MQ=255)

                                                                          cctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAc                                                                         <  4:34036/147‑1 (MQ=255)

                                                                            tGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TAGTCCAG‑ATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACct                                                                       >  3:237316/1‑149 (MQ=255)

                                                                            tGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACct                                                                       >  7:229113/1‑149 (MQ=255)

                                                                            tGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACct                                                                       >  4:175710/1‑149 (MQ=255)

                                                                             gTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACctc                                                                      >  7:111879/1‑149 (MQ=255)

                                                                              tatcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTgag                                                                                                                       <  6:235737/96‑1 (MQ=255)

                                                                                tcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGCTCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTa                                                                    >  1:306488/2‑148 (MQ=255)

                                                                                tcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTa                                                                    >  5:334714/2‑148 (MQ=255)

                                                                                 cAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTaaa                                                                  >  4:201300/1‑149 (MQ=255)

                                                                                  aGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTaaa                                                                  <  8:202793/148‑1 (MQ=255)

                                                                                  aGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAAt                                                                 >  1:337332/1‑149 (MQ=255)

                                                                                   gAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATggg                                                                                 >  4:8491/1‑132 (MQ=255)

                                                                                   gAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATc                                                                >  6:317795/1‑149 (MQ=255)

                                                                                   gAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATc                                                                <  1:281404/149‑1 (MQ=255)

                                                                                     aCATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGt                                                              <  4:282340/149‑1 (MQ=255)

                                                                                      cATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg                                                             <  3:355271/149‑1 (MQ=255)

                                                                                      cATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg                                                             <  2:202987/149‑1 (MQ=255)

                                                                                      cATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTg                                                             <  5:167707/149‑1 (MQ=255)

                                                                                       aTGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGGGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGt                                                              >  5:197100/1‑147 (MQ=255)

                                                                                        tGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTaa                                                                   >  6:339698/1‑141 (MQ=255)

                                                                                        tGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGca                                                           <  8:179595/149‑1 (MQ=255)

                                                                                          cGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGcaca                                                         <  8:115475/149‑1 (MQ=255)

                                                                                           gggCCAATGCGCTTACTGATTCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAg                                                        <  6:287016/149‑1 (MQ=255)

                                                                                               cAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGctc                                                    >  3:122313/1‑149 (MQ=255)

                                                                                                aaTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGctct                                                   <  3:301587/149‑1 (MQ=255)

                                                                                                 aTGCGCTTACTGATGCGGATTTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCGATATGACTGCGAGATTTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTg                                                  <  5:220628/149‑1 (MQ=255)

                                                                                                 aTGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTg                                                  >  5:146511/1‑149 (MQ=255)

                                                                                                  tGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCTATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTgg                                                 <  3:147827/149‑1 (MQ=255)

                                                                                                  tGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTgg                                                 <  4:308328/149‑1 (MQ=255)

                                                                                                    cgcttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGc                                                <  7:3285/142‑1 (MQ=255)

                                                                                                    cgcttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCAATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCg                                               <  4:223003/143‑1 (MQ=255)

                                                                                                    cgcttactgATGCGGAATTACGCCGTAAGGCCGCACATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCAGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGc                                                >  4:49240/7‑148 (MQ=255)

                                                                                                      cttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCg                                               >  1:258878/5‑147 (MQ=255)

                                                                                                      cttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGaa                                             <  5:54368/145‑1 (MQ=255)

                                                                                                      cttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGaa                                             <  3:86087/145‑1 (MQ=255)

                                                                                                         actgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTcc                                                                                          <  7:283003/100‑1 (MQ=255)

                                                                                                            gatgcgcAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCTATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGCCTCTGGCGAAGATCg                                        >  2:331808/4‑148 (MQ=255)

                                                                                                            gatgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATGTGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCg                                        >  1:390860/4‑148 (MQ=255)

                                                                                                             atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCg                                        >  8:296757/3‑147 (MQ=255)

                                                                                                             atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGc                                       <  5:250567/146‑1 (MQ=255)

                                                                                                             atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGc                                       <  6:172388/146‑1 (MQ=255)

                                                                                                             atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCa                                      <  7:258116/147‑1 (MQ=255)

                                                                                                             atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCa                                      <  2:277136/147‑1 (MQ=255)

                                                                                                             atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCa                                      <  6:369222/147‑1 (MQ=255)

                                                                                                             atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCa                                      <  6:114815/147‑1 (MQ=255)

                                                                                                               gcgGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGCTCGCAg                                     >  8:262616/1‑148 (MQ=255)

                                                                                                               gcgGAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAg                                     >  7:16338/1‑148 (MQ=255)

                                                                                                                 ggAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGaaa                                  >  7:358976/1‑149 (MQ=255)

                                                                                                                 ggAATTACGCCGTAAGGCCGCAGATGAGC‑T‑TGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGaaa                                  <  2:355143/149‑1 (MQ=255)

|||

CGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑T‑GACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGT  >  NC_000913/566388‑566681

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