BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	566,535	2 bp→CA	100%	pseudogene (160‑161/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,535	0	T	C	97.6%	5893.8 / 46.6	2122	pseudogene (160/312 nt)	peaD	pseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand): ref base T (28/25); new base C (1101/946); total (1143/979)
Fisher's exact test for biased strand distribution p-value = 8.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
*	NC_000913	566,536	0	G	A	97.6%	6005.3 / 42.0	2121	pseudogene (161/312 nt)	peaD	pseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand): ref base G (30/23); new base A (1110/955); total (1142/979)
Fisher's exact test for biased strand distribution p-value = 7.81e-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 2122 total aligned reads displayed.

GGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGT  >  NC_000913/566389‑566681

||

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

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

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

                      tGCTGAGTCTTATCCGTGGAAATCAAACGCGCACTACTGGCGGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCgtg                                                                                                                            <  3:352336/149‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                     gTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCt                                                                                                              >  5:363952/1‑148 (MQ=255)

                                     gTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCt                                                                                                              >  1:336432/1‑148 (MQ=255)

                                     gTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTg                                                                                                             >  6:114090/1‑149 (MQ=255)

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

                                        gAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAGTTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAAc                                                                                                          <  8:207586/149‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                    gcACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAaacaac                                                                                              <  6:340485/149‑1 (MQ=255)

                                                       ctACTGGCTTGTTACCAACCTGTATCAGCACCTGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGCGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTc                                                                                           <  6:235820/148‑1 (MQ=255)

                                                          ctggctggTTACCAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCt                                                                                         <  6:87565/148‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                                     ccAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGAtcc                                                                                                                <  5:248820/114‑3 (MQ=255)

                                                                     ccAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg                                                                              <  4:84379/148‑1 (MQ=255)

                                                                     ccAACCGGTATCAGAACCTGCGGGCCAATGCGCGTACTGATGCGGAATTACGCCGTAAGGCCGCAGAAGAGC‑TTGTCCATATGACTGCGATAATTAACCGTGGGGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGt                                                                             <  4:306797/149‑1 (MQ=255)

                                                                      cAACCTGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGCCTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCgg                                                                              >  2:330535/1‑147 (MQ=255)

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

                                                                       aaCCTGTATCAGAACATGCGGGCCAATGCTCTTACTGATGCTGAAATACGCCTTAAGGCAGCAGATGATC‑TTGTCCATATGAATGATAGAATTAACCTTGGTTAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAg                                                                           <  5:333511/149‑1 (MQ=255)

                                                                         cctGTATCAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGAc                                                                         <  7:231735/147‑1 (MQ=255)

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

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

                                                                               tcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGCCC‑TTGTCCATATGACTGCGAGCATTAACCGTGGGGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTa                                                                    >  3:42706/2‑148 (MQ=255)

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

                                                                               tcAGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTa                                                                    >  5:113166/2‑148 (MQ=255)

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

                                                                                 aGAACATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGGGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTaaa                                                                  <  5:53731/148‑1 (MQ=255)

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

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

                                                                                    aCATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGt                                                              <  4:174683/149‑1 (MQ=255)

                                                                                    aCATGCGGGCCAATGCGCTTACTGATGCGGAATTACGCCGTAAGACCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGt                                                              <  3:234889/149‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                   cgcttactgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGc                                                <  5:21869/142‑1 (MQ=255)

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

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

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

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

                                                                                                     cttactgATGAGGAATTACGCAGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGAGGTGAGCCGATCCCTCAACAAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATAGTGAACACGCTCAGGCg                                               >  4:280104/5‑147 (MQ=255)

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

                                                                                                        actgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGa                                           <  7:331289/147‑1 (MQ=255)

                                                                                                          tgATGCGGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCg                                                               >  2:237517/1‑126 (MQ=255)

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

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

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

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

                                                                                                            atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCa                                      <  3:122300/147‑1 (MQ=255)

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

                                                                                                            atgcgGAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCa                                      <  6:172373/147‑1 (MQ=255)

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

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

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

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

                                                                                                                ggAATTACGCCGTAAGGCCGCAGATGAGC‑TTGTCCATATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGa                                    >  6:147889/1‑147 (MQ=255)

||

GGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGTTACCAACTTGTACCAGAACATGCGGGCCAATGCGCTGGCTGACGCGGAATTACGGCGCAAGGCTGCCGATGAAC‑TGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGT  >  NC_000913/566389‑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: `‑`